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BY i 
J. L. SYNGE{ AND W. F. CAHILL 
National Bureau of Standards 


Summary. SEACt is used to give an approximate solution to the torsion problem 
for a hollow square, with accurate determination of lower and upper bounds for torsional 
rigidity. 

1. The torsion problem for a multiply connected cross section. To solve the torsion 
problem for a beam of multiply connected cross section, we seek a function ¥(x, y) which 
is harmonic in the domain A occupied by the material, and such that - 


yY = }r’ on the outer boundary B, , (1.1) 
vy = 3r° + C, on the inner boundaries B, , 


the constants C, being chosen so that for each inner boundary we have 


[ oF an = 0. (1.2) 


By on 


Here r° = x + y’, where x, y are rectangular Cartesian coordinates in the plane of 
the cross section. 


The components of shearing stress are 


a 


‘ éO _ am 
ay (y saa ir”), — Ma ax (y ar), (1.3) 


pa 


where u is the rigidity and a the twist per unit length. The lines of stress are given by 
the equation 
y — 47° = const., (1.4) 
this family of curves including, of course, the curves B, , B, . 
The warping of the section is given by ¢(x, y), which is the harmonic function con- 


jugate to y. 
If we denote the torque by wal, then 


ror f (CR) +(e) “ 
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when J is the polar moment of inertia of the cross section, 
I = | (x? + y*) dA. (1.6) 


2. The torsion of a hollow square. The method used below for solving approxi- 
mately the torsion probl m for a hollow square cross section is essentially the hypercircle 
method, originally developed by Prager and Synge’ for more general problems in 
elasticity. 

Details of the method will be found in a forthcoming book,’ togethcr with numerical 
results obtained by a desk calculator. The present paper is concerned with the application 
of SEAC to these calculations, and only essential formulas will be given to make the 
computations intelligibl 

The method is applicable in principle to any cross section, and is particularly suited 
to a cross section which is a hollow rectangle or a hollow square. The only case treated 
here in detail is the case of a hollow square bounded externally by a square J, of side s 
and inte rnally by & square B, of side 58, as shown in Fig. 1. 

The square is triangulated as shown, a good approximation demanding a fine tri- 
angulation. We speak of approximation of order n when the outer side is divided into 
n equal parts. On account of symmetry, only an octant is needed in the calculations; 
Fig. 1 shows the triangulated octant MNZY for n = 16. 
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Fic. 1. Triangulation for the approximation n = 16. 


Each junction in the triangulation is a “station.” There are two types of station, 


if} 


type P at the centers of the small squares and type Q at the corners of these small 


squares. 





1W. Prager and J. L. Synge, Quart. Appl. Math. 5 (1947), 241-269. 
2J. L. Synge, The hypercircle in mathematical physics, Cambridge University Press. 
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To each station we attach a “weight’”’ w. These weights uniquely determine a function 
f(x, y) in the domain A of the cross section, by the rule that f(z, y) = w at each station 


and f(z, y) is a linear function of x and y in each triangle. 
The general plan of approximation is to find two sets of weights, say w’ and w”, 
such that if f’(x, y) and f’’(x, y) are the corresponding functions, defined as above, then 


the following approximations hold: 
“1 / é 1.2 
(x7, y) ~ — a 
J , ¥) vy » ? (2.1) 
f(x,y) ~ ¢. 
This plan is slightly modified to permit dimensionless calculations, “reduced weights” 


b’, b” being introduced. 
3. Approximation for lines of stress, and lower bound for torsional rigidity. To 
tr’ we have to solve the following difference equations: 


approximate y Ay 
b’(P) = 1 > b’(Q) — dP) b(R) +- 1, 


b'(Q) = 1 > bP) — dQ) b’(R) + 2, 3.1) 
b(R) = —5 DWP) +0144, 
ole M NV 12 
e 


b’(Q) = 0 for Q on By or B, . 


‘(P), b'(Q) are reduced weights at stations P, Q. The interpretation of the first 
3.1) is as follows: The }> is a sum over the four stations Q which are neighbors 


of P. d(P | if P is on MN, and otherwise d(P) = 0. b’(R) is an additional reduced 


weight, i uced to allow for multiple-connectivity. The interpretation of the second 
line of (3.1) is similar. Note that if Q is on B, or B, it lacks some of its four neighbors; 
in that case fictitious neighbors are supplied with zero weights. In the third line of (3.1) 
os is summation along the diagonal AZN of Fig. 1 


The weights b’(P), b’(Q) are symmetric with respect to the lines of symmetry of the 


The approximation to y — 4r° is given by 
mn Se ae Y= pee a 
Y-—arw~soF'(«,y +5— d(C, y), (3.2) 
’ 6 n 3 Nn 
where F is a function which is linear in each triangle and takes, at the stations 
P, Q, the values b’(P), b’(Q), while G(x, y) is a function which is linear in each of the 


four triangles in which the outer square is divided by its diagonals, with G = 1 at the 


center of the squares, andG = Oon B,. 
If we denote the torque by pal, then a lower bound for torsional rigidity is given by 


4 ] /D\ < , _! twp , 
$2971 bP) +2 DVO) +3, O'R), (3.3) 


these summations extending over the whole domain A, and not merely the octant 


MNZY. 
4. Approximation for warping function, and upper bound for torsional rigidity. ‘To 


approximate ¢ we have to solve 
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7 


A 


b’"(P) , : b’’(Q), (4.1) 


| 


A 


b’"(Q) 1S” b/(P) + m(Q). 


Here each summation is over four neighbors. If neighbors are lacking, i.e. when Q is 
on B, or B, , fictitious neighbors are to be added by reflection in By or B, , the images 


having the same weights as the actual neighbors. The number m(Q) is given by 


m(@ 0 if Q is inside the domain A, 
m(Q 2x/s if Q is on By (the outer boundary), 

1.2) 
m(Q —2r/s if Q is on B, (the inner boundary). 


Here z is the abscissa of Q, the origin of coordinates being at the center of the squares. 
In this case the weights are skew-symmetric; therefore 


b’’(P) 0, b”(Q) 0, 1.3) 


if P or Q is on a line of symmetry of the section. 


The approximation to the warping function ¢ is given by 


@ 7 P(e, y)., (4.4) 


-—y, ° ° . @ ver ir 
where F’’(x, 7) is linear in each triangle and has the values b’’(P), b’’(Q) at the stations 


P,Q. 


An upper bound for the torsional rigidity is given by 


AS 1 o 


r . 1 [5n' . . : 
. 29 - bs m( bh’ (Q) _— p> m Q)b’ (Q) : (4.5) 


the summations being along ZN and YM in Fig. 1. 
Q,, Q, 2 Qs Qi2+1 


Pi, Pre Pis : : : Pio 


Q,, Q,. Q, 3 Q, 


x 






































QB) Qgri f+ 





Fic. 2. Enumeration of stations for the approximation (a multiple of 4 


5. Plan for SEAC calculations. For the approximation of order n, the stations 
were numbered as indicated in Fig. 2." Then, using the letters P, Q, R for the reduced 
.dditional complication, to the case where the ratio of 


1The procedure may be applied with little 
the sides s, s’ of the outer and inner squares is not 2:1. If it is possible to triangulate with n squares 


along s and n’ squares along s’, then we make the following notational changes in Fig. 2: 


Qinss ae > Qin n’)/2+1,1 
Qin 4) +1, (n/4) +1 — Qin n’)/2+1,(n'/2) +1 
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weights b’, we have, as in §3, 
O(a, y) Q(y, x) Q(—z, y) = Qa, -—y) = QA—2z, —-y), 
P(x, y) P(y, x) P(—2z, y) = P(x, —y) = P(—x, —7), 


7m 


R = 2 > (P+ Q) + 13° 
P 1(Q ] f Qi. l t @) +1,7 t Qi 41.541) + ] = dU(P,)R, 
Q) (iP t+ Py + Pipes + Pi) +2 — dQiaR, 
U(X,,) {I if X;,; sits on diagonal, 
lo otherwise 
| 1; () n/4)+1,) 0. 


The lower bound for I is given by 


r) SO Des eet } 
(= <i wat a tr 2Q) + 36n od be 


B ALL 


Using P, QO for the reduced weights b’’, we have as in §4 
« z ’ ’ 


P(x, y P(—2, y) = —P(a, —y) P(—x, —y) —P(y, x). 


lor interior points, 


on the diagonal 


P,; QQ; 0; 
on the outer boundary B, , 
2(7 — 1) 
Oy = MPiia + Py) + 
n 

ii 1, 0 

on the inner boundary B, , 
2(7 1) 
@ 4 1 L(P., 4 l a FP. 4.3) —_ ; 


Q,, 4)+1,1 Qi, 4)+1,(n/4) +1 0; 


on the vertical axis, 


Q, = 0. 


pe 
Ox, y Q(—x, y) —(Q(ax, —y) Q(-—z, —y) =,.—Qy, 2), 
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5.10) 


a | 
— 
w 


(5.14) 


(5.15) 


(5.16) 
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‘is given by 


— 

3 
a | 
5 


rT \ 5 1 (n/2) +1 (n o) +1 a es 
i) =a5—73) Le QumQu) + D> Qerer.aM@Qwrrsrad?- (5.17) 
\S /y Of U k=1 k=1 ‘ 


The iteration scheme for the lower bound started with all P;; = Q; 0 and the 
following cycle describes one complete iteration: 
1) Set 2 l 
2) Compute all P in row 7 by formula (5.4) 
3) Have we done last row of P? 
No, go to step 4 


Yes, go to step 7. 
t) Replace 2 b: 
5) Compute all Q row z by formula (5.5 | 
6 Ret Irn to ste | 
7) Compute Rk tula (5.3 
8) Ret to 
After eve ry T iterations [where 7’ approximate ly 20) is chosen at the start of a 


problem] the appl xKiImations for some cho en P. for some chosen o, for R and also the 


t : ry Rey ae 
& od Bis Ls BB pa eee 
“ i : — 2 J 2 
‘ VEN ec ) i os : 
‘ 1 ~ i | AL SS kl Tee a oF + 1 
ot | spt TT ae 
Bead 30] -65 Fos F147 $196} 
| Leo Lrof veil coe) 


Ec oe ] fies 4. - r..| 
— } ees | 4 4 4 i as ded 
07] v6 [vs Tras 1 Lvs | e}-rat -s09 
| ise] | -<e] ] 
= - toy ‘nee | +——i—__i - 

| | | | 
— . ail 7 i t,t. - 

| | 
Ko Re GiGi RO tbo oe eo ea 
ms OR A ca OE A A re 
Meda) Melis) lubed Sick st hed Bt Ha | (Bo BUR. { 


mo Bcd Be 8 BPA Bcd) MLD | As pil | od 
Wet 


Bed Bess hon! Bit Beast Bia 2 ] L ; 33 
oT LL: 4. 1 es 
i 


. T a E T T * T T | 
CoTeToToToTotefofofofo] oo] 
48. The numbers are solutions of Eq. (3.1), 


TaBLe I, ¢ omputatlol of lower bound in approximation n= 
é i 
equivalently (5.1) 5.7), divided by 10 and rounded off. 


R = 30.13171 


— 





#8) 
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sum of the absolute values of the residuals are printed. When this sum becomes sufficiently 
small, all values are printed. 
The iteration scheme for the upper bound is similar to that of the lower bound. 


6. Numerical results. The equations (3.1), or equivalently (5.1)-(5.7), were solved 
as indicated in §5, in the approximations n = 8, n = 16, n = 32, n = 48. This last was 
the highest approximation feasible on SEAC, without resorting to tape storage. 

The solutions for n {8 rounded off to a maximum of three significant figures, are 
exhibited in Table I. To within a constant factor, these numbers are the values of the 
function F’(x, y) fef. (3.2)] at the stations. Addition of the last term in (3.2) to the 
numbers in Table I yields the numbers in Table II, which therefore represent, to within 
a constant factor, the stress function (¥ — 47°). A plot of lines of stress could easily be 


pre pared from ‘Table II. 
Similarly, the equations (4.1), or equivalently (5.9)-(5.16), were solved for the 
me approximations. The solutions for n 18, rounded off to a maximum of three 
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Taste II. Computation of lower bound in approximation n = 48. The numbers represent approximate 
values of 


[((6 X< 48)?/s*] (y — Mr?) 


significant figures, are exhibited in Table III. These numbers represent, to within a 
constant factor, the values of the warping function ¢, and from these numbers a plot 
of level lines of warping could easily be prepared. 
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TABLE III. Computation of upper bound in approximation n 18. The numbers are solutions of I¢q.(4.1), 
or equivalent 5.9 ».16), multiplied by LOO und rounded off, they represent approxim ite values of 
100 » t X 48)/s*] ¢. 

The bounds « torsional rigidity, and the number of iterations required to give 
solutions to nine significant figures (starting from zero values) in the several approxi- 
mations, are shown in Table IV. 

Tasie I\ 
Lower and upper bounds for I'/s* 
Number of 
Iterations 
n Lower Bound Upper Bound (U.B. L.B.) L.B. U.B. 
8 1263 0626 .1320 1773 0057 1147 50 60 
16 1282 4044 1300 3396 0017 9352 180 270 
2 1288 4044 1294 2496 0005 8452 680 1070 

18 1289 7704 .1292 8697 .0003 0993 1400 2200 

In the approximation n, the number of simultaneous equations to be solved is 3n°/16 
in the case of the lower bound, and one less in the case of the upper bound. Thus, for 


n 18, there were 432 and 431 equations respectively. 
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By KURT EISEMANN 
(IBM Corporation, New York) 


Objective. The most common problem in numerical computation is the solution 
of simultaneous linear equations. The most frequent obstacle is ill-conditioning of the 
matrix of coefficients, resulting in substantial loss of accuracy. The present paper 
describes a method of overcoming this difficulty as soon as encountered—with a mini- 
mum of computational effort, preservation of reasonable accuracy, and emergence of a 
well-conditioned system. 

The present method applies only to cases in which the equations describe a well- 
behaved system but are poorly specified in the sense that one or a few of the given equa- 
tions are almost a linear combination of other equations. The method does not apply to 
eases like the Hilbert matrices, in which every row is nearly a linear combination of other 
rows. 

To avoid encumbrance by irrelevant details and to fix ideas, we shall illustrate by 
a matrix A of order n 5 throughout. 

The mathematical key. With high-speed computers, the preferred method of 
solution is some variant of the Gauss elimination procedure. After k stages, A has been 
reduced to a form R**’ which, for k 2, is illustrated by 


Ll fic Tis Tia Mis 





In general we can write 
(“ ee where \R,, = upper triangular k X k, 
0 Ss (iS=S” =(n—-K X(n— Bb. 


Retracing backwards the steps of elimination, matrix R‘" is retrieved from R” as 


7 ! 
0 I, | 

KR =T°R®, where T? =|0 by 1 , 
0 ls O 1 


LO ls O O 1 


the elements /,, of column 2 being precisely the elements s‘}’ in the leftmost column of 


*Received Aug. 10, 1956. 
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S“. Continuing in this manner, matrix A is expressed in the form 





; Au RP a FT 6 FR a LOR”, 
where L i --- T™ is illustrated, for k = 2, by 
(1 ) [lu 
lin 1 110 les | . & 
Lb? = 77? = 1h, 01 1}0 be 1 J = |ts bes 1 
lk, OO 1 }}0 de O 1 ln ly, O 1 
lI, 000 11/0 lL, 0 0 1 . & © @ 1 
Similarly, in general, 
7 poy pe ix, © where L,, = lower triangular k X k 
i. J IT=(n-—k)X(n—-k) 


As seen above, the elements /;, are furnished directly by the elimination procedure: 
the leftmost column of S*”’ always becomes, at the next stage, the new column k + 1 
apr”. 

It is noteworthy that matrices L“’, R“’ and the factorization A = L“”’R“™”’ represent 
a generalization of the usual formula A = LR = L™ R™ in which L and FP are postulated 


to be triangular. 


We shall subsequently also need P inverse of L*”’, which can be written 
' 'P,, O , .. Py, = lower triangular k X k, 
I ) where 
ae Je IT =(n—k) X(n— hk). 


In the sequel we shall omit superscript k, it being understood that all matrices refer 
to’stage k. Instead, superscripts will henceforth denote row numbers. 

The source of trouble. Now ill-conditioning reveals itself by the fact that at some 
stage in the elimination procedure an entire row of S takes on very small values. 

Consider the case where all elements of row S‘ are very small. Elimination will then 
lead, at the present or some later stage, to division of this row by one of its very small 
current elements in order to provide a coefficient of value 1. This is equivalent to multipli- 
cation of row 7 by a very large number. Accumulated round-off errors are thereby greatly 
magnified. It is the high relative errors infecting this row that introduce the sudden 
substantial drop in accuracy, subsequently propagating to the entire solution vector. 
The situation could be saved if it were possible to obtain row S‘, which is the result of 
many computational steps, free of the high accumulated round-off errors. We shall show 
how to accomplish this with a comparatively small amount of labor. 

The remedy. 5S‘ is equivalent to the result of subtracting a linear combination of 
rows A! --- A* from A*. We would first like to determine approximately the constants 
of this linear combination, indicating in what way S° is obtainable directly from the given 


original matrix A. 
: -1 
From A = LR (stage k understood) we obtain R = L~'A = PA. We can therefore 


obtain S* (¢ > k) from 
S'~R' =P'A= (Po ,)A = Dir A’ 4- +e +} DA’ + A’. 
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Since the elements of the left-hand side are very small, this relation indicates that 


. . . . 1 k 
Dit *** Di, are the negatives of the constants of linear combination of rows A’ --- A 


approximately furnishing A’. 
We therefore need the elements of P* for i > k, i.e. of P;, . They are found by con- 
sidering 


P, O\fL, 0\ _ ff 


PL i. LO 
Pes I I 02, I 


and equating the expressions for block matrix 0,, . For row i ( > k), this decomposition 
yields P3,L,, + Le; 0, i.e., 
Dial + Piolo, +--+ + Dinter as —I;, ’ 
Pialos -+- ed Pieler _ —I;. ’ 


The coefficients of this system are already available from preceding reduction steps. 
However, they bear accumulated round-off errors. Solution of this system would there- 
fore only yield values p{; which approximate the true p,;; . Fortunately, that is all we 
need. 
The preceding system, which is only of k-th order and triangular, is therefore solved 
recursively for the approximate solution p/, --- pi, . In a computer, the program for 
the “back solution’’ of the elimination procedure, already available within the machine, 
is effectively drawn upon for this purpose. 

As shown above, the exact p;; , applied to rows A’, would yield the row S* of very 
small elements, caused by cancellations in the most significant digits of each element. 
Instead, we take a linear combination with values p/; and replace the original matrix 


TOW : by 
A" = ph A'+--- + phA*+ A’. 


Similar cancellations of significant digits will again take place, leading to a row A” 
with small elements. This time, however, we shall take care to calculate with double- 
precision, furnishing small elements of A’* that are exact. 

Now we multiply through by a large number M to convert all elements to O(1). 
Since the small elements were exact, no error is introduced or magnified by this step. 
The usual slump in accuracy has thus been completely prevented, maintaining the 
latter at its normal level. 

The first k elements of the resulting row are now removed, and the row is thereby 
reduced to the form of R*, by application of the elimination algorithm, utilizing the 
available rows R' --- R*. It is in this reduction that the bulk of accumulation of round- 
off errors occurs. Now, however, the final R* will have elements of O(1), so that the 
effect of rounding errors will be no more appreciable than for all other normal rows. 

Throughout k stages of reduction, row 7 (> k) has no effect on values in all other rows 
of R. Hence, apart from row 7, the entire matrix RF is retained as is. 

We may look upon the reconditioning of matrix A in yet another light: When two 
matrices have their elements scaled such that corresponding row vectors have equal 
length (say all length 1), a comparison of their determinants provides an effective 
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i] 
i) 
ra) 


comparison of ‘condition’. Consider this done for the original matrix A (containing 
A‘) and for the modified A’ (containing A’’ times the large number M). After complete 


reduction, we have | A |L|-| R | [ [t_, J;; . In the original system, a very small 
l;; gave rise to a small | A |, indicating ill-condition. In the modified system, the element 


l, in the identical position has been brought up to O(1), carrying | A’ | to a comparatively 
high magnitude, which indicates favorable condition. 

Summary. We summarize the auxiliary computations on encountering ill-condition- 
ing: (1) recursive solution of a triangular linear system of order k; (2) double-precision 
recalculation of row A”, followed by a “left shift’’; (3) elimination of the first k elements 
of row 2 

It is noteworthy that removal of ill-conditioning requires the double-precision 
recalculation of only a single row and that intermediary results for all other rows are 
retained unaltered, thereby salvaging the bulk of computational effort. 

Numerical example. We now illustrate the power of the present method by a 
numerical example. To simplify illustration, we shall desist from “positioning for size”’ 
for optimum divisors, but shall treat rows and columns in their original order. 

Given the system 


1.3827, — 4.732, + 5.392, — 2.842, + 3.972; 1.33, 
5.682 6.252, + 1.40273 + 7.027, + 4.502; —6.04, 
1.9382, + 1.38427, — 2.167, + 3.8la, — 2.622, = —1.75, 
2.852, + 3.097. + 4.412, + 2.3627, + 3.142; 2.34, 
1.3227, + 0.207, + 4.692, — 1.6327, — 5.392; = 4.50. 


Suppose we stipulate that whenever, in any row, the element of maximum magnitude 
becomes less than 0.20, the entire row has fallen below an acceptable threshold of 


accuracy. 
Applying straightforward elimination and rounding all figures to two decimal places, 
we obtain, after k 3 stages, the system R“’x b” ie. A = LR’, with 


the numerical value 
| 1 -3.58 1.08 —2.15 3.01 1.01 
| ] 1.55 1.37 -—0.89 —0.84 
- ae 8 | 1 —1.21 —0.39 1.17], 
| 6.47 11.61 —5.03 





—0.08 —0.01 0.02) 


L®™ = 11.938 8.25 2.76 


2.85 13.29 13.38 1 








(14.32 15.67 11.35 0 1 
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The smallness of elements in row 7 = 5 now reveals the hitherto hidden fact that 
row 5 of A is almost a linear combination of the first k = 3 matrix rows. Writing L” = 
transpose of the upper left-hand 3 X 3 submatrix L,, of L“™, we therefore solve the 


triangular system Lp’ + /' = 0, i.e. 


1 O2)s51 -f 5.68ps52 f. 1 93p53 =a 4 oa 
] 4 O8ps2 +} 8.25ps3 ae 15.67 
2.76ps3 =a 1 1 ao 


still always rounding to two decimal places, yielding the approximate solution pj, = 
—2.86, pis 1.30, pf, = —4.11. These values are utilized in replacing row 5 of the 
original augmented matrix (A, b) by the row 
A” = pi, A’ + pi,A" + pizA°* + A’ 
calculated to double-precision, i.e. the fifth equation is replaced by 
-0.00352, + 0.09542, — 0.02782, — 0.0407zx, — 0.12602; = 0.0367 
We now multiply through by a comparatively large number, say M = 100, and 


reduce the resultant first three terms to zero by subtracting multiples of rows 1, 2 and 
3 of the matrix (R™, b™). At the end of stage 3, row 5 of the augmented matrix (R“, 


b*’’) now becomes 
(0, 0, 0, — 2.26, 0.31; — 2.47), 
in which the numerical value of the element of largest magnitude is now far above the 


permissible minimum threshold of 0.20. 

Completion of numerical solution yields final values for the unknowns which are 
compared in the table below. Deviations 6; from the true values x; as well as max | 4; | 
and >> | 4, | are likewise tabulated. The improvement in accuracy is indeed striking. 


Comparison of numerical values 

















Quantity Correct Usual Present | Quantity Usual Present 

values method method method method 

Ly —2 —0.30 —1.96 by —1.70 —0.04 

L2 0 0.39 0.00 5a —0.39 0.00 

Z3 2 0.77 1.95 53 1.23 0.05 

x4 1 —0.23 0.96 54 1.23 0.04 

Ls —1 —0.31 —0.97 55 —0.69 —0.03 
max |5;| 0 1.70 0.05 
> 14; 0 5.24 0.16 
p 0 87 3 








Relative percentage deviations 100 | 6; |/| 2; | for individual components become 
infinite or unduly large whenever | x; | = 0 or very small. In order to summarize the 
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overall relative inaccuracy by a single number, avoiding at the same time any instability 


in the neighborhood of small | x; |, we define as an effective measure of overall inaccuracy 


p = overall percentage error = 100 >> | 6; |/>> | 2; | 
Numerical values of overall measure p, shown in the above table, likewise compare very 


avorably. 
Closeness of solution. 
equations, we find that the residual vectors for the two sets of unknowns are both very 


Substituting the numerical solution into the given set of 
small and have elements of comparable magnitude. How then can we say that our 
second set of xz; is any “‘better’’ than the first? 

Let x = exact solution vector, 2’ computed solution vector, 6 
vector for the unknowns, r b — Az’ residual vector for the computed solution. 
Substitution for 2’ in the expression for r yields r = Aé. 

In each case we obtain a small vector r, but this cannot serve as a valid criterion 


. , a . 
~~ a SS error 


, 


for closeness of solution. Our prime objective is the reduction not of r but of 6 = x — 2’ = 


=-!1 ° ° . » ° 
A~’r, and this objective has been successfully achieved. 


231 


TWO DIMENSIONAL COMPRESSIBLE SHEAR FLOW IN A CHANNEL* 


BY 
JAMES D. MURRAY** 


King’s College, University of Durham, Newcastle-on-T yne, England 


1. Introduction. As far as the author is aware, no theoretical or numerical solution 
has been obtained for the two-dimensional, inviscid, adiabatic rotational flow of a 
compressible gas through a divergent channel such as shown in Fig. 1. Kramer and 
Stanitz [2] have obtained a relaxation solution for the incompressible shear flow in a 
90° elbow. Mitchell [3] has derived a method for evaluating the non-isentropic rotational 
field downstream of the bow shock wave formed when a supersonic stream impinges 


on a blunt nosed obstacle. 


: : 











X=- 00 








Fig. 


In the present paper, it is intended to illustrate the principal physical features of 
the completely sub- and super-sonic adiabatic rotational flow through a channel, with 
particular reference to the effect of the divergence of the channel on the vorticity and 
the stream lines. 

2. Fundamental equations and boundary conditions. In the two-dimensional 
steady motion of an inviscid, adiabatic, rotational compressible gas through a divergent 
channel (see Fig. 1), the stream function y is defined in cartesian coordinates by the 


equations 


_ ,19¥ — —,79¥ 
q:z all? oy ’ Vy ~ et p Ox ’ 


where g, and q, are the cartesian velocity components, and p is the density. The vorticity 


w is defined by 


— 9% _ 9% 
~ Ox oy ’ (1) 


*Received March 30, 1956; revised manuscript received June 11, 1956. 
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which, from Vaszonyi [4], gives the equation for y as 


4 9 fe) 6 
BM (_- ov) + J Ge ¥) —p * = 0, (2) 


Ox Ox oy oy 


where h, is the stagnation enthalpy. 
Bernoulli’s equation for an adiabatic gas may be written as 


vy — 1)"'p/p + 3¢° =h,, (3) 


where p is the pressure, q the velocity, y the ratio of specific heats, and pp” constant. 


In an adiabatic rotational flow, h, is necessarily a function of y. 

As to the nature of the boundary conditions necessary for a unique solution to these 
equations in the general case, there is still some doubt. In the particular problem of the 
flow through a divergent channel with straight parallel sides far upstream and down- 
stream of the constriction, the boundary conditions are more straightforward. The 
physical assumption is that at infinity upstream and downstream of the constriction 
the flow is parallel to the channel walls. Thus, taking the z-axis parallel and along the 
lower wall in the narrow section of the channel, the condition is df/dx = 0 at x 0 , 
©, the flow is parallel, the pressure, and 


the constriction being x 0. Since at x 
(or the 


consequently the density, is constant across the channel. Accordingly, if h 
vorticity) is given as a function of y, y is obtained from Eq. (2), and is thus known 
round a closed boundary, enclosing the field of flow, which in the elliptic case (subsonic 
0, 


flow) defines a unique solution over the complete field. Alternatively, if dy/dx 
+o, a solution is again 


and y is given upstream at x —o, and Oy /Ox GQ at x 


defined. 
Letting the suffix zero denote quantities at x , It is convenient at this stage 


to introduce non-dimensionalising quantities py) , po , [1 the channel width, and ¢, the 
f 


speed of sound, where « YPo/po . Thus, with only slight change in notation, y, p, p, 
q, w, h, , x and y, now stand for the non-dimensional quantities ¥(pocol1)~, p/Po , p/po ; 


g/Co _ wH/c, , h,/c, , x/H and y/H. In the following, it must be remembered that all 


quantities are non-dimensional. 


In the specific problem, let the initial vorticity w, be a constant, which from Eq. 


(1) gives 


(q- ,= M,-— Woy, (4) 


where J/, is the initial Mach number of the flow on the lower channel wall. The density 


is initially constant and equal to unity, which gives 
Y = Moy — fuoy’. (5) 
From Eqs. (3), (4) and (5) the stagnation enthalpy h, is given as a function of y by 
h, = (1/y — 1 + 4M) — wo. (6) 


The function h,(y) given by Eq. (6) remains unaltered throughout the flow, since h, 
is constant along a streamline. 
= -++o, where the channel width is 2H (the general 


In the parallel flow region at x = + 
case rH is treated in a similar manner), d¥/dx = 0. Denoting conditions in this region 


l 
by the suffix unity, Eq. (2) reduces to 
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a ( 7 v1) wae 
oy p oy Pp 0% 

which on integration gives 
Vi = A+ By - Spiwoy”, 


where A and B are constants. On the lower channel wall at y = —}3, ¥, = 0, and on the 
upper wall at y = 3/2, ¥, = My — 4» [from Eq. (5)], which results in 


vi = (2M, — wo + Bwopi) + 32M — wo + opi)y — Fwopiy’- (7) 
Thus, with the exception of p, , all quantities are known analytically in the parallel 
flow region at x +o, Equation (3) must be satisfied at all points in the field, which 


in dimensionless form gives 


; adi site 1/ _, dy,\’ e 
y= 2 p? ' + A( pi 242) = (y ~ 1) ' + 1M3 in wor ’ 


where y, is given by Eq. (7). This last equation reduces to 
Wo Pt + 2(y = I) "or + 71 [Mw =~ bers _ Ms Fa ' 2(y a 1)~"] (8) 
+ (4M, — w/4)* = 0. 


By Descarte’s theorem, there cannot be more than two real positive roots of Eq. (8). 
It should be noted that Eq. (8) is not the general equation, but is that for a divergent 
two-dimensional channel whose area ratio is two. The general equation for any area 
ratio can be obtained quickly using the above method. 

In Eq. (3), the stagnation enthalpy h, is constant along a stream line (although it 
varies from one stream line to another), and so for any one such stream line 

(y — 1) + 390 = (Y — 1) "ep * + 3a 
From this last equation, it is seen that if q, < qo, 1: > 1, and if q, > qo, o. < 1. Thus, 
the fact that p, > 1 for subsonic flow, and p, < 1 for supersonic flow is also the case 
when w, > 0. It is true for w = 0, which is the well known irrotational flow through a 
divergent channel. Accordingly, the appropriate solution of Eq. (8) is chosen in each 
case: these are now considered separately. 

3. Subsonic flow. In subsonic flow through the channel, the velocity at station 
unity is less than at station zero, and so p, is necessarily greater than unity. Thus, the 
requisite solution of Eq. (8) is that for p, > 1. With the restriction on w, , given subse- 
quently by Eq. (10), it can be quickly shown (by the method of signs) that the other 
real solution (> 0) of Eq. (8) is that in which p, < 1. From Eqs. (1) and (2), the final 


vorticity w, is given as 


dh, 
= "i en = PiWo ; (9) 


@) 


which gives w, > w, . That is, in the rotational flow of a compressible gas through a 
divergent channel, the vorticity will in general increase in value. Figure 2 illustrates how 
w, varies with w, , for several values of M, . There is a discontinuity in w, , which will 


occur when gp > 1 and q < 1 over some section of the channel near z = — ~. In such 
a case the problem is of mixed type, and there is still some doubt as to the boundary 
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conditions necessary for a solution to be defined. This difficulty does not arise when qo 
is everywhere greater than or less than unity. 

In any physical flow, the velocity in the parallel flow regions must necessarily be in 
the direction of increasing x, that is, g > 0 at all points near x = +o, Thus, w, must 
‘from Eq. (7). Together these 


satisfy Wo < M, from Kq. ( +). and Wy Ss 2M) (1 — tp; 
give the following restriction on the initial vorticity w 


wo < 2M,(1 + 49;) <M. (10) 


Equation (10) ensures that there is no ‘back’ flow region (that is ¢ > 0) in the field. The 


importance of this lies in the fact that if w < M, ensuring that q > O near x —@, 
the divergence of the channel may be such that the vorticity is increased sufficiently 
to make q < 0 at some points in the region near x +o: that is wy, violates the first 


inequality in Eq. (10). Physically this could perhaps result in breakaway of the flow, 
the cause of which has always been entirely attributed to viscosity effects. It should 
be pointed out that Eq. (10) is again for an area ratio of two. 

The effect of the divergence on the stream lines is marked, and is probably shown 


most profitably by the final deflection in the central stream line at z —o, that is 
(Wo) y= : (3M — w,/8). Substitution of this value for ¥ in Hq. (7), gives a quadratic 
equation, and the requisite value of y at 2 +o as the appropriate solution. The non- 
dimensional deflection 6 of this stream line at x +o is given by Eq. (11), where y 


is the solution to the quadratic with the negative sign with the root bracket. The de- 
flection 6 must tend to zero as w, tends to zero (irrotational flow), which necessitates 


the negative sign used in Eq. (11). Thus, the deflection is given by 
6 = y — 0.5 = (2a .p [((M, _ Fe — (Mo + w,/4 + 1ws pi —wM, - WP “], GD 


from which it is seen that 6 — 0 as w, — 0. The root bracket in Eq. (11) may be written 


as [(M, — 4a)” + wrp?(4p; — 1)]’, which is greater than (M, — },) when p, > 1, 
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which is the subsonic case. Thus, from Eq. (11), 6 < 0 when the flow is subsonic, for 
any w) and M, satisfying Eq. (10). The variation in 6 with initial vorticity for several 
values of the initial velocity M, , is illustrated in Fig. 3. One important feature of the 
flow arising from the figure is the fact that this deflection decreases with increase in the 
initial velocity for the same value of the vorticity. It should also be noted that the central 
stream line moves to a region of higher velocity at x = + © than the central value in 


this region. 
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In a problem of this type, which is elliptic, a simultaneous relaxation process can 
be carried out over the complete field of flow. The author has in fact evaluated such 
a problem, where the constriction was discontinuous, each wall having two right-angled 
corners. This relaxation solution was compared with a similar one for incompressible 
flow through the same channel, and it was found that in general the displacement effects 
were smaller in the compressible flow, but of the same sign. 

4. Supersonic flow. In the completely supersonic flow through the channel, assuming 
that no shock waves are present, the differential equations are intrinsically different. 


This case is a hyperbolic problem, amenable to solution by characteristics, and con- 
ditions need only be given initially at x — o, However, from physical considerations, 
it is assumed that the flow near x = + © is parallel to the walls of the channel as before. 

In this case the divergence of the channel increases the velocity, and the value of 
p, is the requisite solution of Iq. (8) satisfying 0 < p, < 1. Accordingly, it follows from 
Eq. (9) that a, < wp , and the restriction on w) for completely supersonic flow at all 


points is. irom Iq. (4 
wy < M, — 1. (12) 


Thus, the divergence decreases the vorticity, and the variation of w, with w» for several 
values of A/, is shown in Fig. 2. In this case qo > 1 at all points. 

The deflection in the central stream line is again given by Eq. (11), the negative 
sign being again used with the root bracket for the same reason as in Sec. 3. From Eq. 
(8), with w, and WW, (> 1) satisfying Eq. (12), it can be shown that p, is in fact less than 
0.5 [this can also be seen immediately from Fig. 2 and iq. (9)] and so the root exponent 
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in Eq. (11) becomes [(My — 4)” — wop?(1 — 4?)]?, which, since 1 > 4p? , is less than 
(My, — 30). Thus, Eq. (11) gives 6 > 0 for all w) and M, (> 1) satisfying Eq. (12) in 
the completely supersonic flow (that is, where p, < 1). The variation of 6 with w, and 
M, is shown in Fig. 3. The essential difference between this case and the subsonic case, 


is that the central stream line at x = — © is deflected towards a region of lower velocity 
at x = +o than the central value at this station. Further, the magnitude of the de- 


flection is very much smaller (in fact an order of magnitude different) in the supersonic 
case. One point of similarity is that the deflection decreases with increase in velocity. 

In supersonic shear flow the effect of the vorticity on the deflection of the downstream 
stream lines is an order of magnitude different from that in the subsonic flow. This 
phenomenon substantiates the experimental work of Davies [1], who found that the 
deflection in the free stream of the stagnation stream line when a pitot tube is placed 
in a supersonic shear flow is opposite in sign and of a very much smaller magnitude 
from that obtained by Young and Maas [5] in the low subsonic case. 

5. Discussion of results. By considering the differential equations (2) and (3), 
the variation in the vorticity due to the divergence in the channel is obtained, when 
the flow is, in one case, completely subsonic and, in the other, completely supersonic, 
It is shown that the difference between initial vorticities becomes smaller with increase 
in velocity in both cases, but in the subsonic case the final vorticity is greater than the 
initial value, the opposite being the case in the supersonic flow, as is shown in Fig. 2. 

The effect of the divergence of the channel on the stream lines is characterised by 
consideration of the deflection in the initial central stream line, as is illustrated in Fig. 
3. It is found that the deflection in the supersonic case is of an order of magnitude smaller 
and opposite in sign, from that found in the subsonic case. This result is in agreement 
with the state of affairs in the experimental work on supersonic flow past a pitot tube 
carried out by Davies. 

One final remark is that there is a possibility of ‘back’ flow, which could result in 
breakaway of the flow in the subsonic rotational flow through a channel of this type. 
Restrictions on the initial rotation and channel width necessary for no ‘back’ flow to 
be present are derived. 

The case of a divergent channel with general area ratio can be derived using the 


method described in this paper, with a proportional increase in algebraic labour 


SIBLIOGRAPHY 


1. F. V. Davies, Some effects of pitot size on the measurements of boundary layers in supersonic flow, R.A.E. 
T.N. Aero, 2179 (1952) 

2. J. J. Kramer and J. D. Stanitz, T’wo-dimensional shear flow in a 90° elbow, N.A.C.A. T.N., 2736 

(1952) 

3. A. R. Mitchell, Application of relaxation to the rotational field behind a bow shock wave, Quart. Mech. 
Appl. Math. 4, 371 (1951) 

. A. Vaszonyi, On rotational gas flows, Quart. Appl. Math. 3, (1945) 

5. A. D. Young and J. N. Maas, The behaviour of a pitot tube in a transverse total pressure gradient, A.R.C. 


Rep. 1770 (1936) 


ee 


237 


INTRINSIC FORM OF THE CHARACTERISTIC RELATIONS IN THE 
STEADY SUPERSONIC FLOW OF A COMPRESSIBLE FLUID* 


BY 
N. COBURN 
University of Michigan 


1. Introduction. ‘The purpose of this paper is to express the characteristic relations 
for the steady, supersonic, three-dimensional, rotational and irrotational motion of a 
polytropic gas in intrinsic form and to apply these relations to the study of a class of 
space Beltrami motions. This class of motions is characterized by the fact that one 
family of the characteristic surfaces consists of ~* parallel planes. 

To obtain intrinsic forms of the characteristic relations, two techniques must be 
used. First, the equations of motion, continuity, and energy must be expressed in terms 
of a set of characteristic variables. Secondly, these must be reduced to intrinsic form by 
expressing the partial derivatives of the components of the velocity vecter in terms of 
the directional derivatives of the magnitude of the velocity vector, g, and the sound 
speed, c, and the curvatures associated with the characteristic manifolds. 

The class of fluid motions to be studied is characterized by the fact that one family 
of * characteristic surfaces consists of planes, z2 = constant. If the additional assump- 
tion is made that the bicharacteristics and their orthogonal trajectories determine two 
families of cylindrical surfaces with generators parallel to the z-axis, then it is shown 
that these surfaces are right circular cylinders with oz as axis; the bicharacteristics are 
concentric circles; the stream lines are helices along which g (or /, the Mach number) 
is constant. 

2. The basic relation in terms of characteristic variables. Let x’, 7 = 1, 2, 3, denote 
a Cartesian orthogonal coordinate system in Euclidean three-space and let us denote 


partial derivatives by the symbolism 


In a Cartesian orthogonal coordinate system, covariant and contravariant quantities 
are equivalent. However, in order to use the Einstein summation convention of summing 
on repeated lower and upper indices, we shall use these two equivalent quantities. 
First, we consider the equations of continuity, motion, and energy in a Cartesian 
orthogonal coordinate system. Let p denote the specific density and v’ denote the velocity 
vector, then the equation of continuity is 
v d;p + pow’ = 0. (2.1) 
For non-isentropic flows of a polytropic gas, the equations of motion may be written 


in the form | 


pv’ OY; a 0,(pc’) == 0, - (2.2) 


*Received April 17, 1956; revised manuscript received July 5, 1956. The author is indebted to the 
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where c is the local sound speed which may be defined in terms of the pressure, p, and 


: 2p) 
C y , 
\0p/ 5 


the specific entropy S by 


and y is the constant of the gas law 
p = p’f(S). 


In addition, we shall use as energy equation the Bernoulli relation 


vid ( meee hi | = 0, (2.3) 
"\y - 1 2 


where ¢ v’v; is the magnitude squared of the velocity vector. 
$y definition, the characteristic manifolds of the system (2.1), (2.2), and (2.3) are 
those surfaces along which discontinuities in the derivatives, O;p, 0;¢C, OV; , can occur. 


Using the methods of a previous paper [2], these surfaces can be easily determined. 


However, these methods are not pertinent to the present problem; hence, the procedure 
which follows will be by analogy. By eliminating p between (2.1) and (2.2) and then by 
eliminating 0,c° in the resulting equation through use of (2.3), we obtain the following 
relation 

v'v Cg ) O,v, 0, (2.4) 
where g represents the metric tensor. Thus, the same basic equation (2.4) is valid in 


9 


the non-isentropic as well as in the isentropic case [3]. 


If the symmetric tensor a’” is defined DY 
a vv —cg’', 


lows the characteristic manifolds of (2.4) are determined by 


then for isentropic 
ann 7. (2.6) 


ector orthogonal to one family of characteristic surfaces. It can be 


where 7; is a unit 
shown that (2.6) determines the characteristic manifolds in the non-isentropic case 
[2]. From the known theory of characteristic manifolds [4], it follows that 

v cn’ + Dd’, (2.7) 


where ?’ is the unit ta 
b = (q° —c)”. (2.8) 
We shall require that both b and ¢ be positive. However, we shall not fix the sense 


of t’, n’. This means that we shall be able to use either nappe of the normal or bicharacter- 


istic cones. 


By use of (2.5), (2.7), we find that 
a’’n, = bet’, (2.9) 
a’*t, ben’ + (q° — 2c)t’. (2.10) 


Secondly, we shall introduce a few basic relations from differential geometry [5]. 
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Since n’ is a unit vector orthogonal to ©* surfaces, we may write 
k k 
Om = Si, + NU , 8,.n =un = 0, (2.11) 
. 7° : » 1 
where s;, is the symmetric second fundamental tensor of the ©” surfaces orthogonal 
to n,; and u, is the curvature vector of the n; congruence of curves. Further, let p; denote 
a unit vector field which is orthogonal to both ¢; and n; so that p; , t; , and n; form an 
ordered orthogonal right hand triple at each point. The curvature vector of the p; 
congruence of curves will be denoted by w; , and the curvature vector of the t; congruence 
by m; so that 
Ww; = p* dp 
: _* (2.12) 
k 
nm; = t O,t; . 
l'rom the decomposition of the metric tensor 
ik k i_k i ik 
g =tt+n'n'+p'p, 
we obtain the relation 
ik k k ao es 
g* d;t, = Ut d;t + n'in* 0; + p'p’ jth « (2.13) 
It is easily seen that 
> 7 i_ &k k 
t't’ 0;t, = 0, n'n Oj; = —tu, pp do; = —tw. (2.14) 
Thus, we may express (2.13) in the form 
k a k A id 
g° 0; = —ti£(w +w). (2.15) 
In this paragraph, the relations (2.7), (2.11), (2.15) will be used to express the basic 
relation (2.4) in terms of the rates of change of g and ¢ with respect to displacements 
along ¢’ and n’, and in terms of the curvatures s;, , u; , w; . The following notation for 
directional derivatives will be used 
0 ta 0 ia 0 iy 
71> c ~. =n C — =D od; . 
ot divs on ve . 
Thus, 0/dt represents rate of change with respect to displacement along ?’, and similar 
interpretations are valid for 0/dn, 0/dp. By differentiation of (2.7), we find 


OU; =C ON 4. b 0; t, + ny 0;c + t, 0;b. (2.16) 
Multiplying (2.16) by a” and using (2.9), (2.10), we obtain the relation 
jk ik « dc , al 2 2, dl —— 
a’* av, = ca™ ayn, + ba™ 0;t, + on % ae a0) + (q° — 2c) _ (2.17) 
To evaluate the term a’*d;n, in equation (2.17), we use (2.5), (2.7) and find 
a” = c(nin® — g* — tt) + Pt + cd(tin’ + n't). (2.18) 


Forming the scalar product of (2.11), (2.18), we obtain 


-~ 
bo 
_ 
e 

S 


ik k 2 isk 2 
a’ On, = chu, + b's;,t't — cM* 
where \/* is the mean curvature of the characteristic surfaces 


M* = g's, . 
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To determine the term a’’d;t, of equation (2.17), we use (2.18), (2.14), (2.15) and find 
a’* ,t, = —cbs;,t't® + c’t,w'. (2.20) 


If (2.19), (2.20) are substituted into (2.17) and the left hand side of the latter equation 
is equated to zero, we obtain the desired intrinsic form of the basic relation (2.4) 


dc . ab : 2, Ob ae ea? 
oo( 2 + - ) + (gq — 2°) — + cbt, (u* + w‘) — *°M* = 0. (2.21) 
ot on Ot 
3. The intrinsic conditions for rotational (or irrotational) motion in terms of 
characteristic variables. If e'’* denotes the permutation tensor, then the vorticity 
vector w’ is defined by 


(3.1) 


wo =e du . 


Since the ordered triad p,; , t; , n; forms a right hand system, the following formulae are 


obtained by forming the scalar products of (3.1) with the vectors /; , n; , p; , respectively 


w't; = (n'p* — n‘p') aw, , 
win; = (p't' — p*t') aw, , (3.2) 
wD; (t'n® — tn’) dw 


A lengthy but direct computation by use of (2.11), (2.12), (2.16) shows that (3.2) may 


be written as 


0c F 
vl Ch = = b(n'p np) da;t 
Op 
Ob F 
on ee bp m (3.3) 
op 
0c 0b 
p —— — bs;,tt© — ctu’. 
ot on ; 
The Eqs. (3.3) determine the intrinsic form of the conditions for irrotational motion. 
veal 


Now, we consider a form of the equations of motion which contain the vorticity [6] 


0;h, — T 0;S C400. (3.4) 
Here, 7’ is the absolute temperature, S is the specific entropy, and hy is the stagnation 


enthalpy 


where y is the ratio of the specific heats of the polytropic gas. Since the stagnation 

enthalpy is constant along a stream line (see 2.3), with the aid of (2.7) we find that 
Oho Oh, —- 
—+b—=0 (3.6) 
on ot 

and S satisfies a similar relation. 

3.4) in intrinsic form, the scalar product of this relation with the vectors 


To express (3. 
t; , n; , and p; will be found. Using (2.7), (3.1), and (3.3), the following formulae are 
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obtained 


Oho mm os 0c ob i 4k k ) 
TM (2 o — bali! — tu (3.7) 
Oho ,OS (2 ob + ok k ) 
Je 4 ae Sn a i , ? meh , 8 
on I on b ot on bsjut't om 3.8) 
Oho Os a ‘. - j 
Bot ai aS _ q oe - (b’m, + c*u,)p* — cbn'p'(d;t, — ,t;). (3.9) 
Op Op Op 


Ividently, (3.8) is a consequence of (3.6), (3.7). 
The above equations can be described by saying that they form a system consisting 
of three equations (3.7), (3.9), and (2.21) in the unknowns 


oq oq oq 

ot’ on’ Op : 
where hy and S are prescribed functions which are constant along a stream line. The 
sound speed, c, is determined as a function of the magnitude of the velocity, g, by the 
Bernoulli relation (3.5). Thus, the roles of c and q may be interchanged. In the appli- 
cations, the stream lines are unknown. Hence, the problem is to determine the functions 
h, , S, and q so that two relations of the type (3.6) are satisfied (one equation in the 
derivatives of hy and the other in the derivatives of S) and also Eqs. (3.7), (3.9) and 
(2.21) are valid. 

The curvature term, 

K = n'p'(d;t, amend 0, t;), (3.10) 
in the right hand side of (3.9) will be briefly considered. Through use of (2.11), this term 
may be written as 

K = —(n't' 0m — t'p* an, = —nit* Or + s.it'p*. (3.11) 


If the unit vector field p; is orthogonal to ©" surfaces, as is the case in plane and axial- 


symmetric flow, then 
, i = py, = 0 
OPe = Tir + PM, Pr = pw, =O, 


where r,, is the symmetric second fundamental tensor of these surfaces and w, is the 
curvature vector of the p; congruence (see 2.12). In this case, (3.11) reduces to 


K = —r,n't® + 8,;,t'p*. (3.12) 


For the special cases of plane and axial-symmetric flow, K vanishes. In both cases, 
p; is orthogonal to ~' planes and hence r;, = 0 and 7p; , ¢; are in the principal directions 
of the surfaces orthogonal to n, [7]. The ©* planes are parallel in plane flow; they 


intersect in a common line in axially symmetric flow. 


4. Canonical form of the system (2.21), (3.7), and (3.9). In this section the equa- 
tions (2.21), (3.7), and (3.9) will be written so that one equation contains only the 
directional derivative, 0/dt, a second equation contains only the directional derivative 
d/dn, and the final equation contains only the directional derivative, 0/dp. 

If (3.7) is multiplied by b and the resulting equation is added to (2.21) the following 


relation is obtained. 
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> a ¢ 5 > ” ‘ oh os . 
b* 2 d -c’ bt ae cs —— "g e 1. ( — — 7 8) f 
at b ce bti,w* + b’es;,t’t + cM* +l ry 1 at)" (4.1) 


For q > c, the equation (4.1) is the first of the desired canonical relations. 
second relation of desired canonical system, we multiply (3.7) 


To determine the 
by b and subtract from (2.21). We find that 


db ! 2», GO 4 R 4 xk 4 k 
- -+- (gq 2c) — + cbi,(2u° + w) — co M* + cb’s,,t'l 


_ U( ae -? 2), 
ot ot 


md term of the left hand side of (4.2), we use (3.5), (3.7) and obtain 


Ob Aho P _ os Ob k 9 k 
— 1)jb— —2— + (y 4 i T—-e--— cbs;,t' t —ctu |. (4.3) 

ot Ot Ot on } 
Converting 0h,/dt, 0S/dt into dh,/dn, 0S/dn by use of (3.6) we find that (4.3) reduces to 





ab Oho ON ol ' ‘ , 
Qe —2 — (y + 1)| cT — + be — 4+ cb’s,,t't® + c*bt,u* |. (4.4) 
; | 


Substituting (4.4) into (4.2) and converting dh,/dt, 0S/dt of (4.2) into dh,/dn. AS/dn 
by use of (3.6), we obtain the second canonical relation for g > c. 
. _ oD ' , 1 oO ; ok > 2 * 
bly — 3)q? + 4c’ 2g? — (y + B)c"]b°sj,t't! + (y — 1)b%°M* 
47 
; ; ; k - 
cb} (y 3)q + 4c jiu” — (y — 1)b’ct,w (4.5) 
ze Aho p : 7 ON 
y¥— 3g — & — de] — + [2¢ +i r —: 
< on on 
Evidently, (3.5 the third equation of the desired canonical system. This equation 
can | tt the n 
o”0g oh, ails as ‘ P . os 
q A - = 25 1 OM TCUuw)P + cbk, t.6) 
: Op Op op 
where K is the curvature defined in (3.10 


5. Aclass of space flows. For plane flows, the two families of characteristic surfaces 
to be right cylinders with parallel generators which are perpendicular to 


1 
can be enosen TO 
the case where on family of characteristic s irface S 


the pla e ( the flow. Here, we consider 
are paraitel planes We shall show that this condition defines a family of space flows. 
3 hat the ordered triad of unit vectors Di; l; ,n; 


We recall that in Sec. 3, we assumed that 

forms a right hand system. If the characteristic planes are assumed to be perpendicular 
to the z-axis of an 2, Cartesian orthogonal coordinate system, then n,; can be chosen 
to be sensed either in the positive z-direction or in the negative z-direction. In the first 
case, the ordered pair, p; , t; forms a right hand system when considered as part of the 
ordered triple p; , ¢; , and the positive z-direction. The second case leads to nothing new. 


Hence only the first case will be considered. 
In any plane orthogonal to the z-axis, we introduce two families of parameter curves: 


variable along t; , 8 = variable along p; . Arc length will be positive when measured 
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° . ° . 7 . 1 . 
in the directions of t; , p; , respectively. Every family of ~* curves, a = variable (z = 
constant), generates a surface. Any surface 8 = constant will be chosen 


constant, B 
. . . 1 rT 
variable, lying in © planes, z = constant. We 


so that it consists of ©' curves, a = 
choose the surfaces, a = constant, in a similar manner. The third family of surfaces 
in our new coordinate system will be the planes, z = 6 = constant. The transformation 


equations relating a, 8, 6 to 2, y, z are of the form 
xz = x(a, B, 5), y = y(a, B, 4), g = 6. (5.1) 
In any plane, 6 = constant, the arc length element is 
ds’ = (Ada)’ + (Bdg)’, (5.2) 


where A, B are functions of a, 8, 5. If 0(a, 8, 6) denotes the angle between the z-axis 
and ¢; then the unit vectors ¢; , p; have Cartesian orthogonal components 


t;: (cos 6, sin 6, 0), (5.3) 
p;: (sin 6, —cos @, 0). 


Thus, we find that 


IX A , 
= = A cos 6, Zw A sin @, 
Ja 0a oa 
(5.4) 
~ B sin 6 oy B cos @ 
— = Bsin —=- s 0. 
0p : og 
In the case of curves lying in the planes, z = constant, the sense of the principal 


normal vector is chosen so that the ordered triple consisting of the tangent vector, the 
principal normal vector, and the positive z-direction form a right hand system. The 
curvature may be positive or negative. By use of the Frenet formula for plane curves 
and the fact that the principal normals of the p, , ¢, congruences lie along & , —p: , 


respectively, we find 


‘OP, = UW, = k’t, , — 
yoo ™ = (5.5) 
t’ d;t, = m = —kp, , 
where k’ and k are the curvatures of the curves, 8 = variable and a= variable, re- 
spectively. From the equations (5.2), (5.5), it follows that 
t’ dp; p dt; eo 
R’ —= : ‘ —k nr —+. (e 0 
B op’ A da 5.6) 
Substituting (5.3) into (5.6), we find that 
tae 18 
~ Bap’ P A da 
Again, we note that the integrability conditions of (5.4) are 
00 10A 00 1 dB — 
— =- = ———- (5.7) 


da 2B ap , aB ~ A da 


The relations (5.5), (5.7), lead to the useful results 
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- l ob 
w,t* k a 
AB da ety a 
(5.3) 
unl ] 1 OA 
sp = —k oe oe 
AB op 
By use of (5.4), we see that the matrix of the transformation (5.1) is 
ON 4) y ra) 2 P | 
- A cos 0 Asin é OQ 
0a Oa Oa 
| oa OY Oz | , 
oa ae Bsné@ —Beosé 0 (5.9) 
OD 0p op | 
ahi OY Oz | Ox OY | 
06 06 06 06 06 
Since the matrix of da/dzx ete. consists of the reduced cofactors of the determinant of 
the matrix (5.9), we find 
Oa OB 06 cos 6 sin 0 0 
Ox Ox Ox A b 
Oa OB 06 | sin @ cos @ 0 5 | 
ae | . (9.10 
Oy oy oY ] B 
Oa Oop 06 
Qy3 a l 
Oz Oz 02 
where 
sin G Ov cos 6 0x 
a a ee . gee 
{ do A 00 — 
(5.11) 
cos 6 dy sin 0 0x 
Qo3 : thee os 
: B 00 B do 
Thus, we find by u the chain rule and (5.10) 
0 Oo 4] 4] ~ 
2) , Qi3 7 = Ges ~ — << (5.12) 
O8/ «9 0a op 06 
The a | the coefficients a ~a can be seen from the following arguments 
By use of (5.9). (5.11 ve find that 
1 (dx 0x Oy Oy 
a ( ap : . (Q){ 1) . 
A 0a do Ja d0 = 2 
{ 13 
Q23 = Z (= = 4 SY + (0)(1) }. 
B OB 06 OB 06 - 


denotes the cosine of the 


constant, 6 


If eos (a. 6 
constant, and a 


between the curves, ¢ 


consti 


‘ constant, 6 = 


Aa, cos (a, 


angle between the curves, 8 constant, 6 


unt and cos (8, 5) denotes the cosine of the angle 


constant and a= constant, 8 = variable then 


Baz; = cos (§, 4). (5.14) 
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Again, we note that if we solve (5.11) for dy/06, 02/06, we obtain 


= = Ba, sin 6 + Ads cos 6, 

c (5.15) 
) ; 

2 Aa,; sin @ — Baz, cos 06. 

06 


From (5.4), (5.15), we obtain a set of integrability conditions which must be added 


to the equations (5.7). These are 


2 (Bay;) — Aay 2 = -AS, (5.16) 
2 (Aa,:) + Bass a ~ at :' (5.17) 
2. (Bi) ~ Ate - Ps 2, (5.18) 
= (Aays) + Bass Ss =- “ (5.19) 
From (5.14), we see that if the angles between the coordinate lines, a = variable and 


6 = variable, 8 = variable and 6 = variable, are 7/2, then Aa,, = Ba,; = 0. The above 
relations show that in this case, 0, A, B are all independent of 6 (as is to be expected). 
The relations (5.7), (5.16) through (5.19) imply that the Riemann tensor vanishes in 
the a, B, 6 system. 

To determine formulas (4.1), (4.5), (4.6), (3.6) in terms of the congruences t’, p’, n’, 
we must evaluate the curvature terms. Since, the n’ congruence consists of straight lines 
parallel to the z-axis, the curvature vector of these curves, u, , vanishes. Further, the 
second fundamental tensor of the planes (z = constant) is s;, = 0; also, the mean curva- 
ture, J/* of these planes vanishes. Again, we consider the curvature, K, of (3.10) 


K = n'p'(d;tk — Ot;). 
Since the vector ¢; lies in the plane, z constant 
n'p* d,t; = 0. 


Further, n’d,;t, represents the directional derivative of ¢, in the z-direction. From geo- 


metric considerations, or use of (5.3), it follows that 


By use of the above results and (5.12), the formula for K becomes 


_(,. 2%, 0 a0 
K= (a1 a + 93 9g a0), (5.20) 


Finally, the directional derivatives 0/01, 0/dn, 0/dp become 
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0 0 0 1 a) 
* B 
ot Ja Op 0p 
I (5.21) 

re | al F 0 0 0 

@i3 7 - Gas = el 

on O02 — \ Ja Op 00 


Now, we express the system (4.1), (4.5), (4.6), (3.6) in terms of the metric quantities 
ed Gia 5 Hh 6. and the partial derivatives 0/da, 0/0B, 0/06. By use of (5.8), (5.21), 


we find that (4.1) becomes 








0 ( , 0 Oh,, one Os - 
b—4 =c—mnB+- =, (5.22) 
da b Oa Ja Ja 
Further, by use of (5.21 5.8). relation (4.5) reduces to 
: - ob ob Ob 
- bl (y 3)q t lc || a lox = 
Ja Oop 00 
1 O > 2 : - 
\¥ L)b'cA — In B [(y d)q (¥ we 
Ja “= 
(5.23) 
Oh, Oho Oh, oe : = 
Qi3 f a - fi 24 1 a are | 
da Op 00 
Os Os Os 
a ; t+ a ~ 
Ja OD 00 
By use ol Qj 20 “ ‘}. 5.8 » we find that | 'f) becomes 
Og al on O . O06 00 00 = 
‘ewe | > b In A + cbhB) a, - Gas = 9.24) 
0p rate! OD op Ja ri]e} 00 
Finally, (3.6) reduces to 
oh on oh - 
ca + bA)- ca rec Q. (5.25) 
Ja Op 00 
From (5.25), it follows that the differential equations ol the stream lines are 
Ada dp ? es 
ao, 5.26) 


1a,. + (M° 1)° a 


where M is the Mach number, q/c. Since Ada, Bdg, dé are the are length elements along 
5), cos (8, 4), it follows from (5.26) that 


. Bay, are COs | 


the coordinate lines and Aa, a, 
in general the stream lines are space curve 

6. A family of isentropic Beltrami flows with plane characteristic surfaces (helical 
stream lines, concentric circular bicharacteristics). We shall consider those flows of 


See. 5 for which 


h constant, S constant, a, A3 0. (6.1) 


From (3.4), it follows that in any isentropic flow for which fy is constant, the vorticity 
vector is parallel to the veloctiy vector. Such flows are called Beltrami flows and are 


more general than irrotational flows. As noted earlier, (5.16) through (5.19) imply that 
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OA OB 00 i 
~ —=-7,=0, (6.2) 
00 00 06 
that is, the curves, 6 variable, are orthogonal to the planes, z constant. Stated in 
other words, every family of © bicharacteristics (6 = variable, a = variable) forms surfaces 
j ; 4 "9 ° 1 : 
with straight line generators parallel to the z-axis; also, every family of ©” orthogonal tra- 
gectories of the bicharacteristics (6 variable, 8 variable) forms similar surfaces. From 
5.23), (6.1), (6.2), we see by use of a simple argument that 
or OB »¢ 
—t = —— = 0. (6.3) 
00 Ja 


rom (6.2), (6.3), we see that, B |, by proper choice of scale factor along the 6 = 


variable curves. Further, (5.7) implies that the orthogonal trajectories of the bicharacter- 


tics are straight lines or k’ 0, 0 O(a) and 
A Bo’ + f(a), (6.4) 


le and f(a) is an arbitrary function of a. Returning to (5.22), (5.23), 


where 6 0 
5.24), we see that (6.1), (6.2) lead to 
q q(B), (6.5) 
M”°* d M° 0 = 
at ya ™ 3 Fy pat = 35. (8.6) 
where M the Mach number. Since A/ is a function of only 8, we see that in (6.4) 
O(a) = ca, f(a) = cx, , (6.7) 
1ere C, , ¢ ibitrary constants. Integrating (6.6) and choosing a proper scale factor 
along the « viable curves, we obtain 
2+ (7 1)M’ (yen . 
cae aa es (6.8) 
lo inter} 6.8). we must determine the significance of A. From (6.7), (6.4) we find 
A =c,(8 + ¢,). 
Since, B , we obtain from (95.4 
x B +c.) sin 0, y —(B8 +c.) cos 6. (6.9) 
From (6.9), we d that, c, + B r, where r is the radial distance variable in polar 
coordinates. | e introduce the polar angle @ @ —7/2 into (6.9) then 
x =r cos ¢q, y =rsin ¢. (6.10) 
Thus, (6.8), may be written as 
2 os ’ )Mn = Crt, (6.11) 
where C is a constant. It is easily shown from (6.11) that for y > —-1, M 
decreases a increases, M | as r approaches infinity, 17 approaches infinity when 


r= ily —1)/C] 
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From (6.9), we see that the bicharacteristics, r = c, + 8 = constant, are concentric 
circles. By integrating (5.26), we obtain 


Ge-p* 
, r 


od = z+e3, r = constant, (6.12) 


where ¢c; is an arbitrary constant. Thus, the stream lines are helices and cut the generators 
of the cylinder, r = constant, in the constant Mach angle y, where 


tan y = (M’ — 1)”. 


By use of (3.3), and the relations (see 5.5) 


‘ dc Ob _ Trae 
K=%, = 8, = eal lia 0, pm=——; 
we find that 
dc db ) j 
te Mea lOa by (6.1 
w a te (i a < mn 6.13) 


Thus, the present class of flows are rotational flows. Evidently gq is constant along each 


stream line [8]. 
By assigning various appropriate values to a3 , @23 , other flows may be obtained. 
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WAVES PRODUCED BY A PULSATING SOURCE TRAVELLING 
BENEATH A FREE SURFACE* 


BY 
H. S. TAN 
University of Notre Dame** 


Summary. The propagation of surface waves induced by a pulsating source travel- 
ling beneath a free surface is investigated by introducing Rayleigh’s dissipative force. 
By defining a disturbance velocity potential Rg(x, y) exp (t#t), it is shown that g(x, y) 
has to satisfy throughout the flow field not only the Laplace equation, but also a differ- 
ential equation resulting from the analytic extension of the free surface boundary con- 
dition. The solution is obtained in the form 


g(x, y) = f(z) + gf), z=axa+1y, 


/) is a complex harmonic function of x, y but is not an analytic function of z 
The wave propagation is found to depend on a parameter 7, which is the ratio of the 
pulsation frequency w of the source strength to the fundamental frequency wo of the 
surface wave produced by a constant strength source travelling at the same speed. The 


OAT, 


case T 0, corresponding to the travelling source of constant strength, does give a 
single undamped harmonic wave train on the downstream side, of wave length 2zc*/g, 
or frequency w)), as is expected. For 0 < + < 1/4, there are four, and for r > 1/4, two 
undamped harmonic wave trains of different wave lengths on the downstream side. 
It is further observed that there exists a critical frequency at r = 1/4, at which frequency 
resonance phenomena occur. Thus violent disturbance is predicted at r = 1/4 by the 
present analysis. No disturbance is found to propagate infinitely upstream. This result 
evidently justifies the imposition of an asymptotic upstream condition of ‘vanishing 
disturbance at infinity’ to replace the effect of the dissipative force in rendering the 
solution unique 

Velocity potential. Consider a two dimensional incompressible fluid of infinite 
pe and extent, beneath whose free surface a point source is travelling at a constant 
speed c. In an (x, y) coordinate system that moves with the source, and under the hypo- 
thesis that the resulting fluid motion is irrotational, one can define an absolute disturb- 
ance velocity potential ®(r, y, 4) through the differential equation 


P.. + ®,, = 0 (1) 
and the boundary condition over the free surface Y(y = 0) 
?, — c®, = —q@gY, 


Y, —<¢r, = ¢,. 


4, — Wh,, + b,, + gh, = 0. (2) 


When the strength of the travelling source fluctuates harmonically, we may write 
the velocity potential in the form: 
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P(x, y, t) = Re(a, y) exp (wt) (3) 
and obtain the following differential system for ¢: 
Gzz + vy = O, (4) 
vo, — 2irve, —veo t+ re = 0, y = (, (5) 
where 
tT = we/g, yv=w/g. 


It is a well known fact that the solution of this problem is indeterminate, as free 
harmonic surface wave trains of appropriate wave lengths can be added without dis- 
turbing either the differential equation or the boundary condition. This indeterminacy is 
removed by introducing an artificial dissipative forge —u V ® and then letting » approach 
zero in the result. This artifice was first introduced by Rayleigh’, and has since been 
successfully employed by various investigators*’*’® for studying surface waves set up 
by moving sources and sinks of constant strengths. 

3y adopting this scheme, the differential equations (1) and (4) are not affected, but 


the boundary conditions (2) and (5) are replaced by 


®,, — 2ch,. + c’o., + p®, — wet, + g?, = 0, (6) 
wy, — Lirr(1 — iB)e, — (1 — 218)o + 7’e.2 = 0, (7) 
respectively for y = 0, where B = u/2w. 

Let us now write the spatial part of our velocity potential (with moving source 
located at ¢ = & + in) in the form 


g=In|z-—¢|—In|z— &| + Gi, y;é, 0), z=xz+y, (8) 


where G(z, y; £, 7) is harmonic on the lower half plane. Noting that on y = 0 


0 ; ; 0 acai ss 
rs Injz—¢] = ——In]z2- Fl = 9e- ",;" 
oY oY (9) 
Q” 3” - 
—Inijz-—¢/=—lmnlz-¢| 
02 ,' : | Ox" sis 


and substituting (8) into (7), we obtain the following differential equation for G on 
y = 0: 

vG, — 2irr(1 — i8)G, — v'°(1 2i8)G + 7°G., = 2v =. In|z— §|. (10) 

oy 
As both sides of (10) involve functions that are harmonic in the lower half plane, (10) 
may be required to hold for the entire lower half plane by analytic extension. It can 
indeed easily be shown that a general solution of the Laplace equation is given by 
f@+9@ 
Now we have the well known relation 


fo) 


@- pt =if expl—ice—- Olde, ytn<0 (11) 


from which the following expression can readily be deduced: 
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%»—in|z—-§| = -» [ exp [«(y + m)]fexp [ix(x — 8)] 





oy (12) 
+ exp [—ix(x — £)]} de. 
Combining equations (12) and (10), we arrive at 
7. exp [ky +n] __ a 
aii [ Te + fal — B) — Ipe +P — ae) OP [tate — Ol de (13) 


ica / ss exp [x(y +] _ 
© do ow — [2711 — 78) + lox + (1 — 278) 


exp [zx(x — &)] dk 





which indeed satisfies both the differential equation (4) and the boundary condition 
10), and is complex harmonic, in the form f(z) + g(2) as expected. 
Finally; we put the disturbance velocity potential in the following form: 
l 


ig —* + G(x, y; g, } exp (wut). (14) 
$ 


p~ = aI 


x jx 


2. Wave propagation. From the disturbance velocity potential (14), it is easy to 
see that the only term that contributes to the propagating wave is G(z, y; £, ») as given 
by (13). Now both integrals of G are singular; principal values must therefore be con- 


sidered. 
The integrand of first integral has singularities at the roots of equation 
rx — [1 — 2r(1 — 78)Jvx + e°(1 — 278) = 0 
( at 
= {1 — 2r(1 — 78) + (1 — 47 + 4778 — 47°”)? Wv/2r’. (15) 
When 'e) 0, the e roots are 
Kyo = [1 — 27 + (1 — 47)" )p/27’. (16) 
Similarly, the integrand of second integral has singularities at the roots of equation 
x — [1 + 27(1 — 78) vx + (1 — 278) = 0 
r at 
[1 + 2r(1 — i8) + (1 + 47 — 4778 — 47°8°)'”? Jv /27’. (17) 
When 6 U, these roots are 
K3 4 (1 + 27 + (1 + 47)" )p/27’. (18) 
TABLE | 
R IJ 
Ki Ke K3 Ka || Ki K2 Ks K4 
0 | g/t 0 {g/t 0 lot 0 |O- O 


Ml (+) (+) | (+) (4+) 0+ 0+ /0- 0- 
Yes) (+) (+) | (+) (4+) I (4) (-)] 0-— O- 
ywe ey (=) C—7 | CH) CO [or (3 | . Ge 
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A careful investigation of (15) and (17) shows that as long as 6 ~ O, the following 
results hold as 8 approaches zero. 
In the table, S\ mbol 0+ indicates the following limiting behavior: 


Lim 9K 0 (19) 


As will be shown later, this limit physically corresponds to the case of an undamped 


propagation of a harmonic wave train of wave length 27/x; . Indeed, by writing (13) 
in the form 
: yp [ r exp [x(y - n) 
Gr: ; | I | exp LK\a t | dk 
T Ld (kK — K(k — ke) 
F exp lk(y t n) | ] 
t | a ~ exp [¢x(e — £)] dk 
P K Ka)(K — Kg : 
y . | | \ ‘ , 1 9 
, | - ) exp [x(y t n) | exp | ’“/sr— £)|] dk (20) 
tT ‘ h Ky K Ko/ ; ‘ 
y ‘f l | , 
| ( ] exp [x(y + 7)]| exp [ix(a — &)] ds 
| + Lr 4 Ty ae I K K,/ ; 
yi(] Lr I, + (+ 47) a 
and putting 2 + im, all the integrals in (20) can readily be evaluated by proper 


1O 
choice of the contour of integration and application ot Cauchy’s integral theorem in 


the complex 2 plane. To fix the proper contour of integration, observe that, in order to 


make the integral meaningful, in addition to the necessary part of the positive real 
axis (m 0,0 </ ©) the remaining part of the contour must consist partly of the 
imaginary axis (/ 0), (as the integral diverges along negative real axis due to the 


presence oi the factor exp [x(y + 7m)| under the integral sign), and partly of a circular 


I 


are of infinite radius over which the integral vanishes. The choice of the propel quadrant 


for the contour I’; for 7; depends entirely on the sign of the difference x — £& When x — & 
1 tm 
—, —, 


wd 

















K A 
O ~ ro) 
Ys 
2 
i; 
(a) zr —E>0U (b) x —E<O 
I for Jy(«1, x2) Tl. for IJo(xs, «s) 


Fic. 1. Contours of integration 
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‘< negative, this part of the contour must be taken in the first quadrant of the complex 


Q plane for J; , and in the fourth quadrant for J, , in order to insure convergence of the 


integral; and vice versa when x — £ is positive (see Fig. 1). Whenever a pole k; (j = 1, 2 
for 1, ,J 8 4 for J.) lies inside its corresponding contour of integration T; (¢ = 1, 2), 


there will be a residue term by Cauchy’s integral theorem. As will be shown later, this 


‘ves rise to a harmonic wave train which propagates in an undamped or damped manner, 


9 
according to whether the pole is on or off the real axis, when 6 approaches zero. 
Taking cognisance of the above remarks, and applying Cauchy’s integral theorem, 


we readily write the general expressions for the different cases as follows: 


a a—-t>O0 
- exp [x(y 4 ) ; - 
| pity | exp [—ix(x — &)] dk (j = 1, 2) 
J0 es 
—2Qri exp [x)(y + 0) — t«K\(e — £)] (21) 
r” exp [—im(y 
+- 2 bn y + mI exp [—m(x — £)] dm, 
Jo m — 2K; 
“= axp [xly 4 ) : . 
] ! 1D) exp fix(a — &)) dk (j = 3, 4) 
4 r - K 
Ori exp [K;(y +m) + x(x — £)] (22) 
*° exp [am ) 
P| y + 9)| exp [—m(x — £)| dm, 
Jo Mm + UK; 
(b) a 0 
*” exp Ky + ) : 5 
I, PL ™) exp [—ix(x — &)] dk (7 = 1, 2) 
h — LN 
271 exp [k (y 4- ") - ix(x — £)] (23) 
r° exp [im(y +n) ‘ 
-+- = = y + 9] exp [m(x — £)] dm, 
Jo m + 24; 
-? exp [«(y + 7) ; - a , 
| pinky | exp fix(a — &)] dk (7 = 3, 4) 
— th 
(24) 


Qi exp [x,(y + 2) + i«i(x — §)] 


-}- [ exp | tet ”)| exp [m(x — £)| dm. 


m — 2K; 


70 


It can be shown that all the definite integrals obtained above are regular, and they 
decrease rapidly with increasing | X — §}; hence they can only give rise to a local disturb- 
ance. Thus it is the residue term only that can lead to infinite propagation of harmonic 


wave trains. Indeed, noting that all the residues are given in the form 
Qi exp [x,(y + n) + ix (a — &)], 


we can readily conclude that an undamped harmonic wave train appears whenever we 
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have lim,.o 
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Sk; 0+. By further taking into account the detailed distribution of 


singularities in the various frequency ranges as listed in Table I, we obtain the following 


results: 


(fn) lon c 


ili 
Thus, on 
local one 
(b) For x 
Ie 


wave train 1s 


Ll 


Thus, on the 


trains propag 


£ > 0 (i.e. upstream of the travelling source): 
Since both x, , x, lie outside of contour of integration [T, for the entire 


range of parameter r, I, has no residue. 


for r < 1/4, both x, , x. lie outside of contour T, , 7, has no residue. 
for 1/4 < r < 1/2, x. lies inside T, , so that 7, has a residue term at x, 
However, as lim,.. 9x. + O, this residue term leads to a harmonic wave 
train that is damped out exponentially with | 2 — & |. For by writing 
K m in, where n > O, we have 
| @) - Kol £)] 

2r7 exp [m(y + 7) n(x t)] exp {—7a[m(x — §) + n(y + )]} 
which verifies our statement 
for 1/2 r, both x; , ka lie outside of contour [, , so that 7, has no residue 


the upstre tm side of the ti ivelling source, the disturbance is always a 


0 (i.e. downstream of travelling source 


For the entire range of 7, both «s and x, lie inside the contour of integration 


lr, , so that J, has two residue terms at «x, and x,. As limg.5 Ik3 , kK; 0, 
there will be two undamped harmonic wave trains of wave lengths 2x 
and 2x/x, on downstream side which propagate to infinity. Indeed, it 


CAN CASILY De hown that 


ii i h 
which varifies our statement 
lor 7 1/4, both x, and 1 lie inside of eontour ir; , 20 that ‘; has two 
residue terms at x, and x,. As limg.o Ik, , K 0, two undamped harmonic 
we trains of wave lengths 27/x, , and 27/x, propagate to infinity on the 


asily be seen as follows: 


downstream side. This ec: 


y 4 1. ik, r - £)] exp [Ko(y - 7 LK (x £)]} 
§ | | a 
exp [iem(y + n) + ma £)| dm 

E \m -7 K m -- wk 
for 1/4 r < 1/2, only k, lies inside the contour T, , so 7, has a single 
residue there. However, as lims.o 9x; # 0, the corresponding harmonic 
rapidly damped out. 
for 1/2 r, both Ky and x. lie outside of the contour T, , so that 7, has no 
residue 


downstream side of the travelling source, there appear multiple wave 


ating to infinity. The number of these infinitely propagating wave trains 


is determined by the parameter r. 
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The case of a travelling source of steady strength can easily be checked with Rayleigh’s 
result. Here the poles are at x, and x; (Table I). Now, for points upstream, x — & > 0, 
x, lies outside the contour T, , x; outside T, , so neither 7, nor 7, has any residue. For 
points downstream, « — £ < 0, however, both «, and x; lie inside their respective con- 
tours I’, and T, , and so both J, and J, have residues. As limg..o $x; , «; = 0, these residue 
terms represent undamped harmonic wave trains on the downstream side, of wave 
lengths 2x/x, and 2x/xs, respectively. However, from Table I, we have x, = x, = g/c’. 
Thus we only obtain a single harmonic wave train of wave length 27c’/g, as expected. 

Summarizing, we have the following statements. 

When a point source of constant strength travels with speed c beneath the free 
surface of a two dimensional, ideal incompressible fluid of infinite depth and extent, 
aside from a local disturbance that travels together with the source, there is also pro- 
duced a fundamental harmonic wave train of wave length A, = 2mc°/g, or of frequency 
rn g/c, which propagates to infinity on the downstream side. There is no disturbance 
propagating to infinity on the upstream side. This was first observed by Rayleigh, and 
agrees exactly with the present analysis (r = 0). 

When the strength of this travelling source fluctuates harmonically, then, aside from 
the local disturbance that travels with it, there are produced on the downstream side 
multiple wave trains propagating to infinity. This phenomenon may be regarded as a 
splitting of the original single wave for constant strength source by the fluctuation of 
source strength. The number of these infinitely propagating wave trains is determined 
by the value of a parameter 7, which is the ratio between the fluctuation frequency 
w of source strength and the frequency of fundamental wave w) explained above 
(r = w/w). When +r < 1/4, there are four wave trains of different wave lengths propa- 
gating to infinity on the downstream side; when r > 1/4, there are only two. No 
disturbance propagates to infinity on the upstream side. 

Inspection of formula (20) reveals an interesting phase of the wave production by a 
travelling pulsating source. It is seen that the amplitude of both local disturbance and 
propagated wave trains contributed by the /,-term exhibits a resonance behavior as 
the parameter 7 approaches its critical value r = 1/4. Thus violent disturbances should 
be expected at r 1/4 according to the present theory. 

The present analysis also agrees with the physical intuition that no disturbance 
could propagate to infinity on the upstream side. This result can evidently be con- 
sidered as a justification of our introducing the asymptotic upstream condition of 
“vanishing disturbance at infinity” as an alternative to Rayleigh’s dissipative forces 
in rendering the solution unique. 
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BOOK REVIEWS 


Randwertprobleme der Mikrowellenphysik. By Fritz 1. Borgnis and Charles H. Papas. 
Springer-Verlag, Berlin, Géttingen, Heidelberg, 1955. xv + 266 pp. $11.50. 


This book is an excellent short treatise on the modern applications of electromagnetic theory to 
antennas, wave guides, and scattering. The well-organized collection of problems appears for the first 
time in book form, although much of it could be found by an exhaustive search through the journal 
literature 

The connecting link between these various problems seems to be the nature of the mathematical 
techniques needed in their solutions. There is extensive use of Green’s functions and their various repre- 
sentations, integral equations, and stationary principles of the ‘‘Levine-Schwinger” type. Computational 
details are often omitted but numerical results are given in graphical form and there is usually a careful 
discussion of the phy Sic il iwnihcance of these re sults 

The very recent and important results on high frequency scattering (work of Keller, Rubinow and 
Wu, Kodis) apparently could not be included in time for publication. These papers must have reached 
the authors in proof since there is a footnote throwing doubt on the formula (8.7 on page 51) for the 
scattering cross section of a circular cylinder. Indeed this result must be discarded and replace d by the 


one found by Rubinow and Wu. It would also seem that the integral equation (5.11) might be discussed 


somewhat more carefully ince the kernel, which is a second derivative of the Green’s function at the 
singularity, does not « t. The proper limiting procedure in this case has been explained by Bouwkamp 
in his review article in the “Proceedings of the Roval Society. 


Mention should be made of the authors’ lucid stvle: their German is surprising] ens to follow 
The printing is of the highest quality and altogether this volume represents a suecessful contribution 
to the literature in the wld 


I. STAKGOLD 


Numerical analy Idited by John H. Curtiss. Proceedings of Symposia in Applied 
Mathematics, Vol. 6. MeGraw-Hill Book Company, Ine., New York, 1956. vi + 3038 
pp SO.75 


The volume contai ll but two of the papers presented at the Sixth Symposium in Applied Math- 
ematics of the Americ \Viathematical Society (Santa Monica City College, August 1953). As the space 


iles out a discussion of the individual contributions, the diversified contents 


available for this revic 
of the volume is indicated by the following list of papers: Computational Problems in the Theory of 
Dynamic Programming (R. Bellman), Some Methods for Solutions of Boundary-value Problems of 
Linear Partial Differential Equations (S. Bergman), Computational Aspects of Certain Combinatorial 
Problems (R. H. Brus Numerical Results on the Shock Configuration in Mach Reflection (D. R. 
Clutterham and A. H. Taub), Some Applications of Gradient Methods (J. W. Fischbach), Some Quali- 
tative Comments on Stability Considerations in Partial Difference Mquations (S. P. Frankel), Approxi- 


mations in Numerical Analysis—A Report on a Study (C. Hastings, Jr., J. Hayward, and J. P. Wong, 
Jr.), The Conjugate-gradient Method for Solving Linear Systems (M. R. Hestenes), Number Theory 
on the SWAC (Ie. Lehmer), The Assignment Problem (T. 8. Motzkin), The Method of Steepest Descent 


P. C. Rosenbloom), Function Spaces and Approximation (A. Sard), Some Computational Problems 
in Algebraic Number Theory (O. Taussky), Machine Attacks on Problems Whose Variables Are Permu- 
tations (C. B. Tompkins), Best-approximation Polynomials of Given Degree (J. L. Walsh), Recent 
Results in Numerical Methods of Conformal Mapping (S. E. Warschawski), On the Asymptotic Trans- 
formation of Certain Distributions into the Normal Distribution (W. Wasow), Error Bounds for [igen- 
values of Symmetric Integral Equations (H. Wielandt), On the Solution of Linear Systems by Iteration 
(D. Young 

W. PRAGER 


(Continued on p. 282 
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ON NON-EXISTENCE OF RECTILINEAR MOTION IN PLASTIC SOLIDS 
AND NON-NEWTONIAN FLUIDS* 


BY 
D. E. STONE 


Naval Research Laboratory 


1. Introduction. Tor the class of materials described by Rivlin’s theory of non- 
Newtonian fluids [1], Mricksen [2] has shown in the case of steady rectilinear motion 
that the speed is overdetermined unless the two scalars which occur in Rivlin’s theory 
satisfy certain definite relations. When the speed is overdetermined, Ericksen con- 
jectures that the only types of steady rectilinear motion possible are either rigid or 
such that the surfaces of constant speed are planes or portions of right circular cylinders. 
We show here that this conjecture is correct for the class of Rivlin fluids characterized 
by the requirement that the stress be unaffected by uniform changes in the speed. This 
requirement is one which is frequently imposed on theories of plasticity and in this 
case produces a theory very similar to the theory of perfectly plastic solids. 

In addition, we consider the case of steady rectilinear motion through closed cylinders, 
wherein the speed is not overdetermined. Restricting our attention to the case in which 
the stress and speed are well behaved, we show that certain conditions on the behavior 
of the sealars which oecur in Rivlin’s equations must be satisfied in order that such a 
motion exist in every closed, smooth cylinder. 

In the case of steady non-rectilinear motion through a tube of elliptical cross-section, 
Green and Rivlin [3] have obtained approximate solutions for a particular class of 
materials which are included in Rivlin’s theory. Here, there is a rectilinear motion 
parallel to the axis of the tube on which is superimposed a vortex-like motion in each 
of the four quadrants of the tube’s cross-section. 

2. Equations of motion. In Cartesian tensor notation, Rivlin’s equations are 


t, —pé,, + FAI, WD d,; + FLT, UD dy d,; (1) 
liz.; + ph; = p(dv,;/dt + 9;,,0;) (2) 
v,; =0 (3) 


where /,, is the stress tensor, d,; is the rate of deformation tensor, defined in terms of 
1 III = det. d,,; , 


the components of velocity v; by d;; S (v;,; + ;.;), 211 = — d,;d,; , 

h, is the body force per unit mass, p is an arbitrary hydrostatic pressure, 7, and F, 
are scalar functions of the indicated arguments, and p is the density of the material 
which is assumed to be constant. 

The analytic equivalent of the requirement that the stress be unaffected by uniform 
changes in the speed is that the function f,; (d,,) = F, (IL, WDd;; + F. (UI, WDdad,, 
satisfy the relation [,; (ad,:) = f,; (dk:) for alla > 0. In this event, it can be shown 
that the scalars F, (II, IIT) and F, (II, ILI) are of the form G, (III/—II*”) [—2I1I-'”7] 
and G, (IIL/—II*”*) [—2I1]~' respectively. Requiring that Eq. (3) hold, we have for 
the case of rectilinear motion where each particle moves parallel to the x* axis of a 
rectangular coordinate system at constant speed, v1, = v. = 0, vs = 2f(x', x”); so that 


*Received June 27, 1956. 
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the quantities II and III are given by II = —[f,; + f,2] and III = 0. Restricting 
ourselves to the case in which —JJ > 0, we then have from Eqs. (1) and (2) respectively 
! 0 O fi; | 
-_ G0) | “| 
; Pp 63; || 4 a 811)” || 0 0 J 42 
175 Ff. OT 
(4) 
7 P i] 
t fie Safe 0 
G.(0) f F f? 0 1] 
; iJ .2 2 
[—2II] |" hs 
| O 0 fits 
and 
t:;, + ph; = 0. (5) 
Assuming that the body forces h; are conservative so that h; can be replaced by ® ,; 
and carrying out the differentiation called for in Eq. (5), we obtain the following three 
equations which describe the motion 
(p p®) , G.(0){(f7:[—211}- 9.1. + OS af.2[—211").2} (6a) 
p — p®).. = G0) {(f sf 2[—211])").. + (f2[—211]").2} (6b) 
p — pd G,(O){(f 211}-'””).. + (f.2[—2I11]}"'””) 2}. (6c) 


Since f is a function of 2’ and 2’ only and G,(0) and G,(0) are constants, it follows from 


Eqs. (6) that p p® is of the form 


If we assum<¢ that the rate of dissipation ol energy, Lda; , 18 non-negative in the region 


LLilit 


of motion, then it follows from Eq. (4) that G,(0) is positive. When this is the case, the 


curvature of the curves of constant speed is constant, since the bracketed expression 


in the right-hand side of Eq. (6c) is the expression for the curvature of the curves 


f(z’, ] constant [4]. We note here that this result depends only on the form of 
F, (11); namely, A{—U “ where K is a positive constant. 

3. The case G,(0 0. When A = 0, we obtain a solution for any function f(x’ , x 
for which the curv: "4 constant are circles of constant radius. For example, 
the function f(z", 2 e+ (2 — x’)'? — 1] gives a solution in the region bounded 
by the planes 7 and x’ . Outside of this region the curves f(z", x) clearly 
intersect so that tl ution is meaningless. It is easily seen in such a motion that the 


material cannot adhere to the wall of a closed rigid cylinder. 

When A 0, the curves f(z’, x”) = constant are straight lines (not necessarily 
parallel), and any function which has this property will give a solution. The most common 
example of motion of this type is the case of simple shear wherein the curves of constant 
speed are parallel straight lines. 

4. The case where G,(0) ~ 0. We show here that the curves of constant speed are 
parallel or concurrent straight lines, where we have assumed, as in the case G,(0) = 0, 
that —II > 0. 

Let us make the coordinate transformation f‘ = f‘(z', x”), where the curves f' = 
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constant are the curves of constant speed and the curves f? = constant are their 
orthogonal trajectories. The components of the metric tensor for this system are given 
by g" oni (af'/dzx')? -- (af'/da7)?, g”” sie ” = 0, g”” “= (af?/dx')? +- (af?/dx’)’, g"' — 
1/9, 9° 1/goo , and Gyo = go, = O. In terms of these quantities, the covariant com- 
ponents of a unit normal to the curves f'(x', x”) = constant are given by 


n, =(g'')"'” and n, = 0. (9) 
The contravariant components are given in terms of the covariant components by the 
relation n' q'’n; , so that 
n' =(g')'” and n’* = 0. (10) 
In this coordinate system, Eqs. (6a) and (6b) are given by 
(p — p®),, = 4G,(0) 0/af' (log gz.) (11a) 
(p — p®).. = 4G,(0) 0/df? (log g,:). (11b) 
Integrating these equations and equating the results, we obtain 
log g1:9o2 = H,(f') + H.(f’) (12) 
so that the product g,,g22 1s given by 
tities << (13) 
Now, let us make the coordinate transformation u' = u'(f’), uw? = u?(s*) which is 
defined by the relations [df'/du'}? = e~*” and [df?/du’]? = e-“*”” so that the com- 


ponents of the metric of the new system are given by 


Hy (f*) -H4(f?) . 
és dire ' Ges = Gxt" |, «aNd Gs = Gy, = 0. (14) 


In any coordinate system, Iq. (6c) is given by 


n'; = A(2)'’?/G,(0) (15) 
and in this system, it becomes 
4(a'')'*a?* day,/du' = A(2)'”?/G,(0). (16) 
Letting B A(2)'*/G,(0) and substituting Eqs. (13) and (14) in Eq. (16), we obtain, 
since a" 1/a,, and a” 1/dz. , 
(a22)~'”? daz,/du' = 2B. (17) 


Integrating, we obtain 
Az, = [Bu' + P(w’)f. (18) 

Now, u' and uw’ are coordinates in the plane, so that we can make use of the requirement 
that the Riemann tensor vanish. Since, from Eqs. (13) and (14), a,,@2.2 = 1, this require- 
ment is satisfied, as is easily shown, if and only if the following partial differential equation 
is satisfied 

d°a,,/du” + 8’az,/du'" = 0. (19) 
From Eqs. (18) and (19), we obtain 


0°a,,/du”” = —2B’. (20) 
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Integrating twice, we obtain 


» 99 1 2 
Ay; Bru  Oi(u')u? + Q.(u'). (21) 
“* . 1 . 
Since a,,a | for all values of uw and u” , we may write 
0 a,,/0Uu (dw ) aa /du\ou). (22) 
Making use of eq 18) and (19), this relation becomes 


OBIBuw + PQ’) A&P) /du 24B/Bu' | P(w)] ‘[dP(u’) ‘du’ |’ Q. (23) 


Since a,, and ay. are never zero, we multiply by [Bu' + P(u’)]’, so that Eq. (23) reduces to 


6BIBul 4 P(u*)] d’P(u°)/du 24Bl[dP(u’) du’ | = (). (24) 


This relation is clearly satisfied for all uw’ and wu’, if and only if B = 0 or P(u°) isa constant. 
However, when P(u constant, it follows easily from qs. (18) and (21) and the 
condition a; do. | that B = O. Thus, Eq. (24) is satisfied if and only if B = 0. 

When B 0, it follows from Eq. (17) that a,, is a function of u° only; hence since 
must also be a function of vw only. We accordingly make the coordinate 


A; Ago lL. Gis 
, : ' ;' ’ = an 
the relations du /dt Ll and du-/dv (Ay») ‘72 The com- 


transformation defined by 
ponents ol the metric tensor in this case are given by 


b, lh | Tu"), b, b ; 0, and Oss = [b**] a (25) 


In the v', v° coordinate system the condition that the Riemann tensor vanish is satisfied 


if and only if the following differential equation is satisfied. 


2d*b,,/dv?- 1/b,,{db,,/dv* |. (26) 
The two solutions to this differential equation are b,, constant and b,, (Co? + DI) 
Making a coordinate transformation of the form w' av + B,,w a,v + B,, Wwe 
can find a, and £; such that the components of the corresponding metric tensor are 
b,, 1, b 1, and b (w")”, b. 1 respectively for b,, constant and b,, 
[(Co* + D)°. When b |, the w', w coordinate system is Cartesian and when b,, 
[w’}*, it is polar where the curves w’ constant are straight lines in both cases. Since 
each coordinate transformation made in this section has been of the form ¢' t'(s'), 
t° i“(s°), and the curves w' constant are either parallel or concurrent straight lines, 


it follows that the curves of constant speed are either parallel or concurrent straight lines. 


5. Existence conditions. 
motion that the speed is not overdetermined when the two scalars which occur in Rivlin’s 


Ericksen has shown for the case of steady rectilinear 


equations are proportional. For this case we show that there are additional restrictions 
which must be imposed on these scalars in order that steady rectilinear motion be 
possible in every closed cylinder, when both the speed and stress are well behaved 
throughout the region of motion. Specifically, we require in every case that the rate of 
work fid' be non-negative; that the stress exist and be finite in the region of motion; 
that the speed of particles on the cylinder wall be zero; and that the speed be continuously 
differentiable in the region of motion. 

For the case of steady rectilinear motion, parallel to the 2° axis of a rectangular 


coordinate system where the speed is not overdetermined, Rivlin’s equation for the 
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stress 18 
fo o fil 
I tas II pl 8 | + PACD 0 0 Sal 
It. fa 0 | _ 
| i: tihe o | 
+eF (11) I fafo fs 0 |. 


| 0 0 fi tse 

As an immediate consequence of the assumed behavior of the stress and speed, it follows 
from Iq. (27) that /’, (11) is finite at all points where —II = f’, + f?, does not vanish 

We now consider the behavior of /,(II) at those points in the region of motion 
where both partial derivatives of { vanish; namely, when II vanishes. That such points 
exist is guaranteed by the assumptions that f is continuously differentiable throughout 
the region of motion and is zero on the cylinder wall and by the theorem which states 
that any continuous function which is constant on the boundary of its region of defi- 
nition is either constant or has both a maximum and a minimum in this region. If we 
exclude the case where f is constant and make use of the requirement that the stress 
exist and be finite throughout the region of motion, then it follows from Eq. (27) that 
lim,;.0 F,(ID)f., and limy;. F,(11)f.. exist and are finite. When this is the case, it 


follows that 


{flim F\(1Df.,]? + flim F,(IDf,2]7}'"? = C (a constant). (28) 
11-0 Ilo 
IXqg. (28) can be written in the equivalent form 
lim (—IT)'?’F,(11) = C. (29) 
11-0 
When C is zero, then lim,,;.0 F, (IDS, 0 for i = 1, 2 when C + O, then lim,;.o 
F,UD/f.; constant for i = 1, 2. Multiplying the expression F,(II)f,; by (—II'”’) 
{ IT) - we obtain. 
F,(11)(—1D'”-f,(-—1ID7”. (30) 


Making use of Eq. (29) and taking the limit as II goes to zero, we obtain 


lim f.(—II)™'” = constant. (31) 
II-0 
xq. (31) represents the requirement that the components of the unit normal to the 
curves of constant speed must exist at points where II vanishes. An indication that this 
condition ean not be satisfied for C ¥ 0 is given in the ease of radially symmetric Poiseuille 
flow [5]. In polar coordinates, the differential equation determining the speed, f, is 


df /drF({df /dr?) = —Dr (32) 


where D is a constant. By symmetry, and since f is continuously differentiable, df/dr 
clearly vanishes at r = 0. In this case, —II = (df/dr)’ so that Eq. (32) cannot be 
satisfied unless lim ,,.9 F,(II) (—II)'” 0. When this limit is zero, Eq. (32) can be 
satisfied. For example, when F({df/dr]’) is given by K(df/dv)~'” a solution to Eq 
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(32) is 
f=Ar+B (33) 


where A and B are constants. 

The reader should note here that this result does not justify requiring that the 
constant in Eq. (29) be zero, since it is assumed that the flow is radially symmetric, 
and we do not show that this is the only type of Poiseuille flow. 

In section 2 of this paper, we showed that few types of steady rectilinear motion 


were possible when F,(II) is of the form K(—II)~'’*, in which case the constant, C, 
defined in Eq. (29) is non-zero. Accordingly, if we expect steady rectilinear flow to be 
possible in every closed smooth cylinder, when the stress and speed behave as assumed 


in this section, then the scalar F,(II) must be expressible in the form 
F,(1I) = K(-—ID7'? + Y(ID (34) 


where K is a non-zere constant and Y is a function for which lim,,.o (—ID)'” Y(II) = 0, 
and is finite at all points where II does not vanish. 

6. Conclusion. If there are classes of Rivlin fluids for which steady rectilinear 
motion is possible in every closed smooth cylinder when the stress and speed satisfy 
the assumptions listed in the previous section, then these materials must satisfy certain 
definite conditions. These conditions are: (1) the stress be affected by uniform changes 


in the speed; and (2) the scalar 7; (II) must be expressible in the form A II)? + 
Y(II) where A is a non-negative constant and Y(IT) is a function which is finite when 
II + 0 and for which limy; 4» ( I1)'”’ Y(IT) 0. 
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ENERGY DISSIPATION AND LINEAR STABILITY* 


BY 
AUREL WINTNER 
The Johns Hopkins University 


1. Let A be an n by n matrix, A* its Hermitian adjoint, and A, the matrix 


(A + A*)/2. In the real field, this means that A = A, + S, where A, is (real and) 
symmetric and S is skew-symmetric. Hence, if C(x, y) denotes the (real) bilinear form 
belonging to ( A, A, , ---), then, since S(x, x) vanishes identically, A(x, x) = A, (x, 2) 


holds for every vector x with real components. The restriction to the real field is made 
only in order to simplify the notations (and, incidentally, it does not involve any re- 
striction at all, since n can be replaced by 2n). 

If A is a continuous function A(/), where 4, S t < o, and if x(¢) is any solution 


vector ol 


z’ = A(bz (1) 


dx/dt) which is distinct from the trivial solution x(4) = 0, then, as is well-known, 
0 holds for every £. The components of x(/) can be assumed to be real, since A (¢) 
is supposed to be real, hence the real and imaginary parts of a complex solution x(¢) of 
1) are solution vectors. Let (1) be called stable if | x(é)| < const., as t — ©, holds for 


the length of every solution of (1) (and for some constant, which depends on the inte- 
eration constants determining the solution). On the other hand, let 
2” = F(Oz (2) 
be called stable if not only | «(f)| < const. but also’ | 2’(é)| < Const. 
holds for all solutions of (2). 

2. Tor a fixed ¢, let A \(t) denote the least, and uw» = u(t) the greatest, of the n 
eigenvalues of the symmetric part’ A, A,() of the coefficient matrix of (1). Then 
the (unilater: restriction 

aT 
—o < lim sup | u(t) di << @ (3) 
t ed W¥to 
is sufficient r that (1) be stable. With 
r min A(y, y), uw = max A(y, y) (4) 
=1 lyl=1 
where ¢ is fixed in A = A(t) when the min and the max are formed). I found this criterion 


as a by-product of considerations on non-linear generalizations of (1). 


quite a while ago 

*Received Aug. 17, 1956. 

1Concerning various types of linear stability in the scalar case, cf. A. Wintner, Quart. Appl. Math. 
13 (1955) 192-195. 

The point is that only this, the symmetric, part of the coefficient matrix of (1) enters into the 
stability criterion (3). For a manifestation of the same situation in a different result, cf. A. Wintner, 
Amer. Journ. Math. 71 (1949) 362-366. 

2A. Wintner, Amer. Journ. Math. 68 (1946) 557-559. The result and the method of proof were 
rediscovered by T. Wazewski, Studia Math. 10 (1948) 48-59, and by others. Cf. R. Conti, Rivista di 
Mat. Univ. Parma 6 (1955) 3-35, where also further references will be found. 
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Since the proof is very short, and since certain variants of it will be needed when 
dealing with (2) below, let it be repeated here: On the one hand, (4) is clear from the 
identity A(y, y) A,(y, y) (which is equivalent to the definition of A,). On the other 
hand, scalar multiplication of (1) by 2(¢) shows that (| x(4|*)’ 
r(t). But if the trivial solution «(/) # 0 is excluded, then division 


2 is identical with A(x, x), 


where A A(t), x 
by | a(0) ? is allowed, and so (4) shows that 


A(t) S (log | x() |)’ S wd. (5) 
Finally, if (3) is assumed, then it is clear from (5) that log | .(¢)|, hence | ..(4|, must stay 
bounded as t —> 
The simplest instance of (3) is u(Z) S 0, that is, the assumption that the symmetric 
part A,(¢) of A(é) is non-positive definite at every ¢. But more than the stability of (1) 
can then be asserted.* In fact, every | x(¢)| must then have a finite limit (2 0 )ast— @, 
since (5) shows that log | 2(f)| is a non-increasing function if u(t) S O throughout, and 
so the assertion follows from the existence of the lower bound 0 for the length | x(é)|. 


This deduction is mentioned here because a more elaborate form of it will be needed 


below. 
3. The following considerations will replace (1) by (2). For the sake of concreteness, 
it will be assumed that the matrix F F(t) is symmetric (that is, that / Fr, for 


every /), so that (2) is a (non-conservative, but reversible) linear Lagrangian system, 
with n degrees of freedom, having the kinetic and potential energies | x’ |°/2 and — F(x, x) 


respectively, and the latter contains ¢ explicitly (there is, of course, no loss of generality 


in the normalization | x’ |°/2 which, if all masses are positive and constant, puts all of 
them equal to 1). The assumptions will be placed on the (time-dependent) eigenvalues 
of F ’, . Actually, it will be clear from the proofs that the symmetry restriction 
F i ean be omitted if the assumptions are placed on the svinmetric part, F(t), of 
the coefficient matrix of (2). 
The energy along a path x r(t) (which actually must be thought of in phase 
space, rather than in configuration space) is the scalar h h(t) defined by 
Qh rad — F(x, 2). (6) 
But h(i const. does not hold (unless /'(¢) Const.). There is’ however a simple 


explicit rule for the dissipation of the energy which, when subject to appropriate assump- 
tions, can be put to about the same use as the conservation of the energy or its “virial” 
formulation in the classical results (Dirichlet, Jacobi, Hill, Liapounoff) dealing with 
stability in the conservative case. In the present case, the explicit form of the law of 
dissipation is particularly simple’; it states that, along any solution path x x(t) of 
(2), the energy h h(t) must vary so as to render the relation 


dh D(x, x), where D dF and D D(t), 


an identity in ¢. 

‘In this regard, ef. A. Wintner, Quart. Appl. Math. 8 (1950) 102-104. 

5For the general (not necessarily linear) case, cf. A. Wintner, Phil. Mag. 39 (1948) 
6Cf. the reference in the preceding tootnote 

If F(t) is just continuous, then (7), as it stands, is meaningless, but it becomes meaningful if it is 


(22-728. 


meant in terms of Stieltjes integrations, with reference to the ‘‘measures”’ dh(t) or D(t). For, on the one 
hand, a Riemann-Stieltjes integral | f(t)dg(t) exists whenever g(t) is continuous and f(t) is continuously 


differentiable (hence, of bounded variation) and, on the other hand, not only f(t) = x(¢) but also f(t) = 
x’ (t) is continuously differentiable if x(t) is any solution of (2), if F(é) is just continuous. 
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If the coefficient matrix F of (2), instead of being just continuous, has a continuous 
lerivative (for the sake of shortness, this will be assumed from now on), then (7) can be 


vritten in tl 

h’ F'(x,x), where h h(t), ’ = Fb (7 bis) 
the existence of a continuous derivative h’ is part of the statement). In fact, scalar 
multiplication of (2) by x’ leads from (6) to (7 bis), since x’-2’” = (| 2’ |°)’/2. 

4. The f ving fact is a corollary: If the quadratic form belonging to the (symmetric) 
matrix I’ (t n-positive [non-negative] definite at a fixed ¢, then the energy h(t) of 
every solutir tor x a(t) of (2) is a non-increasing [non-descreasing| function at 
th 

In fact, if a a(t) denotes the least, and 8 = B(t) the greatest, eigenvalue of F’(t) 

fixed ¢, 1 7 bis) shows that 
a(t) | a(t) |? = W’() Ss BY) | 2d |’. (8) 
But (8) contains the assertion, since the assumption is that 8(4) Ss 0 [or a(f) = O}. 

Corres] ling to the notation in Section 1, let \ = A(/) be the least, and wp = yp(2) 

the great eigenvalue of F F(t) itself; so that (even if x is not a solution vector) 
A(t) | x |* S F(x, x) S w(Y |x|, where F = F(d). (9) 

Since u(f) < O for A() > O] means that F(¢) is negative [positive] definite, an adaptation 

of the proof of criterion (3) for the stability of (1) will now lead- to the following criterion 

for the stability of (2): 

i) Lf the coefficient matrix F(t) of (2) ts strictly negative definite for large t, and if F’(t) 
posilive definite at every f ced large t, then (2) is stable. 

It is understood that by the first assumption of (i) is meant the existence of a positive 

tant ¢ satisfying 
max F(y, y) wp = pt) S —e < 0. (10) 


lvl=1 


On the other hand, the second assumption of (i) means that 8(4) S 0 holds throughout 


and so, since (8) implies that h’(4) < 0, the energy A(d) tends either to a finite limit h(@ ) 
or to . But (10) prevents the second contingency and, what is more, ensures that 
h( 0. In fact, since (6) and (9) imply that 2h = | 2’ |’ — pl 2° |, it follows from 
10) that 

Qe"h(t) = | x(t) |? + [ 2’ |’, (11) 
if, without loss of generality, € is chosen to be less than 1. But (11) shows that h(d) cannot 


become negative. Hence there exists a finite h(~) 2 0. It follows therefore from (11) 
that z(t) and x’(/) stay bounded as ¢ —» ©, This proves (i). 
5. Suppose, for a moment, that (2) is of a single degree of freedom. Then (2) can be 
written in the form 
zr’ +w()r = 0 (12) 


if —w F < 0 for all ¢, as in (i). Suppose that 


wlh+ao as tao, (13) 
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Then one would expect that the successive waves of every solution a(t) of (12) (waves 
the widths of which must tend, ast — @, to 0 by virtue of Sturm’s comparison theorem) 
will compel the graph of x(t) to be such as to satisfy 


aij)—-0 as to om. (14) 


But it is known that this turns out to be wrong, since even’ lim sup | x(d)| © is com- 
patible with (13), and that (14) follows from (13) only under supplementary conditions 
on the regularity of the growth of w(t), the simplest such condition being the monotony 
of w(t), 


duw(t) = 0. (15) 


The known proofs’ of the sufficiency of (13) and (15) for (14) depend very much on 
Sturm’s separation theorem and cannot, therefore, be adjusted to an extension to the 
non-scalar case of (2), the case of more than a single degree of freedom. It will now be 
shown that the more primitive approach, used above, can readily be modified so as to 
fill in this gap, as follows: 

(ii) If the greatest of the n eigenvalues of F(t) lends to rn ast— o~, and if F(t) is 

) 


non-positive definite at every fixed t, then (14 holds for every solution vector of (2 


Needless to say, the assumptions of (ii) reduce to (13) and (15) in the scalar case 


(n 1). It is also clear that the assumptions of (ii) imply those of (i). But it was seen in 
the proof (i) that there exists a finite limit h( » iO); 

The first assumption of (ii) is that, if « is an arbitrarily small positive number, and 
if tis large enough, say t > ¢, , then —y(t) < 1/e. It follows therefore from (9) and (6) 
that 

2h(t) = | a(t) + | x(t) |"/e if t. K<t< om, 
Since f - 0, and since there exists a finite h(~) = O, this implies that 
Qeh(o lim sup | 2(2) [°. (16) 

But the constant / and the expression on the right of (16) are independent of e’ 
and « can be chosen arbitrarily small. Hence the assertion, (14), of (ii) follows from (16 


6. The following variant of (i) will now be proved 
If F(2) is strictly negative definite, and if the n eigenvalues of F’(t) (or, what is the 
same tl ing, if both functions 


min F’(y, y), max F'’(y, y (17 


luwl=1 y ! 


of t) are absolutely integrable, then (2) is stable. 
Note that (iii), in contrast to the unilateral assumption of (1), does not assume that 
the (internal) dissipating action is steadily decreasing. In fact, the absolute integrability 





For an explicit example, cf. A. Wintner, Journ. Appl. Phys. 18 (1947) 941-942, where a simple 
proof for (14) under the supplementary assumption (15) is also given. 

The restriction (15) is a particular case of a more inclusive condition, formulated by G. Armellini 
(1935) and proved by L. Tonelli and G. Sansone (1936), a condition (on the “regularity of the growth” 
of the coefficient function) which, when added to (13), is sufficient to ensure (14) for all solutions of 
(12); ef. the detailed presentation in Sansone’s Equazioni differenziali nel campo reale, 2nd ed. (1949), 
vol. 2, pp. 52-67, where (on p. 53 and p. 61) further references will be found. 


8Cf. the pre ceding footnote. 
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of both functions (17) of ¢ is compatible with both of them coming arbitrarily close to 
—o and ~ aslt— ~, 

Actually, the last criterion can greatly be improved, since the assumption placed on 
the first of the two functions (17) can be omitted entirely, and the negative pari of the 
second of the functions (17) does not matter either (this agrees with (i), where that 
negative part is arbitrary but the positive part is missing). In other words, all that is 


needed is 


~ 


| y(t) dt < @, (18) 


where, if 6(/) has the same meaning as in (8), 
a max (0, 8) (19) 


(for a fixed 1); so that the situation is as follows: 

(i bis) Jf F(2) ts strictly negative definite, and if B(t), the greatest eigenvalue of F’(t), 
salisfies the unilateral restriction (18), where y(t) is defined by (19), then (2) ts stable. 

The proof of this extension, (i bis), of (i) requires a combination of the proof of (i) 
with an adaptation of a device used in another context”. First, it is seen from (19) and 
(8) that, if a /-value is called of the first or of the second kind according as h’(t) S$ 0 or 
h’(t) > 0, then | A’(O)! S g(| x(d|? for t-values of the second kind. But since F(é) is 
assumed to be strictly negative definite, it follows, as in the proof of (i), that (11) holds 
for every ¢ and for a fixed « > 0. Hence, for t-values of the second kind, | h’(t)| is majorized 
by a constant multiple of y()Ah(d. But (11) also shows that A(2) is positive for all ¢ if the 
trivial solution, x(t) 0, of (2) is excluded. Hence, the preceding inequality, when 
combined with the assumption (18), shows that 


[ | d log h(t) | < (20) 


where denotes the sequence of all /-intervals which consist of the ¢-values of the second 
kind. 

It is clear from (20), where h(t) > 0, that log h(t) tends to a finite limit, hence h(?) 
to a finite positive limit, if 4 tends to © so as to be confined to the set, 2, of the ¢-values 
of the second kind. On the other hand, if ¢ tends to © on the set of the ¢-values of the first 
kind, then h(t) < 0, hence h(t) must then tend either to a finite limit or to — ©. But the 
second case is impossible, since h(t) > 0, by (11). Accordingly, h(t) stays bounded 
whether ¢ tends to © on ¥ or on the complement of 2. In view of (11), this proves that 
z(t) and z’(t) are bounded as t > , as claimed by (i bis). 

7. It is seen from (19) that the criterion (18) for the stability of (2) can be thought 
of as the analogue of the criterion (3) for the stability of (1) (except that, for (2), the 
strict positive definiteness of —/'(t) is assumed in (i bis), in order to complement the 
strict positive definiteness of the kinetic term | x’|*, of (6), the Hamiltonian of (2)). 

The traditional connection between (1) and (2) starts with (2), and, if n is the number 


of components in (2), it replaces n, x and F by 2n, (a, x’) and 


iul™ * (21) 
FO 


°For a slightly weaker version, cf. the reference’. 
0A. Wintner, Amer. Journ. Math. 69 (1947) 5-13. 
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respectively (J denotes the n-rowed unit matrix). But there is also another formal con- 
nection b ( 1 (J nd +), a connectlo vhich starts with (1) 


Professor W. Prager was good enough to call my attention to a paper which he 


published in a not ly accessible peri \dical'’ and in which, starting with a scalar 
equat ! ne tori / he is le | by d rentiation of (3 to the consideration of scalar 
equations (of third o1 irth order in a hich are of second order in 2’ or x’’. An inter- 
pretation of this pt in the phase space of (2) suggests a corresponding procedure 


for (1 
The . th I traighforward enough, 1 quite inst uctive. In fact, it turns out 


that what leads fro | 2) (for any n precisely the (matrix) equation of Ricca 
F(t) 1’(é) + A*(é). (22) 


In fact, if A(t) is d rentiable, then (1) implies that 2” 1’x + Az’, and so (2) follows 
from (1) if F is defined by (22) 

Needless to say, the case (22) of (2) is only a necessary condition for a solution x(t) 
of (1), mply becau the general soluti of (1) de pends on 7n, whereas that of (2) 
depends on 2n, int on constants. The pitfalls of this situation can be illustrated as 
follows 

Let F(t) be any ? n matrix which is continuous for large positive t and has the 


; { $ xy] 


property that, whe fixed, its syinm 


Tric part 


F,() is non-negative definite 23 
(for inst: , that / svmmetrie and non-negative). Then it is known” that, if x(¢ 
is any solution vé f (2), the graph of | x(t)|* in the (¢,s)-plane always turns its 
convexity toward tl <is and that. if z,(t), --+ , ,(f), asi (tl), «++ , 2, (2) is a complete 
set of linearly independent solutions of (2), then it is possible to choose x(t), 
to be bounded vectors (hence, by convexity, such that (AI, + 5 | Sell tend to 
finite, non-negative ts ast — o), whereas 2,,,(¢), r.,(t) are unbounded (which 
implies that holds for any solution vector which is linearly independent of 
Zs 


2 } ») tial) 


Let tl e app! 1 to the particular ease (23). It then follows that, if (22) sa 
(23), the ase (22) ol 2) has exactly n, il early independent, solution vectors which are 
bounded as t > _ But none of these n solutions will in general be a solution of (1) 
also. In fact, if A(¢ T for all ¢, then (22) satisfies (23), but (1) reduces to x’ = x and 
has. therefore no bounded solution distinct from x(t) = 0 

uW. Prager, Rev. de la Fac. des Sciences de |’Univ. d’Istanbul, 1 (1936) 37-43. 


12See the paper referred to under 3, 
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4 


Summary. A somewhat limited design procedure for elastic-perfectly plastic struc- 
tures was developed previously [1]**. It is extended here to provide upper and lower 
bounds on the minimum weight of three dimensional structures and is specialized to 
safe one and two dimensional structures in which either direct stresses or bending stresses 
are negligible. The generalization also includes the influence of body forces. In principle, 
therefore, such troublesome factors as the weight of the structure itself or centrifugal 
“forces” may be designed for in a direct manner. Radially symmetric plane stress and 
plate bending examples are solved to demonstrate direct design procedures. 

Introduction. Although design rather than analysis is the real problem in machines 
and structures, far more research effort is spent on analysis. The reason is primarily 
the specific nature of the problem posed and thé greater possibility of obtaining an 
unambiguous solution. In what follows the material of construction is idealized as 
elastic-perfectly plastic. This first approximation to the behavior of real structural 
metals beyond the elastic range and the accompanying techniques of limit analysis are 
suitable for problems in which load carrying capacity is of primary importance. In 
addition to providing more realistic answers for ductile materials, the great advantage 
of limit analysis over elastic analysis is its relative simplicity. The assumption of perfect 
plasticity opens the possibility of direct design in the strict sense as opposed to pre- 
liminary guess and repeated analysis with the end point not necessarily the best that 
ean be achieved. lor problems of such difficulty that direct plastic design is not feasible, 
it is possible to obtain bounds on an optimum design. Such bounds provide a sound 
basis of comparison with any proposed structure and indicate the gains, if any, which 
ean be achieved by further refinement. 

An economical structure ordinarily will be one of least weight within the restrictions 
of the procedures of fabrication. I’'ramed steel structures, for example, often are best 
constructed of beams of constant cross-section between joints. Heyman [2], Foulkes 
[3], Prager [4] and Livesley [5] have studied minimum weight frames of this type. The 
present paper dealing with beams, plates, sheets, and space structures is in the spirit 
of Michell [6] where the minimum weight is sought without regard to problems and 
costs of manufacture and construction. Prior studies of plates in bending from this 
point of view have been made by Hopkins and Prager [7] [8], Freiberger and Tekinalp 
[9] 

A start on a broad theory was made [1] by the establishment of a criterion for absolute 
minimum weight design for structures which are subjected to direcé or membrane 
stresses and for relative minimum weight in the case of beams or plates in transverse 
bending. These criteria which include that of Michell as a special case suffer from the 
same disadvantage. They cannot be satisfied in all, and probably not in most, problems 
which arise. As part of the extension of this earlier work, the matter of existence of 

*Received Oct. 4, 1956. The results presented in this paper were obtained in the course of research 
conducted under Contract Nonr 562(10) by the Office of Naval Research and Brown University. 

**Numbers in square brackets refer to the bibliography at the end of the paper. 
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solutions will be explored briefly, but major emphasis will be on the inclusion of body 
force and the obtaining of bounds. 

2. The problem and its relation to limit analysis. The theorems of limit analysis 
[10] are concerned with whether a body of elastic-perfectly plastic material will collapse 
or not under given loads. The term collapse is used here to describe conditions for which 
plastic flow would occur under constant loads if the accompanying changes in the 
geometry of the body were disregarded. The boundary conditions are assumed to be 
of the stress type. That is, each of the components 7’; of the surface traction is prescribed 
except where the corresponding velocity component is prescribed to be zero. The behavior 
of the elastic-perfectly plastic material is characterized by a positive yield function f 


which is a function of the stress components o;; (2, j 1, 2, 3) referred to rectangular 
Cartesian coordinates x; . The material behaves elastically under states of stress for 
which f < k’, where k is a constant. Plastic flow can occur under states of stress for 
which f = k”. States of stress for which f > k° are not permissible. 

The yield condition /(c;;) k* can be considered to define a surface in stress space 


in which the components of stress are used as coordinates. Plastic flow can occur under 
a state of stress represented by a point on this surface. In the same space, a plastic 
strain ra/e e?; can be represented by a vector with components proportional to the 
components of the plastic strain rate. For a perfectly plastic material the vector represent- 


ing the plastic strain rate ¢?; associated with a plastic state of stress o;,; is in the direction 


of the exterior normal to the yield surface at the stress point o,; . If the yield surface 
has a singular point where there is not a unique normal, then the strain rate vector 
must lie in the fan bounded by the normals to the yield surface at adjacent points [11]. 


The rate D of dissipation of energy per unit volume of the material due to plastic 


action is given by a;,;<« For a given yield function /, the plastic flow rule enables the 
dissipation rate D to be expressed as a function of the plastic strain rate only, 

D 0 5; €;; D(é;). (1) 

The term safe statically admissible stress field will be used to denote a stress field 


a;; Which satisfies the equations of equilibrium 


Oi; + F, 0 (2) 


7 


throughout the volume V of the body, satisfies the boundary conditions 


o,n; = T (3) 


‘ 


on the surface S where the tractions 7’; are prescribed, and which is below yield every- 


where in V, f(o;;) hk". In (2), F,; are the components of the body force per unit volume 
of the material and in (3), n; are the components of the unit normal to the surface S. 
Repeated subscripts in the same term denote summation 

A velocity field u,; will be called kinematically admissible if it represents a possible 
plastic flow in the body and if it vanishes on those portions of the surface of the body 
where the surface traction is not prescribed 

The following theorems have been formulated [10] for an elastic-perfectly plastic 


body under surface tractions 7’, 


Theorem 1. If a safe statically admissible stress field o;, can be found then the body 
will not collapse under the loads Ts 
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Theorem 2. Collapse must impend or have taken place previously if a kinematically 
admissible velocity field u* can be found such that the rate at which work is done by the 
external forces equals or exceeds the dissipation rate, that is for u* continuous in V 


[ veévav s [ ratas+ [ ratav. (4) 
V 8s Vv 
In (4), the rate of dissipation D is calculated from the strain rates 


é; —_ Aut; + M5 .a)> (5) 


which are treated as purely plastic. 

As the basic problem of design for minimum weight we consider a volume Vy within 
which a structure loaded by a given force distribution is to be placed. The loads 7’; are 
to be supported by some distribution of a given elastic-perfectly plastic material in the 
volume V, . The problem is to determine a structure or distribution of the material 
which is just at the point of collapse under the loads 7’; and which involves minimum 
volume of the material. We shall assume that the material is homogeneous so that the 
problem of design for minimum volume coincides with the problem of design for minimum 
weight. 

It is apparent that for the problem to admit a solution it must be possible to design 
at least one structure in Vy capable of supporting the loads 7’, . Also the loads 7’; will be 
reached by some loading program and it is presumed that the minimum volume structure 
would not collapse at any previous stage of the loading program. 

We remark that in design for minimum weight, modifications to the approach below 
can be made to allow for non-homogeneity if it is prescribed that different materials are 
to be used in different specified regions of the volume Vr . 

3. Bounds on the minimum volume. A design based on any safe statically admissible 


stress field o°, in Vy for the given material will support the loads 7’; . The body formed 
by filling the volume Vs where o;; is non-zero with the given material will not collapse 


under the loads 7’; , by Theorem 1 above. In general the stresses a7; are not the actual 
stresses in the body Vs subject to the loads 7; . By definition, the magnitude of the 
volume V,,, of a design with minimum volume is bounded from above by the magnitude 
of Vs, Vn < Vs. In practice it is convenient to use stress fields for which f(o?,) < k’. 
lor such fields V,, < V, . 

In order to bound V,, from below the second theorem of limit analysis is applied. 
Since the minimal design is just at collapse under the loads 7’; , then by Theorem 2 for 
any kinematically admissible continuous velocity field ui defined in Vr , 


[ D(é:;) dV > [ Tu dS + [ Fuk av, (6) 
Sei Ja To 


where A is the loaded surface. The body forces /’,; are assumed to act only if material 


is present. The equality sign holds in (6) only if u‘ is a possible collapse mode of the 


: 


minimum volume design. Inequality (6) can be written 


/ A(u') dV > i) Tut dS, (7) 
Vea A 


where 


A(ut) = D(é;) — Fyuj . (8) 
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A lower bound on V,, follows from (7 
V,>Vi= | Tui dS/max A), (9) 


where max A (u 3; the maximum value of A in V7 which includes V,, . We note that 
A must remain finite ’, if (9) is to provide a non-zero bound. For a material obeying 
° 7 


either the Tresca o1 Mises yield condition this requires the velocity components 
h 


u; to be continuous throughout V,7 and to satisfy the condition of incompressibility. 


For some materials discontinuous velocity fields will be permissible if the material is 
such that the work « per unit area of the discontinuity surface is zero. Soil which 
is unable to take tension provides such an example when the discontinuity in velocity 


is associated with the formation of a tension crack [12]. 


The lower bounds V, determined from kinematically admissible velocity fields wu; 


00 low to be of val 


by (9) may prove to be ie in bracketing V,, . The less restrictive state- 
ment (7) may then provide useful information by employing the following procedure. 
For a given kinematica \ admissible velocity field u; in Vz, the value of the modified 
dissipation function A (u;) will in general vary throughout V,; . Material is placed in 
the regions of V, where A takes its largest values until a stage is reached at which 
material occupies a region V, of V7 and 


. . 


| A(ui) dV = | Twi ds. (10 


vV; 7A 


The surface bounding V, will be a surface on which A has a constant value. The values 
of A at points in the volume V, will not be less than the values A assumes at points in 
V, not in V, . It follows that the volume V,, satisfying the inequality (7) cannot have 
less volume than } V,, > V, . Equality between the lower bound V, and the lower 
bound V, given by (9) can only occur if A (u;) is constant in the region V, so that in 
general the bound V, is an improvement on the bound JV, . 

4. Planestress. The case of plane stress follows the general approach of the prev ious 

I 


section. It is required to design a plane sheet of variable thickness which is just at the 


e 1 


point of collapse under loads acting in the plane of the sheet and which has minimum 


volume. The sheet is formed by distributing a given material at points on or at a small 
distance from a p e middle surface A The variable thickness of the sheet at any 


point of A; is denoted by h. Forces 7h which are independent of h act in the plane of 





A, on specified lines L in the plane A, . The body forces F; per unit volume of the 1 


also act in the plane of A 
The middle plane A ; is taken to lie in the plane x 0 for definiteness. The assump- 
tion of pl ne stress requires the stress col ponents « to vanish. The remaining com- 


ponents o a, — |, 2) are functions o and T2 ONL and for equilibri im satisfy the 
equations 


(o.sh) gs + Fh = 0. (1] 


The stress boundary conditions on the lines Z where the forces Th are prescribed 
require 


(og ofl )7 — 7 on (12) 


where n, are the components of the unit normal to L in the plane x, = 0. 


ae 
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A kinematically admissible velocity field u; representing a possible plastic flow of a 
sheet with middle surface A, is such that u, and wu, are functions of x, and x, only. The 
rate D of dissipation of energy per unit volume can be found from the velocity components 


1, , U. since the stresses o;, are zero, 


D = castap = D(€as)- (13) 


Tap 


As in Sec. 3, a design based on any statically admissible stress field o%,, will provide 
an upper bound for the minimum volume V,, . The stress field o4,, must satisfy Eqs. 


11) and (12) for some distribution of thickness h, and must be at or below 
yield, f(o° < k*’. The volume J, , 
V.= | h dA, (14) 
JAT 
} uppel bound on V. 


or any kinematically admissible velocity field defined by the components uf , uj 
over .1, it follows from Theorem 2 that since the minimal design of thickness h,, is just 
at the point of collapse under the loads 


. 


[ Auth dA > [ (Pahut al, (15) 
Jar Ju 
where 
A(u*) = Deis) — Fauée . (16) 
The lower bound V, on V,, follows from (15) as-before, 
V2 ¥,.% [ (Thus dL/max A(uS), (17) 
Ju 


where max A (u2) is the maximum value of A in the plane area A; . 
As described in Sec. 3, the inequality (7) can be used to improve the lower bound 
(9) and similar remarks apply here with respect to inequality (15) and the lower bound 


(17). However in this case an upper limit may have to be placed upon the thickness h 
in order to avoid the somewhat unrealistic line flanges or similar reinforcements of 
finite cross-sectional area and therefore infinite height. 


2 


5. Membranes. In analogous fashion to the case of plane stress, the theory of Sec. 3 
can be applied to membrane design. The membrane has a prescribed curved middle 
surface A, and supports prescribed loads which in this case may act normal to the 
surface A, as well as in the tangent planes to 4; . 

6. Special stress and velocity fields. A special case arises if it is possible to design 
a body V. in Vy, which is compatible with a collapse mode u‘ such that A (u{) is constant 
throughout V, which includes V. . In the absence of body force this is the problem 
treated in the previous paper [1]. As uf is the actual collapse mode for V, under the 
loads T; , 


. - 


A(us) dV = | Tui dS. (18) 


For any other design V, in Vy just at the point of collapse or capable of carrying the 


© 


loads T,; , by Theorem 2 above, 
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[ A(ui) dV > | T wu; dS. (19) 
¥V, 7A 
As A (u‘) is constant throughout V; , comparison of (18) and (19) shows that V. < V,. 
It follows that a design compatible with such a collapse state is a minimal design V r 
The term compatible means that there exists an equilibrium state of stress which is 
consistent with the plastic strain rates and which satisfies the boundary conditions on 
the surface traction and satisfies the yield condition everywhere in V, . 

A similar result applies to the plane stress and membrane cases. I’or plane stress, 
a sheet design V, of thickness h, compatible with a collapse state for which A given by 
(16) is constant over the middle surface A, is a design with minimum volume. An 
example which uses this result is given later in Sec. 8. 

In most eases it will not be possible to design a body compatible with a collapse 


state A = constant. This may be seen by analogy with small strain linear and non-linear 
elasticity. For a positive yield function f(o;;) which is homogeneous of degree n in the 


stress components o,;; and which defines a smooth yield surface, the plastic strain rates 
are given by 
P of / 
Oe ae (20) 
Oa;, 
where \ is a non-negative scalar function of position. In the absence of body force, the 


a p . : 
dissipation rate A j or De! is given Dy 


D(e;;) 05 5€., navf. (21) 
The body is at yield everywhere, f(o; k*. It follows that for D to be constant, A 
must be constant throughout V. . For n 2, relation (20) is then formally equivalent 


to a linear stress-strain relation between elastic strains e:, and stresses Cis» the function 
Af playing the role of strain-energy density. In general a solution of the equations of 
elasticity will not have constant strain-energy density, and therefore in general a collapse 


> 2 the analogy is with 


state with D constant will not exist for the body V..lorn 
non-linear el isticity and the function Af Is equivalent to the complement of the strain- 
energy density. 

The result V V. also follows from (18) and (19) if A (uj) is constant throughout 
V.and smaller elsewhere in V, . A stronger result is the following. 

A design occupying the volume V, compatible with a collapse mode such that the 
value of A at any point in V, is not less than the value of A elsewhere in V; is a minimum 
volume design. 

A velocity field which is discontinuous across a surface may be considered as the 
limiting case of a continuous velocity field. The approach of Sec. 3 and the result of the 
previous paragraph therefore apply to discontinuous velocity fields. The term involving 
the dissipation of energy at the discontinuity surface can be taken into account ex- 
plicitly as in Ref. 9. Discontinuous velocity fields may prove useful but only if the total 
rate of energy dissipation at the discontinuity surface is finite. 

When the formal approach of this section is of no practical value, the more general 
approach of Sec. 3 must be used to obtain either the minimum design or bounds on its 
volume. 

7. Bending of plates. A thin homogeneous plate in bending considered as a three 
dimensional body provides an example in which A cannot be constant throughout. 
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Here the results of the previous paper [1] are generalized to include body force. In bending, 
the strain rates €,, in the plate vary linearly throughout the thickness and the rate of 
dissipation per unit volume is given by 
a Pa 2 | v3 | 
Die.) = Oo apt ap _ h Dy ’ (22) 
l 

where D, is the value of D at the plate surfaces 7, = + h/2. The modified dissipation 
rate A defined by (8) is therefore given by 


9 7 
A(w) = 2 | #1 De — Fw, (23) 


where w is the rate of (normal) deflection of the plate. For a plate, with middle surface 
A, and thickness h, , which is on the point of collapse in a deflection rate pattern w* 


[ pwd + [ Pyw'h, aA : [ D(é.,) dV = / 1Dth, dA, (24) 
JA J Ae / Ve Ae 


where p is the transverse load per unit area. For any other plate, with middle surface 
1, and thickness h, , just at collapse or capable of carrying the load p, 


[ po dA + [ Fach dA < [ 4Dih, dA. (25) 
JA JAs Ae 
From (22) 
, wh, : 
dD; Di | ; (26) 
s 
and it follows from (24) and (25) that 
[ (GD; — Fwy. dd < [Dit — Pyw h, aA. (27) 
J Ae JA be 


We now suppose that the plate with middle surface A, is any neighboring plate to 
the plate A, in the sense that h, = h, + 65h where 6h « h, . The middle surfaces A, and 


1. are then the same, A, = A, , and if the second order term in (27) is ignored 
[ (Di — Faw )sh aA > 0. (28) 
Jac 
[t follows that if Dj — F,w* is constant (and positive) over A, , 
[ h,dA < [ h, dA. (29) 
Ja Jae 


Thus a design compatible with a deflection rate pattern for which A is constant on the 
surfaces of the plate is a relative minimum in the solution of the minimum volume 
problem. 

It is more usual to consider a plate in two dimensions rather than three. The corre- 
sponding result in terms of bending moments 1/,,[9] and rates of curvature k,, can be 
obtained directly or it can be derived from the work above. The rate of dissipation per 
unit area of the middle surface is given by 


h/2 
Nua I. D dry = }Doh. (30) 
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A relative minimum is provided by a design compatible with a deflection rate pattern 


such that 2M. ,k.3/] F’,w is constant over the plate. 

The results for the non-homogeneous sandwich plate are more definite. By a sand- 
wich plate is meant here a plate which has a light-weight core of constant or variable 
thickness H between thin identical face sheets each of thickness ¢ which varies from 
point to point of the middle surface A 
bending moment across the plate is carried by equal and opposite membrane stresses 
in the outer faces. The problem is to determine the variable thickness ¢ of the faces so 


Che core carries no bending stresses and a 


that the plate is just at the point of collapse under the prescribed loads and the face 


sheets have minimum volume. Since H is considered to be large compared with ¢, H > t, 
the strain rates in the faces of the plate can be considered to be constant throughout 
t. It follows from the general result of Sec. 6 that a sandwich plate designed 


the thickness 1 
to collapse in a mode w such that A is constant in the faces of the plate is an absolute 


n 


minimum volume design 

8. Plane stress example-rotating disc. In the previous paper [1] the formal approach 
of Sec. 6 was used to obtain minimum volume designs for a circular annular dise loaded 
by uniform forces per unit length on the inner and outer edges of the disc. Here an 
illustration is given of the way in which centrifugal “forces” can be taken into account. 

\ circular annular disc of inner radius r; and outer radius ry is rotating with constant 
angular velocity w about a perpendicular axis through the center of the disc. The outer 
edge of the disc is loaded by a uniform tensile force 7’ per unit length and the inner 
edge is stress free. The problem is equivalent to that of a stationary disc with outward 
radial body force pwr per unit volume, where p is the mass density of the material. The 
sumed to obey Tresca’s \ ield condition, lig. l, with yield stress 


9 


material of the disc is as 


a 





Fie. 1 Tre ’s yield condition for plane stress with radial symmetry. 
‘or equilibrium the principal stresses o, and satisfy 
d , J aS . . 
(ho,) — F — pw Th, (31) 


dy r 
where A is the thickness of the disc. In terms of the radial velocity u the strain rates 
are given by 


e. = du/dr, & = ur, (32) 
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and for A constant, 


A o,€, + o€) — pwru = constant. (33) 

Che condition A constant restricts the position of the stress point on the yield 
hexagon of lig. 1. For example, the stress point cannot lie on the side AF because the 
normality condition requires e, 0, and the condition (33) cannot be satisfied in view 


of (32). It is readily found that only the stress states represented by points C and F, 
B and F, can be associated with a velocity field for which A is constant for a finite range 


or fr. 














. 

In the present problem the stresses are tensile so that the stress points B and C are 
prohibited. The stress point /# is also inadmissible as it would require inward radial 
velocity at the inner edge r = r; . The remaining stress point F provides the solution. 
At the stress point F, ¢, = ao» = oo and equilibrium requires 

dh pu” 
— = ——Ph, (34) 
dr Jo 
The region r; < r < 19 of the disc is at the stress state F so that 
T oie : . 
h = — exp (x — 2) <2 SX, (35) 
do 
where x = (pw’/2e))'” r. In order to complete the design, a line flange must be added 
at the inner boundary, lig. 2a, to satisfy the condition of zero stress at r = r,. The line 
flange is of infinite height but of finite cross-sectional area of amount 
2 ’ 
rh,/Q — 225), (36) 
| |>w 
a «flange T 
| V7 host 
A f ; \ o 
_ | WO 
— Tr 4 
| V4 
V4 
om @§ =<] 
- b—— |, —_— 
| i ’ r — 
| fo [" , 
; (a) l (b) 


Fic. 2 Minimum volume design for rotating annular disc under exterior tension (a) of unlimited thickness 
and (b) of restricted thickness. 


where h, is the value of h given by (35) at x = 2; . The speed of rotation of the disc is 
restricted by the condition 2? < 1/2, that is pw’r? < oo. 
It remains to associate a velocity field of constant A with the design. For the stress 
point F, 
du ,u . 
oo\— +—] — pwru = constant (37) 
dr F 
and 


¢, = du/dr > 0, « = u/r> 0. (38) 
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The line flange is at the stress point A so that du/dr < Oatr = r,; . Taking du/dr = 0 
at r = 7; gives with Eq. (37) 
“u= [exp (2° — 25)/(1 — 225) — 1], (39) 
x 
where A is a positive constant. The restrictions (38) are satisfied for r > r,; provided 


that zx? < 1/2. 
The volume V,, of the minimum volume design is given by 


V pw . ' 
sieeall = *" ( +* ‘i, “ > Es — ( 
nT exp (2 x)/U ) 1, (40) 


and is an increasing function of 2, . 
The somewhat unrealistic line flange can be avoided by placing an upper limit H on 


the thickness h, Fig. 2b. The full thickness H is taken for r; < r < r2 and the stress state 
in this region is represented by points on the side AF of the yield hexagon, a, = a , 
0 <a, < o,. The edger r, of the full thickness region is found from the condition 
that (c,h) is continuous at r = r, . The collapse mode must be such that du/dr = 0 for 
r,; <r <r.,80 that u is constant in this region, u B, and 
A B(oo/r — pwr iT <a %. (41) 
Thus A is a decreasing function of r forr; < r < r,. Forr, <r < ro the radial velocity 
is determined by (37), (38) and the condition du/dr = Oatr =r. , 
6 
u — [exp (a° — x2)/(1 — 222) — 1] So Ms (42) 
z 
The positive constant C is chosen so that u is continuous at r = r, . The value of A is 
constant for r > r, and equal to the value given by (41) at r = r, . A lower limit must be 


placed on H and an upper limit on the speed of rotation w for the design to be possible. 
The design is clearly the minimum volume design but this can also be seen from the follow- 
ing. The design is compatible with a collapse mode for which the value of A in the body 
is not less than the value of A elsewhere in the region where material may be placed and 


is therefore a minimum volume design. 


9. Sandwich plate example. The problem of this section illustrates the direct design 
approach of Sec. 6. A simply supported circular sandwich plate of constant core thick- 
ness 7 is loaded by 
is required to find the variable thickness ¢ of the two identical face sheets so that the 
plate is just at the point of collapse under the pressure p and the face sheets have minimum 
volume. The weight per unit volume of the face sheets is denoted by y and the body 


uniform pressure p over its upper sur!.ce, which is horizontal. It 


force is taken into account in the analysis. We first outline the design procedure for a 
material with a general yield function and then consider the particular case of the Tresca 


yield condition. 
We write the yield condition on the radial bending moment M and the circumferential 


bending moment WN in the form 
F(M,N) = M,, (43) 


where M, is given by 


M, = oft (44) 


‘ 
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and og» is the yield stress of the material in tension or compression. The yield function 
F is homogeneous of order two in the moments M and N and will be assumed to define 
a smooth curve in the M, N plane. The von Mises yield criterion for example provides 
the special form F = M* — MN + N*® = Mj. For equilibrium the moments satisfy 
the equation 


d 


F dN 
;(rM) — — r+ 2yir = 0 45 
dr ) dr + 7 + Y ’ ( ) 
where r measures distance from the center of the plate. Measuring the deflection rate 
w in the downward direction, the radial and circumferential rates of curvature «x, \ are 
given by 
dw 


t= “— = 


dr’ 


] dw 


r dr (46) 


The moment-rate of curvature relation associated with the yield condition (43) requires 


dn = ak /38. a 
OM ON 

rom the discussion in See. 7, a design compatible with a collapse mode w such that 
the modified dissipation function A (w) is constant throughout the plate is a minimum 
volume design. The radial and circumferential strain rates in the upper face of the plate 
are } Hx and — 3 HX respectively and are equal in magnitude and opposite in sign 
in the lower face. We therefore seek a deflection rate w such that 


‘4 d?w 
A(w) p( H me : ; a aw) — yw = constant, (48) 


where the form of the dissipation function D is known from the yield condition of the 
material. At the simply supported edge r = r, of the plate the deflection rate is zero 
and at the center dw/dr = 0 since 2 is finite there. Also for equilibrium 1 = N at the 
center of the plate so that « = \ at r = O since the yield condition (43) is symmetrical 
in M and N for an isotropic material. In general, Eq. (48) will be non-linear and analytical 
solution of (48) subject to the conditions on w and its derivatives will be difficult. A 
possible numerical approach is the following. Numerical values are arbitrarily assigned 
to the positive value of w and the common negative value of d’w/dr* and (1/r)dw/dr 
at the center, 


‘d?w 1 dw 
(w),-o a> 0, ( z = |- = —b < 0. (49) 
dr r=0 r dr r=0 


The constant on the right hand side of (48) is then determined by substitution of the 
values (49). With this value of the constant, Eq. (48) is integrated numerically in the 
range 0 <r < r, to give a deflection rate which will be denoted by W. If Wo is the value 


of W atr = r, , then the deflection rate w given by 
w= W— Wo (50) 


satisfies Eq. (48) with a different constant on the right hand side and satisfies all the 
conditions on w and its derivatives. 
When the collapse mode w such that A (w) = constant is known, the rates of curvature 
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and the normality condition (47) 


16), The yield condition (43) 
1 N j 1 terms of the thickness t. The 


Y In terms of \/ and henee i 
en so that MW and N satisfy the equilibrium equation (45). At 
i, U so tis taken 0 be zero atr \lso the resultin: 


iP equation 


ngular at r U and the solution must be such that ¢ is finite at 

| a minimu ime design can be obtained for a mate rial with 
iT) (43 . 

e Tresca n, which does not define a smooth yield 


in consider 


e red here the curvature rates x, \ 
e condition A tant requires 
\ 
H ( ! l dw\ J 
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2 dr/ 


= J a;r, P 
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Creep and f aciure ¢ nelals at high len pe rature. Proceedings of a Symposium held at 
National Physical Laboratory, Teddington, England. Her Majesty’s Stationery 


Office. London, and British Information Services, New York, 1956. iv + 420 pp. 


$5.60 

This be is a summary of the proceedings of an international meeting concerned with the creep 
and fr I sat high temperatures The meeting was held at the National Physical Laboratory, 
Tedding England Mav and June 1954. There are twenty-two papers falling into four categories: 

Section I, Deformation Processes In Simple M iterials; Section II, Creep Resistance of Complex 
Materials: Section III, Theory of Fracture; Section IV, Tertiary Creep and Fracture. 

These proceedings | le an excellent immary of a cross section of the present ideas (on an 
inte rnatior s if t creep and fracture it elev ited tempt ratures In the two years since this 
conferet { e have some significant advances in the understanding of grain boundary motion. 


,0oHN ‘TRUELL 


By Aris Phillips. The Ronald Press Company, New York, 


19 20) } 87.00 

Che i ich form one-| the book, deal with the one-dimensional problem of 
the p! i ) d frames, the plastic bending of beams with and without axial forces, the 
c d frame m the moment-curvature relation, and an elementary 
theo f curved rs. In the preface the author states that the book is written from the point of vie w of 
the s \ pproxl pages Chapters 2-6 ire devoted to the listing 
( hi! l I probl 

I I R in relations are stated in Chapter 7 and are applied in the r naining 
two cl lving combined s. Except for a brief discussion of plane strain, 
the p mplifying f ial or radi metry. The treatment throughout 
is nd is at time hat inelega 


1T>] 


Mall VM } Volumes I : 1 II of Handbuch der Physik). Kdited by 
Ili Spl \ y Berlin-G ne Heidelberg. Vol. I: 1956, vu + obt pp 





817.2 \ 1 | | 520 pp S°0. 05 

The volumes contain the following articles. Volume I: Fundamental Concepts of Classical Analysis, 
Ordi Differential Equations, Theorv of Functions (J. Lense, Technische Iochschule Miinchen, 
89 py Partial Di ntial Equations (J. Lense, 30 pp.)—Elliptic 1 inctions and Integrals (J. Lense, 
27 pp. : 1 I ) ( Mathematical P ies (J. Meixner, Technische Hochschulk Aachen, 
71 pp \ Problems (F. Schligl, Universitat Kéln, 135 pp.); Volume II: Algebra (G 
Falk, T nische Ho Aachen, 116 pp Geometry (H. Tietz, Technische Hochschule Braun- 
schweig, 81 py Funct Analvsis (I. N. Sneddon, University College of North St iffordshire, 151 
pp. Numeri ind UI ical Methods (L. ‘ itz, Universitit Hamburg, 122 pp.) Modern Com- 
puters (H. Biickner, Gen 1 Electric Compar Schenectady, N.Y., 28 pp.). A critical discussion of 
these irticles, or ever mere listing of section he dings, would require more space than is available for 
this review. Since, however, Sneddon’s article, which is in English, may interest a greater number of 


in German, its scope is indicated by the 


readers of this Quarterly than the other articles, which are 
following headings and subheadings: Integration and Abstract Spaces (Introduction—Banach Space)— 
Integral Transforms (The Laplace Transforms—The Fourier Transforms The Mellin Transform 
The Hankel Transform—Finite Transforms—Approximate Methods of Evaluating Transforms)— 
Hilbert S; Schwartz’s Theory of Distributions—Variational Methods in Functional Analysis. 
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ON THE ANALYTIC FUNCTIONS OF ORDER n' 
(APPLICATION TO THE PLANE ELASTICITY) 


BY 
SAFFET SURAY 


Ankara University and Harvard University 


|. In the course of an investigation on a problem of plane elasticity we have been 
led to consider the so-called “analytic functions of order n’’. These functions, introduced 
for the first time by P. Burgatti [2]’, are, by definition, expressions of the form 


n—1 
,(z,2) = > 2" f.(2), (1) 
0 
where f,(z) (v = 0,1, --- ,n — 1) arearbitrary analytic functions of the complex variable 


n indicates the order of the function. The properties of analytic functions of order n 
seem to be in close connection with those of polyharmonie functions of two variables. 
On the other hand, the latter, especially the biharmonic functions occur in many problems 
of physics and mechanics; and therefore, their study is particularly important for appli- 
cations. In what follows, keeping in view this point, we shall establish certain properties 


x +7yina domain D and Z is the imaginary conjugate® of z. The positive integer 


of analytic functions of order n and indicate some applications to two-dimensional 
elasticity. 

An immediate property of the analytic functions of order n is the following. 

The product of any finite number of analytic functions of order n, , m2, +++ 5 Mm 3 


2, , Qn, 5 °** , Qe, defined respectively in domains D, , D, , --- , D,, is an analytic func- 
tion Qy = Q,,,,---Q,, of order N = ny + ny + +++ + 2, — (m — 1) defined 
in the domain D’ which is the common part of D, , D., --- , D, . 


We shall mention two special cases: (i) the product of an analytic function of order 


n, 2, by an analytic function of order 1, Q, , i.e. by an analytie function in the ordinary 
sense /’(z), is an analytic funetion Q* = F(z)Q, of order n; (ii) the mth power of Q, is 
the analytic function (2,)” of order m(n — 1) + 1. 


The above results follow very readily from the form of expression (1). 

2. Connection with real polyharmonic functions. A real polyharmonic function 
of order h, H of the variables x, y is, by definition, a solution of the equation (d?/dx* + 
a°/ay’)" H = 0 or, in terms of the variables z and 2, 0°*H/d2"dz" = 0. Burgatti showed 
that the real and imaginary parts P and Q of an analytic function of order n, 2, satisfy 
the above equations, and called them “conjugate polyharmonic functions’’*. These 


‘Received June 25, 1956. This work was sponsored in part by the U. S. Office of Naval Research 
under Contraet N5ori-07634. 

2Numbers in brackets refer to the bibliography at the end of the paper. 

‘In the following any quantity % will represent the imaginary conjugate of z. 

‘This designation is justified by the fact that Almansi’s representations [1], [2] for two such poly- 
harmonic functions contain two sets of harmonic functions which are conjugate two by two. 
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functions can be written in the following form 


P=19,+ 9) =4 - &Sf2+2’h), 

a (2) 
Q ti, ~ & > eS 2" f(2)] 

a are eer 


Burgatti’s definition is concerned with the conjugate polyharmonic functions of order 
n which are the real and imaginary parts of an analytic function of the same order n, 
2, . However, it should be remarked that the number n, i.e. the order of the analytic 
function ©, , is generally different from the order of harmonicity n’ of the polyharmonic 
functions which constitute its real and imaginary parts. It is also clear that the number 
n’ is less than or equal to n, and that n’ is the highest power plus 1 of the product z Z 
existing in the terms of 2, . For later reference, we shall call these numbers n and n’, 
“order of analyticity’’ and “order of harmonicity” respectively of the function Q, . As 


examples we can cite the functions 


which are both of order of analyticity n, and have as order of harmonicity | and n 
respectively. The coefficients a, are arbitrary complex constants. 

3. Differential relations between two conjugate polyharmonic functions. Let 
P and Q be two conjugate polyharmonie functions of order m which we assume to be 


the real and imaginary parts of an analytic function of the same order m, ©,, ; that is, 
Siw P + iQ. The equation 0”Q,,/02 0 or (0/dx + 70/dy)"(P + 7Q) = 0 is equivalent 


to (1). If 6;” and 63” designate respectively the real and imaginary parts of the operator 
(d/dx + 7A/dy7)”, then we have 
a T 165" 0/0x + 10 Oy) re (3) 


The expansion of the right-hand side of (3) by the binomial formula furnishes for 6;” 


and 6,” the following expressions: 
0; : (— 1) (7") O /OXx OY 
—— \2y, 
[N m 2 if m is odd and A m/2 if m is even], 
, 
>» (— 1)’ | = rae y””** 
—_ mere Be OL oY 
eat 2v+ 1 
[N m 1)/2 if mis odd and N (m — 2)/2if mis even], where (%) is the 
number of combinations of m distinct objects taken k at a time and (>) is assumed to 
be equal to 1. Then writing (6;"" + 76 P + iQ) = 0, we obtain the partial differential 
equations 
gi P ig Q, i" P — = 5," Q, (4) 


which constitute for two arbitrary polyharmonic functions P and Q of the same order 
m, a necessary condition that they be conjugate or, in other words, that they be the 
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real and imaginary parts of an analytic function of order m. In this respect the relations 
(4) correspond to those of Cauchy-Riemann for two conjugate harmonic functions’. 

4. Real analytic functions of order n. We shall next establish the general form 
of the real analytic functions of order n which occur in certain applications. The problem 
is to determine functions g,(z) in such a way that the expression 


> 29.2) 


vel 


Wn 


be real. It is known a priori that the functions g,(z) are polynomials in z of degree n — 1 
with suitable coefficients. The latter can be calculated in the following way 


Wr a (aiz" | + ae? + ees tat 2+ aj) 
1 
+ 2" *(aje" | + aye”? + + aj) 
he 
1 2(a,,_ 2" + ay” + + ay-1) 
+ (ai2""" + qe2" + + a), 


or, more compactly, 


n n 
a,z 
1 


vel yp 


rv np 
~ 


, (5) 
where a; is real and a; a;. The expression (5) is the general form of a real analytic 
function of order n; it can be observed that this expression is a special polyharmonic 
function and that the order of analyticity and the order of harmonicity are necessarily 
equal for such a function. 

5. Definition of the transformation T,. Application to the analytic functions of 
ordern. Let p(z) and g(z) be two polynomials in z both of the same degree k and having 


no common tactor, 


kh 7 k 
sé a,2’, q(z) = p Ds b,2’, 


and replace in an analytic function of order n, 


P\z) = 


n—1 


> 2’f,(2) 


vO 


‘~ 
~ 


> and 2 by p/q and p/@ respectively. This substitution will be called a “transformation 
07 


7’,”’. The effect of this transformation on Q, is expressed by 


Theorem 1. If Q,(z, Z) = b ' 2” f,(z) és an analytic function having as order of 


analyticity and order of harmonicity the same number n, then 
p/qG) is an analytic function having as its two orders the same 


Proof. When in Q,(z, 2) the variables z and Z are replaced 
” ‘ “gee > «(2m) (2m) . °. - /5 
§ One has also m == $1 +6, , for VV" = (0/dx + 
«(m) +. (m) (2m) Ls 2m) 
(0) — 203 = 0] T O02 
form (6(°” + 63°") H = 0. 


ia/dy)™ (d/dx — i0/dy)™ 
; and the condition that H be a polyharmonic function of order m takes the 


the expression (p)"~* Q,(p/q, 
number k(n — 1) + 1. 
by p/q and j/@, respectively 


(af + 183”) 
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this function becomes 
n l 


2,.(p/q, b/@ = >, (b/O'f.(p/@ 


r=0 

n-1 

=1(""* LPT’ $.0/9; 
y=0 
hence we obtain ‘ 
n—-l 

g)"'2,(p/¢, p/D = be a’ f\p/q. (6) 
In (6) the product p’ g"~’*” is a polynomial in 2 of degree k(n — 1) and f,(p/q) is 


a function f*(z) of z, then the right-hand side of (6) is of the form 


1) 


| hs 2’F,(z) (7) 


) 


which represents an analytic function of order of analyticity k(n — 1) + 1. Now writing 


F(z) 2" = Gz , (7) becomes 


- 2*2°G(z). (8) 


From (8) it follows that the new analytic function also has the number k (n — 
as well as order of harmonicity. Hence the proof is complete. 
in general changes the two orders of 2, ; these 


Special case. The transformation 7 
1, so that the linear transformation 7’, , 


orders remain unchanged if and only if k 
(az + b)/(cz + d) where a, b, c, d are arbitrary complex constants, gives from 2, another 
analytic function having the same number 7 as both orders. 

6. Effects of the transformation T, on a polyharmonic function. Let us now consider 


a real polyharmonic function of order n of two variables and suppose that it constitutes 


the real or imaginary part, say real part, of an analytic function of the same order n. 


The effect of the transformation 7, on such a function is expressed by the following 


theorem 


Theorem 2. If H(z, 2) is a real polyharmonic function of order n, then (qq)""* H(p/q, 
p/q) is a polyharmonic function of order k(n — 1) + 1. 
Proof. Writing that H(z, 2) is the real part of an analytic function of order n, 2, , we 
have 
H(z, 2) 1/2, + 2); (9) 


then replacing in (9) the variables z and Z respectively by p/q and p/g, and multiplying 


both members by (qq)"~", we obtain 


(q@)""H(p/q, B/D = 34(G""2,(p/a, B/D) Q"" + (q""2,(G/4, p/@)7""}. ~~ (10) 


By Theorem 1, the term (g)""* 2,(p/q, #/@) in the right-hand side of (10) is an analytic 
function having as both orders the number k(n — 1) + 1. The multiplication of this 
function by g”* gives another analytic function of the same order. On the other hand, 
the term q”’ 2,(5/9, p/q) (@)"" is the imaginary conjugate of the preceding analytic 
function. It follows that the right-hand side of (10) represents a polyharmonic function 


of order n. This completes the proof. 
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19: 


qr 
~I 


The special case mentioned above for Theorem 1 holds also for Theorem 2; that is, 
the linear transformation (az + b)/(cz + d) conserves the order of polyharmonic func- 
tions. Theorem 1 may be considered, in the case of two variables, as a generalization of 
the similar theorems previously established by Lord Kelvin [5], Levi-Civita [4], P. 
Burgatti [2] and Nicolesco [8]. 

7. Nets of principal lines attached to analytic functions. Let p and —26 be the 
absolute value and argument of an analytic function of order n, 2, , and write 


n—l 
Q, = pe = D1 2'f,(2); 
v=0 
the values of p and @ are given by p = (2,,)' and 
0 = dare tgi ™—* = 44 log 2. (11) 
i > Q ow 


sen vin 


[f the function 2, contains real factors (constant or not), then the value of its argu- 
ment is evidently independent of those factors and is always equal to that of the complex 
factor which 2, contains. Suppose, for instance, that Q, is the product of an analytic 


function of order r, 2, , which contains no real factor, and a real analytic function of 


order s, w, ; it is clear, by Eq. (1), that r = n — (s — 1) and that Q, = a,. 2,-¢.-1). 
Since w, is real, we have arg 2, = arg 2,_;,—:). In the following, the order of harmonicity 
h of 2,_;.-:) , that is the order of harmonicity of the complex factor contained in Q, , will 


also be called ‘order of the argument of ©,”, and we shall use the notation —2é6, to 
designate an argument of order h. 

The application of the transformation 7, to the argument of an analytic function 
2,, gives a simple property which may be observed as a corollary of one of the preceding 


theorems 


If —26, (z, 2) is the argument of an analytic function of order h, Q, , then —20, 
p/q, p/q) is the argument —286, of an analytic function whose order of harmonicity is 
less than or equal tol = k(h — 1) + 1. The proof of this result is easily obtained when 


Theorem 2 is applied to both terms of the first expression of 6 in (11), or when Theorem 
1 is used for the second expression in (11). 


We now consider the family of curves (C,) in the plane (x, y) whose direction angle 
with the z-axis at each point is equal to 6, —20@ being the argument of an analytic func- 
tion 2, . The curves (C,) with their orthogonal trajectories constitute an orthogonal 
net of plane curves which we call ‘‘net of principal lines attached to Q,”’. As an example, 
when h 1, we have an isometric net of plane curves. We shall be concerned with 


certain properties of nets of principal lines which can easily be deduced from those of 
analytic functions of order n. 

Theorem 3. Let (C,) be a net of principal lines attached to an analytic function of order 
n, 2, ; it is possible to determine real analytic functions w,, of arbitrary order m such that 
the function w, . 2, admits (C,) as a net of principal lines. 

Proof. Let 2, = pe~?*** = >-"=) 2" f,(z) be the given analytic function and let 


m m 


Oo, = > z ary 


v=1 p=l 


be a real analytic function of order m; the argument of the analytic function of order 
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n+ m Seer w®, 2, Will be arg Qy4m—1 arg QQ, - 20, ; if (C,) is assumed 
to be a net of principal lines attached to an analytic function of order N = n + m - 


then 


Special case. When n 1, we have an isometric net (C,), and Qy is of the form 
Q, F(z) 2 >> ate” pla (12) 
‘ M | 
where F(z) is an ordinary analytic function whose argument is the harmonic angle — 20. 


We shall end this seetion with the extension of a property of the ordinary analytic 
functions to the analytic functions of order n. This will be the subject of the following 
problem. 

Probl m. Consider two analytic funclions of the same orde r Determine the conditions 
under which the arquments of these functions are two symmetric angle s with re spect to the 


x-axis. Let 


() ST sf ) «()* ST ar 
rs hat ~ J ? _— . —_~ J * 


be two analytic functions both of order n + 1, and let 


a ' 
5; og (2 Q log (Q%,,;/0%.:)] (13) 
~ al 
be their arguments. The condition that the angles (13) be symmetric with respect. to 
the x-axis is equivalent to Q . OF, Q . OF, a real quantity. Hence, we have to 


form the product Q,,, Q*,, which will be a real analytic function of order 2n + 1, w».,., 


(Sects. 1 and 4). Therefore, 


Ww) S* BP (2) f¥*(z > >. @i pict 14 
homed heed JH / » had = heed , 
V ) Mu l I 
Equating now the coefficients of 2”"(m 0, 1, +--+ , 2n) in the two members of (14), we 
have 
> Bae y ® a ae! m 0 
or iii (15) 
> Fe 2. a> 12 . . © ane mam > mM. 
When m varies from 0 to 2n, (15) furnishes a system of 2n | equations between 2n + 2 
unknown functions f, and f* (u,v = 0, 1, --- , n). If one takes, for instance, f,/fo and 
2n for 


f*/f* as unknowns of the problem, the number of the unknowns reduces to 
2n + 1 equations. Suppose that the unknowns /*//* are eliminated from the system 
fo will remain for n equations. On the other hand, Iq. (15) 


(15), the m unknowns f, 
fo or f*/f* will satisfy 


are obviously symmetric with respect to f, and f* ; then either f, 
the same conditions. We do not want to enter into the detail of this discussion. A solution 


of the problem is known a priori; it consists of the isometric case, i.e. when the angle 6 
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is a harmonic function. In order to indicate an example of a non-isometric case, we shall 


give the following functions 


nilon 1 

’ ’ =” k+1on i+1 

2... = Fe) > > ¢,, 27°" ' 
kh 1 | 1 


n+l ntl 


’ -n +l n-k+l1 
.. 1/F(z) > > tas” # : 
i 1 


(16) 


kewl 


where the coefficients C,,; are arbitrary complex constants. 

8. Applications to plane elasticity. Since the stress deviator of a two-dimensional 
elastic field is an analytic funetion of order 2, Q, = Zf,(z) + fo(z), the above develop- 
ments on the analytic functions of order n may be used, as far as the geometry of the 
elastic field is concerned, for the problems of two-dimensional elasticity. The net of the 
lines of principal stress constitutes a net of principal lines attached to the analytic 
function ® which represents the stress deviator of our elastic field; therefore, the theorem 
established for the nets attached to analytic functions of order n can be applied to the 


case ol plane elasticity. 


1. We consider first the special case of Theorem | for which & = 1. Thus we have: 
If Q, Sf t- f(z) ¢s the stress deviator of a two-dimensional elastic field in a domain 
D, then (62 + d). Q.) (az + b)/(cz 4+ d), (42 + b)/(é2% + d)] is the stress deviator of another 


two-dimensional clastic field in the domain D! into which D is carried by the transformation 
1 

B. Considering now the special case of transformation 7’, in the property concerning 
the invariance of the argument of an analytic function of order n, we can state: 

Let 0(2, 2) be the direction angle of a plane elastic stress field defined in a domain D; 


Lhe angle 


s the direction angle of another plane elastic stress field defined in the domain D' into which 
D is carried by the lransformation ‘ie . 

C’. Let us now take Theorem 3 in the special case of isometric nets. For plane elasticity 
this theorem takes the following form which states, as a special case of the general 
property, a well-known result. 

Let (C,) be an isometric net of plane curves defined by means of the argument of an 
analytic function Q, F(z); it is possible to determine real analytic functions of order 2, 
v» , such that the plane elastic stress field defined by the analytic function Q, = w, Q, accepts 
C’,) as a net of lines of principal stress. 

The corresponding analytic function Q, of order 2 is determined by Eq. (12); in fact, 


taking m 2. N = 2in this formula we obtain 
Q, = F(z)(a\zz + abz + aiz + a3), (17) 


where the coefficients have the same meaning as in Sect. 4. This question has previously 
heen discussed in different ways by P. I’. Nemenyi [6], U. Wegner [11], W. Prager [9], 
S. Siiray [10], R. Legendre [3] and P. F. Nemenyi and A. W. Saenz [7]. 

D. The general problem discussed above has been solved in the special case of plane 
elasticity by Nemenyi and Saenz [7]; these authors also gave an example of a non- 
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isometric case of this problem. A more general example may be deduced from the ex- 
pressions (16) if one takes n = 1: ; 

Q, = F(e)(aysFZ + ayo% + ioe + Ao), 

QF = 1/F(z2)(G,:22 + dig z + doo), 


where the coefficients are arbitrary complex constants. 
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THE REFLECTION COEFFICIENT OF A SURFACE OF RAYLEIGH 
DISTRIBUTED IMPEDANCE* 


BY 
H. 8S. HEAPS 
Nova Scotia Technical College, Halifax, Canada 
Abstract. Formulas are obtained for the probability distribution P(k) of the ampli- 
tude k of (z — c,)* (z + ¢,)~* in which ¢, and ¢, are real constants. z is a complex number 
of the forme + z, where c is a complex constant and z, is a complex variable of Gaussian 
distribution of amplitude and uniform distribution of phase (see Fig. 1). 





[ 





L 
=C O S 
Fic. 1. Geometrical configuration of the vectors z + ¢; and z — ce. The parallelogram containing dx; 
and dx, represents the area swept by the end point of x; as x; ranges between x, and x; + dx; and Xe 
ranges between x. and x2 + dx, . The end point of z ranges over an equal area. 
The above relation between k and z arises when c, = cs in connection with the re- 


flection of a plane wave of sound from a plane surface and relates the power reflection 
coefficient k to the acoustic impedance ratio z. 

The present paper expresses P(k) in terms of the modified Bessel functions J, and 
1, ifc is real, and in terms of the incomplete gamma function if c is complex. Asymptotic 
formulas are included to allow calculation of P(k) beyond the range of tabulation of 
the modified Bessel and incomplete gamma functions. 

1. Introduction. The concept of surface impedance is widely used to describe the 
behavior of a surface with regard to the reflection of acoustic or electromagnetic radiation. 
The normal acoustic impedance is defined to be the ratio of the acoustic pressure at the 
surface to the particle velocity normal to the surface [1-3]. Division by pc, where p is 
the density of the propagating medium and c is the sound velocity, results in the acoustic 
impedance ratio. Similarly the electromagnetic impedance is defined to be the ratio of 
the tangential electric field at the surface to the tangential magnetic field in the direction 
at right angles [4, 5, p. 80]. The impedance ratio is then obtained upon division by 
(u/e)* where uw and ¢ are respectively the permeability and dielectric constant of the 


propagating medium. 





*Received August 21, 1956; revised manuscript received October 22, 1956. 
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In the acoustical case the velocity potential y of the field at the reflecting surface 


depends upon the impedance ratio ¢ through the relation 


V/(dW/dz ¢/(tkp), 1) 
where i, is the wave number of the incident radiation and z is the direction of the out- 


ward normal to the surface. In the reflection of a plane electromagnetic wave the com- 
ponent of electric he d parallel to the surface and at right angles to the direction of 
incidence satisfies | q. (1) at the surface. In either instance the impedance ratio is related 
to the power reflection coefficient / of the surface with regard to a plane wave incident 
by the relation 
| C COS a | + ¢ cosa ys 
where a is the angle of incidence measured from the normal to the reflecting surface’ 
Graphs of the dependence of / upon ¢ for various complex values of ¢ are given by 
Morse [2, p. 367] 
Suppose that the surface impedance is not uniform but that its variation is in- 
significant over a distance equal to the wavelength of the incident radiation. Then the 
reflected radiation received at a distant point is predictable by ray theory, and it 


consists of a plan e reflected with a power reflection coefficient given by (2) in 


which ¢ is the impedance ratio of that portion of the surface that may be viewed at 
the angle of reflection a from the receiver [6]. As the receiver moves in space the variation 
of received power 1S & measure olf the variation of | | C COS Q@ | + ¢ COS @ over 
the reflecting surfa 

lor manv non-uniform reflecting surtaces the impedance of any portion cannot be 
measured directly since the placing of apparatus near the surface affects the local sound 
or electric field. In general ¢ has both real and imaginary components while (1 C COS a 
(1 + ¢ cos a) may often be measured in terms of amplitude only and not with regard to 


phase. It is consequently of some practical interest to consider the manner in which the 


probability distribution of amplitude and phase of ¢ affects the distribution of | 1 COS 
al*|1+ cosa 

The present papel derives the probability density of the amplitude k of | z ( 
ze in which ¢c, and c, are real constants. z is a complex number of the forme + z 


where c is a complex constant and z, is a complex variable of Gaussian distribution of 
amplitude ind uniform distribution of phase. The geometrical configuration of the 
vectors which represent these complex numbers is as indicated in lig. 1. 

2. General formula for the probability density of k. If z has Gaussian distribution 
of amplitude z,; and uniform distribution of phase 6 the joint probability density of z, 
and @ is 


|/mo )2,; exp | z1/a ) (3 


in terms of the standard deviation o of the amplitude. The factor 2, appears in (6 since 
2,dz,d0 represents the area swept by the end point of z,; as 2, and @ vary between z, and 
z, + dz, and @ and @ + dé respectively. 

Denote z + ¢, by x, and z c, by X, (see Fig. 1). Possible values of x, and x, are 


those that permit a real solution z of the equation 
| | 


te 2 — 2¢.x7. cos 6 
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where c, denotes c, + ¢ and @ is the angle between x, and — c, . Equation (4) has two 
solutions in z if 
la, —a|<lea| <2,+ 2, (5) 
and no solution otherwise. As x, varies between z, and x, + dx, and x, varies between 
and x, + dx, the end point of z, and hence of z, , ranges over an area of 
| 27-3 
Ly, Xol[4ajp x2 ¢ ZZ, =~ Le) - dx, dx “ 
[aie yy de, a 
[42323 — (c; — xi — x3) ] , d(x}) d(x). 
Thus the joint probability density of 2{ and x; is 
l wo’) [4ajx> _ (ce; —2- x) | , 4 — Zi a). (7) 


Because the element of area (6) includes pairs of solutions of (4) it may be supposed in 
7) that the phase of x, and x, lies between zero and z. 


Now put a ka? and d(#2) = x? dk. The probability density of k is P(i) where 
Pk |/ro | [-4/ "ri — (c; — x; _ kx})"] by? exp (— 2 /a°) d(x}), (8) 
with integration over the range of 27 from ec, (1 + k*)-? to ce, (1 — k*)7?. 
The substitution a e,(1—k) "fl +k - 2h: cos ¢] reduces (8) to the form 
Pil [ce; | 1 k |~*/xo"] (1 + k — 2k cos] exp (— z{/o”) do. (9) 
Before proceeding with the evaluation of (9) it is necessary to express z, In terms 
of . Ife = c exp (78) and z = z exp (zy) then 
c +2° — 2z cos (y — B). (10) 


Let c./c, be denoted by a. Then 


2 2 2 ¢ . 1 
COs ¥ (2; _ Ge 9 )(2¢,2) (11) 
fa+k-— ki + a) cos ¢]e,(1 — k) ie 
and, since y lies between zero and 7, 
siny = k*(1 + ae, | 1 —k|7'2 sing. (12) 


Substitution of the above expressions for cos y and sin y into (10) enables 2} to be 


written in the form 


ce + [4ak + (1 + b)\(a’? + A]ei(1 — &) 
— Qe,(a + k)\(1 — k)™" cos B (13) 


— 2k*(1 + ae,(1 — k)?R cos (6 — @), 
where ? and w are given by 
R e(1 — k)? + ci(a + hk)? — 211 — Aa + Adee, cos B (14) 
and 


cosw = [(k + ae, — (1 — Ke cos BIR", (15) 
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sinw = |1—k|cR'sinB. (16) 
Thus (9) becomes 
P(k) = [c3 | 1 — k |~*/o*] exp (— c’/o”) 
x exp {— [4ak + (1 + D(a? + KH] — WK ’et/o"} (17) 
< exp [2cc,(a + k)(1 — k)~’o™* cos 8] 


Pg | ai+k- 2k* cos ¢) exp [b cos (6 — w)] dg, 
where b = 2k? (1 + a)e, R(1 — k)7c”. 
Reduction of the formula when c is real. In the case that c is real and positive then 
8 = 0 and w = 0 or a. The case of c negative follows by interchange of c, and c, but 
does not arise when c denotes an acoustic impedance. The integral in (17) may now be 
expressed in terms of the modified Bessel functions I,(b) and J,(b) as #[(1 + k),(b) — 
2}:*T,(b)] [7, p. 181]. The probability density P(k) is thus 
P(k) = [ce3 | 1 — k |~°o *] exp (— e'/o’) 
x exp {— [4ak + (1+ Aa + 4] — &) “ci /a°} (18) 
xX ; exp [2cc,(a + k)\(1 — ko male + k)I,(b) — 2k*7,(b)], 
where b = 2k? (1 + a)e, [(a + k)e, — (1 — k)c] (1 — k)~*o™. This latter form of b 
includes the effect of the cos w factor which may be + 1. 


For values of b between 0 and 10 the functions J,(b) and J,(b) are tabulated in 


[8, p. 226]. For larger values of 6 it is convenient to make use of an asymptotic form 
for P(k). Such is necessary, for example, in order to evaluate P(1). For large values 


of b the asymptotic formula for J,(b) [7 p. 203] may be applied in the form 


I.(b) ~ (2rb)~*(exp b)[1 — (4r° — 1)/8b + (4r° — I(4r° — 9)/128b"), (19) 


| through the relation 7,(b) = (—1)’J,(—b) when b is negative. 


when b is positive and 


The resulting asymptotic expansion for (1 + k)I,(b) — 2k*J,(b) is 


1 + k)I,(b) — 2k*I,(b 2rb |*(1 — k*)*(exp | b | )F(k, 5), 20) 
whe 
F(k. b) ~ 1 — 3/8|b| — 15/1280? + (1 + Bi — k*)77/2 | b (01 
+ 3(1 + A — k*)-?/160". 
In the above equations k? denotes the positive square root of k when b is positive and 


the negative square root when b is negative. 
Substitution of (20) into (17) permits P(k) to be written as 


P(k ca? | 2eb |-4(1 — k*)? | 1 — k | °F (Kk, b)[exp (— c°/o°)] (99 
(22) 
x fexp [— (a — k°(1 + k*)*ci/o")} {exp [2(@ — k*)(1 + k’)e,e/o°)}. 
The equation (22) applies in particular when k = 1 and reduces to 
P(1 1/8)r73(c2/e,0)(1 + a) [exp (— c?/o”)]{exp [— 34(1 — a)’e:/o°}} (23) 
“0 


x fexp [(a — lec/o*}}(1 + 20°/c;). 
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For sufficiently large values of c it is permissible to put b = —2k* (1 + a) (1 — ky 
c,c/o and hence to write (22) in the form 


P(k) = 39 bio" [1 + @kecl* | 1 — kl] 1 +h | ORG, +) OG (24) 


in which the final brackets are as in (22). It is to be noted that if k > 1 then b is positive 
and P(k) contains the factor | 1 — k* |~*. However, if k < 1 then b is negative and P(k) 
contains the factor | 1 — k* |~*”?. Thus for large c there is greater probability of finding 
a value of k slightly less than unity than of finding a value slightly in excess of unity. 
This is in agreement with the fact that ce — ec. <<ce+¢,. 

For sufficiently large values of c/o and c/c, the value of b is large compared to unity 
and the order of magnitude of P(k) in (22) is governed by the product of the exponential 
terms. The dominant term is exp (—c’/o’) and the product is negligible save for values 
of k near to the solution of the equation 


2a — kK) + k)ec/o? — (a — K)7(1 + 2)? /0? = c?/o”. (25) 
The solution of (25) is 
k = (¢ — ae)" + ¢)” 


which is also the discrete value of k obtained when o = 0. 

4. Reduction of the formula when c is complex. For complex values of c the integrals 
appearing in (17) cannot be evaluated in a closed form, and for purposes of numerical 
calculation it is desirable to express them as rapidly convergent series. The required 
integrals may be written as follows 


a® 


| exp [b cos (@ — w)] do = [ exp (b cos ¢) do 
J0 0 


(26) 
+ [ exp (b cos ¢) do 
and 
[ cos @ exp [b cos (¢ — w)] dd = cosw | : cos ¢ exp (b cos ¢) dd 
+ cosw | 3 cos @ exp (b cos ¢) dd (27) 
0 


— (2/b) sin w sinh (b cos w). 
Following the method of Hadamard [7, p. 204] put ¢ = 2b sin?(3¢) and let b, denote 
2b sin*(4w). Then 
P exp 6 (2n)! mor ; 
| exp (b cos ¢) d@ = ae » > is i "i exp (— 0) dt 


n=0 


(28) 


exp b (2n)! 
= Opi Le Higinpe Ven + 4, bi), 


where y(n + 34, b,) denotes the incomplete gamma function. Similarly 


ps _ exp b =“ (2n)! 2 = ; 
| COS ? exp (b COs >) do _ ap) (n!)? apps = eae ; y(n +3 2? b,). (29) 
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Thus the integrals (26) and (27) may be expressed in the form 


[ exp [b cos (@ — w)] do 
0 


(30) 
exp b < (2n)! ‘ : 
— —s-sazn [y(n + 4, 61) + y(n + 3, 5.)] 
(am)? 2 G2 b | , sid 
and 
[ cos ¢ exp [b cos (6 — w)| do 
v0 
b at a. (31) 
; exp 0D (Zn)! Zn , , 
—( COS w) , 6 Ey - y(n + 4, b,)--4 (n 5, b.) 
(2b)? 2. Gl)2"b" On — {i 21 Oi) sr ym F 25 92) 
— (2/b) sinwsinh (b cosw), 
where b, = 2b cos*(3w). 
The function y(n + 4, b,) may be found from Ref. [9] over a wide range of values 
of n and b, . Although b” appears in the denominators of (30) and (31) the series are 
: : . . n+} 
convergent for small values of b since y(n + s, b,) tends to zero as b; */(n : For 


large values of b, or b, the value of y may be found from the asymptotic formula [10, 
p. 13 


oO] 


,; 
y(n + 4, b,) ~ Tin + 3) — OY? exp (— J,). 


When b, and b, are both large compared to unity the first terms in the series (30) 


and (31) may be combined with (32) to give 


| exp [b COS (@ — Ww | do ~ 2n(2rb) exp b (33) 
| cos @ exp [b cos (6 — w)] de ~ (cos w)2n(2rb)* exp b (34) 
and hence from (17 
P(k ~ [ l _ } Oo exp . a C go) 
xX exp {-— [4ak + (1+ A(a +h] — k) ei /o'} i 
oo0) 


X exp [2cc,(a + k)(1 — hk) ao ~ cos B] 
x 2(2rb)*(exp b |] -+- k- 2k cos w]. 


5. Concluding remarks. For given values of the parameters c,, c,, ¢, 7, and 8, as 
k ranges from zero to unity, the range of b is from zero to infinity. Thus, in general, in 
order to plot the graph of P(k) it is convenient to use Eqs. (30) and (31) when k is small 
and Eq. (35) when k is close to unity. The case of k > 1 has no physical meaning when 
R in (2) denotes a reflection coefficient. It may be remarked, however, that the quantity 
|z — c, |” | z+, |~* is also related to the problem of the detection of a constant signal 
C, upon a noise background z. It represents the ratio of the power of the signal minus 
noise to the power of a multiple of the signal plus the same noise. In this instance there 


is no physical restriction to values of k less than unity. 
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— NOTES — 


ONE-DIMENSIONAL DIFFUSION WITH THE DIFFUSION COEFFICIENT A 
LINEAR FUNCTION OF CONCENTRATION: REDUCTION TO AN 
EQUATION OF THE FIRST ORDER* 


By D. H. PARSONS (University of Reading, England) 


Abstract. The problem considered by Stokes', of one-dimensional diffusion from 
an initially sharp boundary between two semi-infinite columns of liquid, the diffusion 
constant being a linear function of concentration, is discussed. It is shown how the differ- 
ential equation may be reduced to an equation of the first order. Some properties of the 


solution are investigated, and the method of obtaining numerical solutions is considered. 





The differe ntial equation to be satisfied IS 


1 ac le 
2£(p J+y—=0, (1) 
dy : ay : dy 


with bou darv conditions 
c— ¢, when y > — 
( > when 7 > CO 


) 


where c is the concentration, c, and ce, the initial concentrations of the two columns, D 


a 1Y i . ° . ° 
the diffusion coefficient, and y at™*, where x is the distance from the boundary in the 
direction of diffusion, and ¢ is the time (see Stokes, loc. cit.). Let D* be the mean concen- 
tration, and let b = ratio of diffusion coefficient at concentration c, to that at c, (6 ¥ 1). 
Then, since we assum¢ that D is a linear function of ¢, we may write 
] bh jile - ¢, -C j 
D D*¥i 1-4 “ie acai ’ (2 
, 1 + 0 | (Cy - Co) 
and, r ti sul tution in (1), al d also writing y oe. z, the equation to be 
satis be 
] h L(p Lc.) — el jd | — pb | /de\? de 
=> | — ——= = | +- 2z 0, (3) 
ee oa aU Co) ) | lz 1+ b g\C; — C \dz dz 
with tl liti c, when ( Cc. when z 
Let 
l 
es TT (4) 
(1 + 0 
Cy Co a c 
C= + 7 eee 4 (5) 
— }) 


*Received March 12, 1956. 


IR. H. Stokes, Transactions of the Faraday Society, No. 358, Vol. 48, Part 10, October, 1952. 
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Then, substituting these expressions in (3), the equation to be satisfied reduces to 


d’w/duw? + u dw/du — 2(1 + w)(dw/du)? + u(3 + u’)(dw/du)*® = 0, (6) 
or finally, writing p = — dw/du, 
dp/du + up + 2(1 + u’)p* + u(3 + u*)p*® = 0, (7) 
an equation of the first order. It is clear from (5) that c = c, when w = log (b7') and 
c = c, when w = 0. Consequently, from (4) and (5) the new boundary conditions are 
w — log (b+) when u > — @, 


w—+ QO when u—- ©, 
We may remark here for reference later on, that we have 


—1 de —2(1 +b)” isl dw/du 


C, —¢,dz 1—b 1 — udw/du 
2(1+ 0b)” _, p 
‘ ot. Na = oe 
1—b 1 + up 
Now the equation (7) is an ordinary differential equation of the first order and the 
normal form. Its solutions depend upon one arbitrary constant, which may conveniently 
be taken to be the value of p when u = 0. It follows from the form of (7) that p cannot 


(8) 


be zero for any real u. For if p = 0 when u = wp , say, then from the form of (7) this 
integral is certainly analytic at the point uv = wp , and is therefore simply the integral 
D 0. p would therefore be identically zero for all u, which makes w, and hence ec, 


constant throughout, a trivial solution which we reject. We may thus suppose p positive 
throughout: for a solution for which p is negative throughout is obtained by writing 
p’ for p and —w’ for u, which does not change the form of (7). Suppose, therefore, 
that p = \ > 0 when u = O. Then the following results, the proofs of which are given 
in the appendix, may be established: 

(i) Every point of the real axis of u is a regular point of the integral. 

(ii) When u increases from zero to ©, p decreases strictly, and p> 0 asu —> @, 
while up increases strictly from zero to a single maximum at a point in the range 
0 < u < 1, this maximum value being < 1/3, and then decreases strictly, and — 0 
as u 
(iii) When u decreases from zero to — ~, p first increases strictly to a maximum at 
a point in the range —1 < u < 0, then decreases strictly, and — 0 when u— — o. 
But, according to the value of \, either up decreases strictly from zero to a single minimum 
value, this value being > —1, at some point for which u < —1, then increases strictly, 
and — 0 as u — — ©; or else up decreases strictly whenever u < 0, 
and — —1 as u— — ©. For all small enough \, the former circumstance obtains. 

(iv) {3 pduand {% up du are convergent, uniformly with respect to \, in any positive 
range of X. 

(v) [°..pduand f°. up du are convergent uniformly with respect to \, if A lies in a 
range such that wp — 0 when u—> —o@., 

Suppose then that we select a positive value of X, and calculate numerically the 
solution of (7) for which p = \ when u = 0. It is essential that the calculations be carried 
on far enough, in the negative direction, to make sure that up passes its (negative) 








300 NOTES [Vol. XV, No. 3 
minimum value, and thereafter increases towards zero; for otherwise we cannot be sure 
that | up | tends to zero and not to unity as u — — ©. If this condition be not satisfied, 
we simply start again with a smaller value of X: for it is proved in the appendix that 
for all \ up to a positive limit, wp — 0 as u > — o. With this understanding, we have 
dw/du = —p, and hence w = A — fj pdu. But since we require that w— 0 when u > @, 
and since the integral is convergent, we see at once that this gives 


a2 


w= | p du, (9) 


vu 


and thus w is known to any desired degree of accuracy, from the calculated solution of 
(7). This latter being, by hypothesis, such that [*.. p du is convergent when u— — o, 
the second boundary condition shows that 


b = exp (-2 | p au), (b < lsincep > 0). (10) 


The constant } is therefore calculated from (10), and clearly depends only on \. From 
(4) and (5) we then have the required solution parametrically, in terms of u, in the form 


z= eae u exp (~| p au), (11) 
_ 1 — exp (-2 p au) 
Snack. J u (12) 





i-e, 1— b . 


b being given by (10): and from the form of these expressions and from (10), we see that 


c—c, when u— —© andz— —”,c—c, whenu— ~, z— o, so that the original 


boundary conditions are satisfied, and we have the required solution of (1). 
We also have, from (4), (8), (9), and (10) 
’ VS/) ’ ’ 


bo 
bo 
3 


3 up 


é —E= 2s 








2 





Z 


z 


—] de ( x) 
aa ——i——— exp | — ) 
dz 1—b 1 + up 1—b1l+wu nied 2 | p du (13) 


C; — C2 


/dz—>0 


and since 1 + up > 0 throughout’, and wp — 0 when u — + ©, we see that zdce 
when z— + ©. 

The family of curves of c against z, or the gradient curves — (c, — c.)'‘de/d 
z, can then be plotted for various values of b, by means of the expressions (11), (12), 
(13), which give these expressions in terms of the parameter u. The method is firstly 
to calculate the solutions of (7) for an arbitrary positive \; then to calculate b from (10); 


Zz against 


and finally z, c and de/dz in terms of u, from (11), (12), (13). 

The advantages of the foregoing method of procedure are two-fold. First of all, 
there is the advantage of calculating a solution for a first-order differential equation in- 
stead of a second-order one. Secondly, the single arbitrary quantity is the value of p at the 
point u = 0, orz = y = x = O. Thus the calculations are begun at a point corresponding 
to the boundary between the two columns of liquid, and proceed symmetrically in both 
directions. Thirdly, the convergence of the integrals required to give b, and w in terms 
of u, has been shown to be uniform with respect to \ when \ does not exceed a fixed 





*See Appendix. 
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limit, H say. And denoting the solution of (7) which takes the value \ when u = 0 by 
p(u, A), it is proved in the appendix that p(u, 4) < p(u, uw) when X < 7. Then if 0 < u, 


u, and uz; < us < 0, we have 
0< | plu, h) du < | plu, H) du, 
and 


0< | ‘plu, h) du < [ ‘pu, H) du. 


It follows that, when calculating the infinite integrals which occur in (10)-(13), if we 
carry the calculations far enough in each direction to attain the required degree of 
accuracy for \ = H, then carrying the calculations the same distance will certainly 
give the required accuracy for any smaller value of X. 

I am very much indebted to Professor Stokes who drew my attention to this inter- 


esting problem. 
Appendix 


We require to prove the properties (i)-(v), listed above, of that solution of (7) which 
takes the value \ > 0 when u = 0. 


(i) The equation (7) has the fixed singularities u = 0, u = + 34. The integral 
such that p = \ > O when u = O is analytic when u = 0. Consequently, it follows from 
a result of Goursat*®, that a point « = uy on the real axis can only be a singularity if 
p— © when u — u . We have shown that p > 0 throughout. But from (7) it follows 
that for all large enough positive p, dp/du < 0 when u, > 0 and | u — wu | is small, and 
dp/du > 0 when uy < Oand|u — uw, | is small. This contradicts either of the hypotheses 


that uw. > 0,p— + © when u — wu, from below, or that u. < 0,p— + © when u — up 
from above. It follows that every point of the real axis is a regular point of the integral. 
(ii) From (7), when u > 0 and p > O, dp/du < 0, and p decreases strictly. But 
p > 0 throughout, therefore p — a > 0 whenu— + ~.Ifa> 0, clearly dp/du —- —o, 
contradicting the hypothesis that p — a. Hence p > 0 asu— @~, 
Next, let » = up. Then, using (7), we have, after simplifying, 


u dv/du = v1 + {1 — wv) — (8+ ur}. (14) 


When u > 0, v > 0, and therefore, from (14), when u > 1, dv/du < 0. It follows that 
when u— ©,v—- 8 > 0. But if B > 0, dv/du ~ —B(1 + Bu — ©, contradicting 


v — B. Therefore, »v > 0 as u — ~. To show that v has only one maximum, occurring 
when 0 < u < 1, we first observe that when u = 0, dv/du = pp = \ > 0. When dv/du = 0 
and u > 0, we have, from (14), d’v/du? = —2v(1 + v)? < 0. Hence there is only one 
maximum value. At this value, from (14), we see that 

Ynez = (1 — u’)/(3 + u*?) = 1/3 — 4u7/3(1 + u’) < 1/8. 


Thus proposition (ii) is established. 





(iii) When u < 0, we may write the equation (7) in the form 
pe —_— (< 2 2) 2 i 5 
p/du 13 +0) [{u(3 + u)p + (1 +w)}* + 1)], (15) 


3. Goursat, ‘‘Mathematical Analysis”, Vol. II, Part II, sect.‘67, pp. 182-184 (Ginn and Co., 1945). 
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and from this we see that when wu < 1, dp/du > 0. Since p > 0, it follows that p — 
y > 0 when u — ©. Ify > 0,dp/du— © asu— ©, contradicting p — y. There- 
fore p— O0as— — ©, Again, when u = 0, dp/du 2* < 0. There is thus a maximum 
value between 1 Oand u 1. To show that there is only one maximum, we observe 


from (15) above that when dp/du 0 and i<u< @ 
d’p/du Dil + dup + 311 + w)p*} 

) ; : 

I {30 + w)p + Qu}? +3 —w] <0: 


and the first part ol the proposition is proved, 


To establish the second part, we again write v up, so that v again satisfies (14) 
above, and if dv/du Qand u < 0, we have d*v/du? 2v(1 + v)*. But when u < 0, 
v < 0, and it follows from the form of (14) that v > | forany u < O: if 1 1 when 
u Uy < 0, the integral is certainly analytic at wv , and therefore v l, p 1/u, 
a solution which we may reject, since by hypothesis p A when u 0. Consequently, 
when u (0. | 0, and therefore if dv/du 0, d°v/du° > 0. Hence there can be 
at most one negative minimum value of v. But from (14) we see that when —1 < u < 0, 
dv/du > 0. Hence the minimum, if it exists, must oceur for u < b 

Since dv/du > O when u 0, and since —1 < v < Owhen u < 0, it follows that when 
u— ©, v either decreases strictly to a limit 6 where | < 6 < 0, or else decreases 
to a minimum and then inereases strictly to a limit e, where 1 < e€ < 0. In either 
case, suppose v tends to a limit 7, —1 < » < 0. Then dv/du ~ n(l + )° u 
asu— ©, contradicting the hypothe is that 2 >». Hence 7 0 or else 7 1. In 
other words, as u decreases from 0 to ©, either v decreases to a minimum value between 
0 and |, for some value of wu less than |, and then increases strictly to the limit 


zero, or else v decreases strictly to the limit 
We have yet to show that, for small enough A, the former event takes place, Since 


dv/du can have at the most one zero when u < 0, and since dv/du > 0 when u 0, it 
follows that if dv/du 0 for any negative u, then certainly v has a minimum, and 
increases to zero as u —> ©. But, selecting any fixed negative number N, where 
N « |, and observing that p 0 is a solution of (7), it is a standard result* that a 
positive number / exists (depending upon N), such that whenever 0 < \ < JL, the 
solution of (7) for which p = XA when u 0 is such that when 
uM N, 0 < p(N) - a a 
VN'+3 
and therefore 0 > (3 + N*)v > 1 N~. rom (14) we see, therefore, that whenever 
\ < L, dv/du < 0 when u N: and the result follows. 
(iv) Denoting the integral of (7) which takes on the value \ when u = 0 by p(u, A), 


we first observe that if \ < yw, then for all u, p(u, A) < p(u, w). For by the usual existence 
theorem, there cannot be two integrals taking a given value for any given value of u, 
so that {p(u, u) — p(u, A)} # 0. But it is impossible for {p(u, uw) — p(u, \)} to be negative 
for any u: for being continuous, and positive when u = 0, it would then be zero for some 


other wu. 


‘E. Goursat, ‘Cours d’Analyse Mathematique”’, tome ITI, sect. 461, pp. 16-18. (Gauthier-Villars, 
1942). 
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Now, integrating the equation (7) between the limits 0 and u, we have 


r p | jup(l + up)” {- (2p° } Sup )} du. (16) 
Suppose that 0 < 4 < M. Then we have proved that, for any \, up — 0 asu— @, 
and that p(u, A) < p(u, M). Hence p = o(1/u), p’ = o(1/u’), up* = o(1/u*): and there- 
fore [* (2p? + 3up")du is uniformly convergent as u— ©, But since p — 0 and p(u, d) 


p(u, M) as u eo, the left hand side of (16) is uniformly convergent to \. Therefore, 
since for all large enough u, (1 + up) > 4, it follows at once that {> wp du is uniformly 
convergent when u — ©: and the same is therefore true of {% p du. 

(v) We may write (16) in the form 


DN p | fup(l 4 up) -| (2p° {- 3up’)} du. 
We have proved that for all small enough A, say A < H, up > 0 as u — —, In this 
event, it follows, by precisely the reasoning of (iv) above, that f°.. updu and f-.. pdu 
are uniformly convergent, 


AN ELEMENTARY SOLUTION OF TWO STRESS CONCENTRATION PROBLEMS 
IN THE NEIGHBOURHOOD OF A HOLE* 


By W. T. KOITER (Technische Hogeschool, Delft) 


1. Introduction. ‘The problem of the stress concentration around a central circular 
tension (or in bending in its plane) has been investigated extensively. 


hole in a LTip 

The theoret olution of this problem has been given by Howland [1]. However, 
Hlowland’s numerical calculations were restricted to hole diameters which do not exceed 
half the strip width, and his method requires a rapidly increasing numerical effort for 
larger hole Experimental investigations have yielded valuable data on the stress 
concentration factor for large holes. On the basis of these experimental results Heywood 


that the stress concentration factor, referred to the nominal stress in 


conjecture LnAt 


the net section, tends to the value 2 if the hole diameter approaches to the strip width. 


This conjecture confirmed by the elementary analysis in See. 2. 

Recently the more complicated problem of a central spherical hole in a cylindrical 
bar in tension has been investigated by Ling [3]. His numerical evaluation has also been 
restricted to a hole radius equal to 1/4 and 1/2 times the bar radius. Here again the 
numerical effort, required in evaluation of the theory, increases rapidly with increasing 
hole radius. According to Ling “the value 1 for the stress concentration factor K (referred 
to the nominal stress in the net section) in the limiting case \ = 1 (where J is the ratio 
to bar radius) can be visualized readily from physical considerations of 


the cylinder” [3, p. 391], and the graph of his results [3, Fig. 2] has been completed for 


the entire range 0 - 
incorrect, as will be shown in See. 3. The correct limiting value of the stress concentration 


of hole radiu 
\ < 1. However, Ling’s argument, which has been cited above, is 


*Received May 18, 1956. 
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factor is 


nek = RSS. TL (1) 
as o — 4y 
where v is Poisson’s ratio. Ling’s graph for the stress concentration factor K as a function 
of X should be corrected accordingly. 

Our elementary solutions for the stress concentration problems for the strip with a 
large circular hole and the cylindrical bar with a large spherical hole are based on beam 
theory and shell theory respectively. It might be objected that both these theories are 
approximate in character. On the other hand, it may be contended that the accuracy 
of beam and shell theories increases indefinitely when the “slenderness” of beam and 
shell increase indefinitely, and the limiting values for the stress concentration factor 
K for \ — 1 should therefore be correct. 

Nevertheless, a more rigorous investigation is considered worthwhile, the more so 
because it may then be possible to investigate the manner in which the limiting value is 
approached. A suitable method is mentioned in Sec. 4 and some preliminary results for 
the strip problem, confirming the result of beam theory, are given. 

2. Strip with a large circular hole. ‘The strip with a large circular hole (6 « R) is 
pictured in Fig. 1, and the upper part to the right of the minimum section is given on a 
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larger scale in Fig. 2. It is obvious from symmetry that no shear force is transmitted 
by the minimum cross-section AB. On the other hand, in order to maintain this cross- 
section in a vertical position, which is again obviously necessary on account of symmetry, 
a positive bending moment M, must occur in this cross-section. 
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An element of the “beam” in Fig. 2 has been drawn in Fig. 3. Obviously the shear 
force is zero and the normal force has the constant value N. 
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dx 























Equilibrium of moments requires 


w &. (2) 


In the neighbourhood of the minimum cross-section z = 0 the height ¢ of the beam is 


approximated by 


2 
x 
t= 6+ 55 (3) 
Hence, Eq. (2) may be written in the form 
dM 1 x 
dz 2 a 3 (4) 


Because dt/dz is small in the neighbourhood of x = 0, we may apply the beam formula 


_ pp fe 1 ppp SY 
M = BITS = EMSS, (5) 


where w is the deflection (positive upwards), h is the strip thickness and # is Young’s 
modulus. Equation (4) now reads 


+ mm £(¢ e) =» ~eWz (6) 


2 6 dz\ dx” es 
The general solution of (6) contains three integration constants. Two boundary 
conditions are supplied by the requirement that w and its derivatives tend to zero for 
x — © because the deformation is negligible for large values of 7/5. The third boundary 
condition is that the slope dw/dx at the minimum cross-section x = 0 is zero. Trying 
a power series solution in ¢-* it is immediately seen that 


2NR (7) 


a. 
Eh(6 + 2°/2R) 
satisfies both the differential equation and the boundary conditions. 
The bending stress at the inner edge of the beam is given by 


aie Et ew. (8) 
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Putting 
H 


(R65) 


(Y) 


frr 


we obtain for the sum of bending stress and tensile stress (due to the tensile force NV) 


at the inner edge and outer edge respectively 
y 9 -2 4 9): 
N —é¢ \ 2 


oo -2 


- = 5 Bic — - s° (10) 
hé (1 + &°/2) hé (1 + &°/2)° 


The maximum stress occurs at the inner edge at x 0; its value is 2N/hé, and the stress 


concentration factor K 2: 

3. Cylindrical bar with a large spherical hole. Tl igure | also depicts a meridional 
section of a cylindrical bar with a large spherical hole (6 < R), whereas Fig. 2 may be 
regarded again as an enlarged picture of the upper part to the right of the minimum 
section. This part is now considered as a cylindrical shell of variable thickness ¢, given 
by (3). An element of this shell of circumferential angle dg is drawn in Fig. 4, where the 
bending moment m and the tensile and shear forces n and q are defined per unit circum- 
ferential angle. Obviously the tensile force n is constant. Equilibrium of moments 














md 
nd) 
| 
qd! 
ia. 4 
requires 
dm l dt 
- q=— 4 (11) 
dx 1 9" dx’ ‘ 


whereas the equilibrium of radial forces is expressed by 


dq 
lo. (12) 
dx ~ 
where og, is the tangential stress in the middle surface. From Hooke’s law we have 
w l ( A ‘ 
€ -=lo, — v=}, (13) 
a ae i Rt 
and we obtain 
, W 7 . 
a. iD rv ° (14) 
7 R Rt 
Finally the bending moment m is given by the shell formula [4] 
1 KHtR dw : 
a. . (15) 


” 121 —y dx 
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[liminating m, q and o, from (11), (12), (14) and (15) we obtain the differential equation’ 


_E & , dw (3 ) n 
12(1 — v”) dx’ ( ) + R?’ _— 2 ve R° (16) 


The boundary conditions are here again that w and its derivatives tend to zero for 


©, whereas at the minimum cross-section 


] ) 
z = 0: = 0, q = 0. (17) 
dx 


Trying again a power series solution in / ', we obtain the simple solution 


2(1 — v)(1 + 2r) n ] 
——, £ — ———, . 8 
5 — 4° E6+2°/2R (18) 


Ww -~ 


which satisfies both the differential equation and the boundary conditions. 
The bending stress at the inner surface of our shell is given by 
| i d°w n | 2v 1 — 3° 
al 2 eDe na _s <i 2» (19) 
Z21—v dz 6R 5 — 4v [1 + (€/2)] 


where € is defined by (9). The maximum stress occurs at the inner surface in the minimum 


section 
n 1 + 2p 
O max °D l + ~ St ee (20) 
bhi a — 4p 
the nominal stress is given by n/6R, and the stress concentration factor is therefore 
- (6 — 41 +») 
5 aa (21) 
2 — 4p 
varying from 1.2 for pv 0 to 15 for» = 6.5. 


4. Alternative approach. The stress distribution depends on the two coordinates 


and z in Fig. 2. Putting 
x = (R5)'”’é, 2 — 09, (22) 


we may rewrite the basic equations and boundary conditions (the generalized plane 
stress equations for the strip problem, and the axially symmetric equations for the 
cylindrical bar problem) in terms of the new independent variables £ and n. Developing 
these equations with respect to the small parameter « = 6/R, and assuming a similar 
development for the stresses, a development for the stress concentration factor may be 
obtained with respect to e. It may well be that such a development does not converge 
for any finite value of « but even then it will represent an asymptolic expansion. The 
first term has already been obtained by this method for the strip problem, with the 
result 


kK = 2+ 0(, (23) 


confirming our elementary result of See. 2 for « — 0. Further terms in the strip problem 
and the similar development for the cylindrical bar problem are now being investigated. 


'The author is indebted to his collaborator Mr. J. G. Lekkerkerker who pointed out an error in 





the original analysis 
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THE ANISOTROPIC TENSORS* 
By G. F. SMITH anp R. 8S. RIVLIN (Brown University) 


1. Introduction. We consider an nth order Cartesian tensor with components 


Opp ity Uh, 1, 2, 3) in the rectangular Cartesian coordinate system zx; (¢ = 1, 2, 3) 


and components a* 
zx; and 2* are related by the orthogonal transformation 


in the rectangular Cartesian coordinate system x* , where 


r* S ? (1.1) 
with 
$..§ bau (1.2) 
where 6;, is the Kronecker delta. 
Then, 
a* = 8,.;.8 S;,;,a 1.3) 
If 
" ) SEEe Se eee ee 1.4) 
for all s;; satisfying (1.2), then a,,;,...;, are said to be the components of an isotropic 


tensor. If the relation (1.4) is valid only for a subgroup {T} of the group of transfor- 
mations defined by (1.1) and (1.2) we shall describe a,,;,...;, as an anisotropic tensor. 
Q;,i,---;, May then be described as invariant under the group of transformations. The 
isotropic tensor is, of course, invariant under the orthogonal group. 
If (1.3) and (1.4) are valid for a transformation s;; , then it follows that 
Qisis = &;,;,8 eigen erry oer 1.5) 


Seng, i 2 inin 


If 83; satisfies the relation (1.2), then we readily see from (1.5) that 


Bisissertn = 8i55:8i nin °° ° Sintec Bisiare’ 
, j inQjyjaeesi (1.6) 


It is shown in Sec. 3 that any tensor which is invariant under the group of transfor- 
mations {T} may be expressed as the sum of a number of terms, formed from the outer 
products of a finite set of tensors, with scalar? coefficients. This finite set of tensors, 


*Received June 27, 1956. The results presented in this paper were obtained in the course of research 
sponsored by the Office of Ordnance Research, U. 8S. Army, under Contract No. DA-19-020-3487. 

tThroughout this paper, we shall employ the term scalar in the following sense. If corresponding 
quantities g and ¢* are defined in the coordinate system z; and in each coordinate system zf , into 
which z; is transformed by the group of transformations {T}, then if y* = ¢, we shall say that ¢ is 
scalar with respect to the transformation group {T} or, more briefly, scalar. 
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each of which is invariant under the group of transformations {T}, is described as a 
tensor basis for the group of transformations {T}. 

It is also shown in Sec. 3 how such a tensor basis can be found for any group of 
transformations {T} for which we can determine a polynomial basis for polynomials 
in the components of an arbitrary number of vectors, which are form-invariant under 
the group {T}. In Sees. 4, 5 and 6 examples of such tensor bases are obtained for the 
monoclinic-domatic and rhombic-pyramidal classes of crystal symmetry and for the 
case of transverse isotropy. 

The results in this paper are obtained for three-dimensional space. However, it is 
immediately evident that similar methods can be used to determine the isotropic or 
anisotropic tensors for a space of arbitrary dimensions. 

2. Genera! considerations. Let P be any scalar polynomial in the components of 


n vectors uS?, ul, --- , u, which is form-invariant under a sub-group {T} of the 
group of orthogonal transformations defined by (1.1) and (1.2). Then, if 
(r) (r) 
u* = t; ju; (r -_ 1, 2; viache nN), (2.1) 


where ¢;; is any transformation of {T}, we have 


(1) (2) (n) (1) * (2) (n) f¢ 
P(u;'’, Mey °** pe) = PP” ut, oo td. (2.2) 


? 
There exists a finite polynomial basis for the polynomials P. Let J, , Jz , --- , Iw 
denote this polynomial basis. Then, 


Pius”, us”, --- us”) = QU, , In, «+: , In), (2.3) 


? 


where Q is a polynomial in the indicated variables. 
Let us consider the scalar polynomial 


(1), (2) ( 
P @ @p.4,:49e te °° he (2.4) 
in the components of the n vectors uf”, uf”, --- , uj”, the coefficients a,,,,...;, in which 
tisfy the relation 
Qivigretin -_ lisisbisis ne bintnBisiee**in (2.5) 
ois bs ibiste bininBistgserin 


for every transformation ¢;; of the group {T}. 
It is readily seen that P is form-invariant under the transformations of the group 


{T!. From (2.1), we have 


* (1), #2) x(n) (1), (2) a, «™) f 
Qi igrerinti, Ui, eh Ne Us, — Qinigess tsiistisis “ Cinintlis Uj, Poa Uj, a (2.6) 
From (2.6) and (2.5), we obtain the relation 
% (1), (2) oo (1), (2) (n) ~ 
Birigeoog UE UE 88s UE me Be, cgeoceelle, Min °° * Wis » (2.7) 


stating the form-invariance of P. 

Since they satisfy the relations (2.5), the 3" quantities a,,;,...;,, may be regarded as 
the components of a tensor, invariant under the group of transformations {T}, in the 
coordinate system x; and in the coordinate systems into which 2; is transformed by the 
transformations of the group {T}. We thus obtain the result: a scalar polynomial in 
the components of n vectors, linear in each of the vectors, in which the coefficients are 
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tensors invariant under the transformations of the group {T}, is form-invariant under 
the transformations of the group {T}. 

This conclusion can be readily generalized to yield the result that any scalar poly- 
nomial in the components of any number of tensors of any rank, in which the coefficients 
are tensors invariant under the transformation group {T}, is form-invariant under the 
transformations of the group {T}. 

Now let us consider a polynomial P in the components of n vectors —. «© ort, 


us" _ of the form 


P me bj.g00058i, Hin °° Mee (2.8) 
which is form-invariant under the group of transformations {T}. Then, 
a Se ere A He) 2.2 ge (2.9) 
Employing the relation (2.1) in (2.9), we obtain 
Deg. cocg thee Ue, °° * Wen De igseecabess tects t;.;.;, Uz, *°° Uz. (2.10) 
Whence, 
Se ees ee ee ee (2.11) 


We thus see that the coefficients in a homogeneous scalar polynomial of degree n in 
the components of n vectors, which is linear in each of the vectors and form-invariant 
under the transformations of the group {T}, are the components of a tensor of rank n, 
invariant under the group of transformations {T}. 

This conclusion can be readily generalized to yield the result that any scalar poly- 
nomial in the components of any number of tensors of any rank, which is form-invariant 
under the transformations of the group {T}, has coefficients which are tensors, invariant 
under the group of transformations {T}. 

3. The tensor basis. Since the polynomial P defined by (2.8) is form-invariant 
under the group of transformations {T}, it must be expressible as a polynomial in the 


elements of a polynomial basis for polynomials in the elements of n vectors which are 


form-invariant under {T}. Also, since P is linear in each of these vectors, the elements 
of the polynomial basis which are non-linear in any of the vectors need not be con- 
sidered. Let Jd), Js, °** ; J p be the elements of the polynomial basis which are linear 
or of degree zero in each of the vectors. Let us suppose that 


ek | ae © Mere ae 3.1) 


i 


Then we rear ily sec that 


o”"P 


D ees ee ee im) 
Ou; OUu;.” +++ OU;. : 
si is (3.2) 
- a (JpJe see Jy) 
Ou; Ou;, *°* OU;. 
The product J,.J/, J. oceurring in each of the terms under the summation sign in 
(3.1) is of degree n in the elements of the n vectors u§”, u{”, --- , u;” and is linear in 


each of the vectors. No two of the factors J, , J, , --+ , J: occurring in a single product 
involve elements of the same vector. Let us assume that J, is formed from the vectors 


(1 g ° . . (e+ +2 a) 
Us, the --» , u,” only, that J, is formed from the vectors u;’ 7 a mee Us; 
(s+2 (t) 


only and J, is formed from the vectors uj*""’, uj’"", ++* , Ui only. Then Eq. (3.2) 
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, a a J, a’ J, 
Vivigtttin — SApgeest a. (i)4.G) +a @) a (oti) a (ots) — 2 
rv OUu;, Ou;, $e du,” due ee ja Ous. (3 3) 
ie F 





au PAUL? +++ aus? 

Since the quantities A,,..., must be scalars under the transformations {T}, we see 
that b,,,;,...;, must be given as the sum of terms formed from the outer products of a 
number of tensors of the type 0”J,/du{’” dui? --- dui, obtained by differentiating 
the elements J, , J2, --- , Jz of the polynomial basis with respect to the vectors from 
which they are formed. The coefficients of these terms are scalars with respect to the 
transformation group {T}. 

4. The monoclinic system—domatic class. Monoclinic symmetry may be described 
with relation to three preferred directions in space. We denote these directions by the 
unit vectors h, , h, and h, . The vectors h, and h, are not at right-angles and the vector 
h, is perpendicular to the plane defined by h, and h, . 

Let us choose as reference frame a rectangular Cartesian coordinate system x, , the 
x, axis of which coincides in direction with the vector h, . The axes x, and x; may be in 
arbitrary perpendicular directions in the h,h, plane. Then, the group of transformations 
{T} associated with the monoclinic-domatic symmetry class consists of the identity 
transformation T, and the transformation T, defined by 


{—1, 0, 0) 
T,=| 0, 1, 0 . (4.1) 
L 0, 0, 1) 
Let us consider n vectors with components us? us, +++ , uf? in the coordinate 


system x; . Any polynomial 


(1) (2) (2 


(1) (1) ) (2) (n) (n) (n) 
Pas, ie ste 5 ee ee oe ee 


> 


is unaltered by the transformation T, and is transformed into 


(2 (2) 


/ 1) ( (2) = ane (n) (n) (n) 
P(—u,’, tle» Uy» ~My ye Mey oO ** » yp he fp Me: 


by the transformation T, . The necessary and sufficient condition that P be form- 
invariant under the transformation T, is that it be expressible as a polynomial in 


r) (r) 


— = i (7,3 = 1,2, --- ,n). (4.2) 


This result follows immediately from the first main theorem of classical invariant theoryT, 
according to which a polynomial basis for polynomials in &, , 7, (r = 1, 2, --- , n), which 
are form-invariant under interchange of £, and 7, , is given by £, + 7, and é,n, + &,”, . 
Taking 


1) (1) (1) (n) (n) (n) 
» > °°* 9 &) mt (ey Me Me 5 o** 9 Me Me te 


Ser 
ier 


tSee, for example, H. Weyl, The classical groups, their invariants and representations, Princeton 
Univ. Press, Princeton, N. J., 1946, p. 36 et seq. 
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“ 


and 

0). 18)... 08) (n) . (mn) 2. (n) 
AM is Mag 8s n.) = (—w 9 U2 »Uzg 5 °** y —U;", Us » U3 ) 
and neglecting redundant elements, we obtain the result (4.2). 

The expressions (4.2) form a polynomial basis for polynomials in the components of 
the n vectors u‘” (r = 1, 2, --- , n) which are form-invariant under the transformations 
characterizing the symmetry class considered. From the results of the previous section, 
we see that the quantities dus” /du’”, dus” /dus” and a? (uf? ul”) /dulf? dul” (r * s) 
define a basic set of anisotropic tensors for the monoclinic-domatic class in terms of 


r) 


which an arbitrary anisotropic tensor for this class can be expressed. We note that 
r) r) lal . 
dus” /dus” defines the vector with components 


5, = (0, 1,0) = (a) (say) (4.3) 
in the coordinate system x; , that dus” /du{” defines the vector with components 
6;; = (0,0, 1) = (6,) (say) (4.4) 


in the coordinate system 2x; and 0°(uj” uj”) /duj”du;” (r # s) defines the tensor of rank 


2 with components 


1, 0, O 
10, 0, 0} = @;) (ay) (4.5) 
0, 0, O 
in the coordinate system x; . Further, any tensor with components a; ;..., in the coordinate 


system 2; , which is invariant under the transformations of the monoclinic-domatic 
class, is expressible as the sum of terms formed from outer products of a; , 8; and ¥;; 
with scalar coefficients which are invariant under the transformations of the class. 

5. The rhombic crystal system—rhombic-pyramidal class. Rhombic symmetry 
may be described with relation to three preferred directions in space, defined by the 
mutually perpendicular unit vectors h, , h, , h; . We choose as our reference system the 
rectangular Cartesian coordinate system 2x; , the axes of which coincide in direction 
with the vectors h; . The group of transformations {T} associated with the rhombic- 
pyramidal class is composed of the identity transformation T, and the transformations 


y . :¥ and c# defined by 


0, 0 il, 0, O 1, O, 0 
T= /|0, -—1, Ol, T,= /|0, 1, 0 and T,= |0, -—1, 0}. (5.1) 
0, oO, 1 = = . £ wf 


We can determine, in a manner similar to that adopted in Sec. 4, a polynomial basis 
for polynomials P in the components of n vectors which are form-invariant under the 


: fm ‘ (1) (2) ( : 
transformations (5.1). If these vectors have components u{”’, uj”, --- , u;” in the co- 
ordinate system x; , then invariance under the transformation T; implies that P must 

. . Li - (r) ( (r ( P 
be expressible as a polynomial P’ (say) in the quantities uj”, u3”, u,” uz” (r, s = 1, 2, 


--» .n). Thus, 


P = Put”, us”, us” uy”). (5.2) 
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Invariance of P’ under the transformation T, imposes on P’ the limitation 


P’(ui”, us? , usPus”) = Pus”, —u$”, us” us’). (5.3) 


It is readily seen that P’ must be expressible as a polynomial in the quantities u{”, 


us”? us? and us” us’. Such a polynomial is obviously invariant under the transformation 
T, and hence under all the transformations of the group {T}. 

Following the procedure described in Sec. 3, we see that du{”/dus” (= a; , say), 
d°(us” us”) /du§” dus”? (r # 8) (= Bi; , Say), O°(us” us”)/dus”? dus” (r ¥ 8) (= ya; , Say) 
define a basic set of anisotropic tensors for the rhombic-pyramidal class. We have 


(a;) = (6,,) = (1, 0, 0) 


0, 0, 0} (0, 0, 0 
(8;;) = |0, 1, O| and (y%;) = |0, 0, 0} (5.4) 
lo, 0, 0) lo, 0, 1 


6. Transverse isotropy. Transverse isotropy may be described with relation to a 
single preferred direction in space defined by a unit vector h. We choose as our reference 
frame a rectangular Cartesian coordinate system x; , the x3; axis of which coincides in 
direction with the vector h. Then, the group of transformations {T} associated with 
transverse isotropy with respect to the direction of h is composed of the transformation 
T., defined by (4.1), the transformation T, defined by 


cosw, sinw, 0} 


T. = |—sinw, cosw, 0| (0 <w < 2n), (6.1) 


SS: By 1 


and the transformation T,T, . 

We can determine a polynomial basis for polynomials P in the components of n 
vectors which are linear in each of the vectors and form-invariant under the trans- 
formations T, , T, and T,T, . Let u{"’, uf”, --- , us” be the components of the n vectors 
in the coordinate system x; . Then, the limitations imposed on P by the condition that 
it be form-invariant under the transformations T, , T, and T,T, are 


»/,, (7 . 1, (r) (r) os ( . (r) 
P(ul?, uf?, uf?) = Pluk? cosw + uf” sinw, — uj” sinw + us” cosa, us”) ’ 
(6.2) 
r r (r) \ 
P(—uy ’ us , U3 )- 
It is seen that the dependence of P on uS” gives rise to no restrictions on the form of 
3 
P and, therefore, a polynomial basis for P is formed by us” (r = 1, 2, --+ , mn) and the 
polynomial basis for polynomials P’ which satisfy the conditions 


) 


, r r) , (r) (r) (r ° (r 
P’(ul?, us?) = P’(uf” cosw + uz” sinw, — uj” sinw + u;” cos w) 


= P’(—uj;”, uz”), 


and are linear in the two-dimensional vectors u“” (a = 1, 2). The relations (6.3) imply 
that P’ is form-invariant under the group of two-dimensional orthogonal transformations, 
proper and improper. Now, the quantities uf? uS” (r,s = 1, 2, --- ,n;a@ = 1, 2) forma 


polynomial basis for polynomials in the components of n two-dimensional vectors uf” 


(6.3) 
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which are form-invariant under the full orthogonal groupt. Since P’ is linear in the 


vectors, we can omit those elements uu‘ u” for which r s. We thus have the result 
that the quantities 
r) s) : . 
5 the Ua (r, 8 Ly 2, *** 2297 2) (6.4) 


form a polynomial basis for polynomials in the components of n vectors u;" which are 
linear in each of the vectors and form-invariant under the transformations associated 
with transverse isotropy about the zx, axis. 


As in the previous sections, we immediately see that a basic set of anisotropic tensors 


for transverse isotropy about the x, axis is defined by 0u;"/du{” (= a; , say) and 
a(us”’ Uy /Ou;" du B:;) (1 F sya 1, 2;2,7 = 1, 2,3). We obtain immediately 
(a;) = (6;;) = (0, 0, 1) 
and 
I 6, © 
(B, (6:;6:; + 52:6.;) = |0, 1, O}. (6.5) 
0. 0, © 


ON A PRINCIPLE OF RECIPROCITY BETWEEN HIGH- AND LOW-FREQUENCY 
PROBLEMS CONCERNING LINEAR DIFFERENTIAL EQUATIONS OF 
SECOND ORDER* 


By AUREL WINTNER (The Johns Hopkins University) 


DP 


following considerations, in which nomenclature and point of view will be 


1. Th 

the same as in two earlier notes’, deal with the oscillatory case of the differential equation 

x’ + fiijxz = 0, (1) 
in which f(t) is a re alued, continuous function for large positive ¢t, and z(t) is any 
real-valued solution distinct from the trivial solution x(t) 0. 

If (1) is oscillatory when f(t) = f,(¢), and if f,( s f(t), then (1) is oscillatory 
when f(? t). This follows from Sturm’s comparison theorem, which implies also 
the following fact: If (1) is oscillatory and if d, = d,(/) d,(f; x) denotes the distance 
ti... — f, , where ? f; x) is the nth zero of a solution x = x(t) = 2,(t) of (1), then 
lim sup d,(f;)/d,(f2) < 1, where n — ~, holds whenever /f,(¢) — f,(¢) exceeds a positive 
lower bound as ¢ - 

3y choosing f,(4) = w, where w is a positive constant, and letting w > © or w — 0, 
it follows that d (f) >JOasn—- wolf 


fow-o- as too, (2) 





+See, for example, H. Weyl, loc. cit. p. 52 et. seq. 


*Received July 5, 1956. 


1A, Wintner, Quart. Appl. Math. 7, 115-117 (1949) and 13, 192-195 (1955). These two papers will 


be referred to as [1] and [2] respectively. 
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and that d,(f{) ~ - asn— o if 
fij-0 as too (3) 
(provided that (1) is oscillatory; this additional proviso is necessary, even if 


S > 0, (4) 


in the case (3), but is automatically satisfied in the case (2), or even if just f(t) = w” 


holds for some positive constant.) For this reason, (1) will be called of high frequency 
in the (oscillatory subease of the) case (3). 
2. If f(d) \/t (for 1 S t < @), where Xd is a positive constant, then (4) and (3) 


[ toa=« (5) 


fails to hold). But it depends on the numerical value of \ (hence, not merely on the 
order of magnitude of f() > O for large #) whether (1) is or is not oscillatory in this 
case. In fact, the remark of A. Kneser, quoted in® [2], shows that the case {(#) = /@° of 
(1) is oscillatory only if \ > 1/4, rather than just \ > 0, is assumed. 

Thus (3) and the negation of (5) together are compatible with the oscillatory as well 
as with the non-oscillatory character of (1) in the case (4). On the other hand, (1) must 
be oscillatory if (4) and (5) hold; in fact, more than this was proved in [1]. 

In what follows, both (4) and (5) will be assumed. This will imply, in particular, that 


are satisfied (and 


(1) is oscillatory even if (3) holds. 

3. There is an extensive literature” dealing with the asymptotic behavior (as t > @) 
of the solutions x(¢) of (1) in the case (2) of high frequency. Let any such theorem be 
referred to as of type (//), and let (L) refer to theorems dealing with the asymptotic 
behavior of the solutions x(t) of (1) in the case (3) of low frequency. 

An inspection of the literature reveals a certain empirical dualism between these 
two types of theorems: As soon as a theorem of type (/) is discovered, its wording and 
proof can usually be altered so as to supply a theorem of type (1), and conversely. But 
as far as the literature consulted goes, this principle of transfer appears only ad hoc, 
since its applicability is not assured a priort. 

The following consideration will dispose of this anomaly, by exhibiting a formal 
fact which underlies the principle of tranfer, (H) — (L) or (L) — (A). 

4. It will turn out that what actually matters in this regard is neither (2) nor (3) 
but a connection between the differential equations (1) and 


+ g(s)x = 0, (6) 


where the dots, in contrast to the primes of (1), refer to differentiations with respect to 


the independent variable s = s(t) defined by 


ds = f(t) dt, (7) 








2Concerning the brackets see the preceding footnote. 

‘cf. G. Sansone, Equazioni differenziali nel campo reale, vol. 2 (2nd ed., 1949), chap. VII, and F. G. 
Tricomi, Equazioni differenziali (2nd ed., 1953), chap. III-IV. These expositions, and the papers men- 
tioned in the first footnote above, contain further references to the literature. 
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( 


and the coefficient functions of (1) and (6) are reciprocals by virtue of (7) or 


dt ds /g(s), (8) 

that is, 
bi l/g or g 1/f. (9) 
Note that, because of (4) and (7), the correspondence s = s(¢) between the given 


t-range and the s-range is a one-to-one correspondence, which is continuously differ- 


entiable along with its inverse, ¢ = ¢(s), and that, by virtue of (4), (7) and (5), 
t— o ifandonlyif s— o, (10) 


It is clear from (9) that (6) is of low frequency or of high frequency according as (1) 
is of high frequency or of low frequency. But the connection between (1) and (6), being 
of a formal nature, will be such as to apply also in the intermediary cases, cases in which 
neither (2) nor (3) is satisfied. 

5. The proof depends, as in [1], on an appropriate application of Riccati’s equation 
(an application which, however, cannot be the same as in [1]; ef. Section 6 below). 

First, if J is a ¢-interval on which a solution x(t) of (1) does not vanish, then it is 
possible to form the logarithmic derivative 

m(t) = x’(t)/x(t) (11) 


at every point ¢ of J, and m = m(t) will be a solution of Riccati’s differential equation 


m’ + m + f(t) = 0. (12) 
In fact, if (1) is divided by 2, then (12) follows from (11). 
Next, suppose that not only x(t) but also 2’(¢) is distinct from 0 at every point ¢ of J. 
Then (11) shows that it is possible to form on J the function 


n(t) = —1/m(s), (13) 
and it is seen from (13) that (12) is equivalent to 
n’ +1+ f(in® = 0. (14) 


Finally, if (14) is divided by f(é) (which, in view of (4), is allowed), then the defi- 


nitions, (9) and (7), of s and g show that (14), where n’ = dn/dt, can be written in the 


form 
i +n? + os) = 0, (15) 


where n = dn/ds. 

This proves that, as a matter of formal coincidence, (1) is equivalent to (6) by virtue 
of (9) and (7). In fact, (15) belongs to (6) in the same way (that is, as the Riccati resolvent 
of the linear differential equation) as (12) belongs to (1). 

6. Actually, the proof is not yet complete. For, on the one hand, Section 5 confines 
t to a t-range J on which neither z(t) nor x’(t) has a zero and, on the other hand, the 
problems of the asymptotic behavior concern themselves with the ¢-range I which is 
a half-line 4; < ¢ < © on which x(t) (hence z’(¢) as well) has an infinity of zeros, since 
(1) is oscillatory. In fact, the oscillatory character of (1) is implied by (4) and (5), and 
(5) cannot be omitted, since otherwise (10) becomes false. This difficulty can however 


readily be disposed of, as follows: 
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Instead of a single real-valued solution x(t) of (1), choose for (1) two solutions, z = x(t) 
and x y(t), which are real-valued and linearly independent. In view of the super- 
position principle, the (complex-valued) function 

2(t) = x(t) + iy(t) (16) 
is a solution of the (real) differential equation (1). But 
ot) XO and 2(t) #0 (17) 


now hold on J when J is the entire half-line under consideration (47 S tf < ). In fact, 
since x = x(t) and x = y(t) are linearly independent solutions of (1), their Wronskian 
— yx’ is a non-vanishing constant; hence neither x and y nor x’ and y’ can vanish 


c "g 
at the same é, and so it is clear from (16) that both inequalities (17) are satisfied at 
every lf. 

Owing to the first of the inequalities (17), the logarithmic derivative 


m(t) = 2’(t)/2(t) (18) 


exists on the entire half-line 7. On the other hand, the second of the inequalities (17) 
assures that (13) defines the function n(¢) on the entire half-line J if m(t) is defined by (18). 
Finally, since x = z(t) is a solution of (1), Riccati’s deduction of (12) from (11) and 
x(t) of (1) remains unaltered if (11) and (1) are replaced by (17) and the 
2(t) of (1) respectively. Since the transition from (12) to (15) had nothing to 


the case x 
case r 


do with the reality of m(t), the proof is now complete. 


ON THE ROTATION OF AN INCOMPRESSIBLE CONTINUOUS 
MEDIUM IN PLANE MOTION* 


By WALTER NOLL (Carnegie Institute of Technology) 


We consider a plane motion of an incompressible medium with homogeneous density 
p = const. We assume that the plane region occupied by the medium is of the following 
type: at each instant it is bounded by a finite number of smooth curves. If it is finite 
all but one of the closed boundary curves move rigidly. If it extends to infinity all finite 
closed boundary curves move rigidly. We label this motion by I. Now we consider 
another motion, labeled by II, which is obtained from I by superposing a uniform 
rotation with angular velocity # = const. about a fixed axis perpendicular to the plane 
of the motion. 

It will be shown that motion II is dynamically possible whenever I is, and that the 
stresses in II differ from those in I only by a scalar stress given explicitly in (10) below. 

G. I. Taylor [1] proved this theorem for the special case of the irrotational motion 
of an inviscid fluid caused by rigid cylindrical bodies moving in an infinite container. 
[It has been shown recently that Taylor’s result continues to hold for linear viscous 
fluids [2] and for a special type of non-linear viscous fluids [4]. The main purpose of 
this paper is to point out that Taylor’s theorem remains valid for any material whatsoever 
(it may even be an anisotropic solid). The essential assumptions are that the medium 


be incompressible and that the motion be plane. 


*Received August 17, 1956. Written under Contract Nonr 228(09), Office of Naval Research. 
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A motion with velocity v is dynamically possible if there is a stress tensor S such 
that the continuity equation, the equation of motion, and the constitutive equation 
of the material are satisfied with v and S. The first two of these equations are in our case 
divv = 0 (1) 

div S + pf = pv, (2) 


where f is the density of the extraneous forces. In general, the constitutive equation will 
serve to determine the stress S by the motion (and possibly certain initial data). How- 


ever, in an incompressible medium the motion determines the stress only up to a scalar 


stress pl (/ unit tensor) so that only the extra stress 7 = S + pl is completely 
determined. It follows from the principle of isotropy of space (see [3], p. 19) that the 
constitutive equation is not altered by the superposition of any rigid rotation. 


We describe motion ITI in a frame of reference which rotates with the superposed 
rotation. Then both I and II will have the same velocity field v,; = v,; = v. Therefore, 
by the principle just mentioned, the extra stresses must be the same for both motions: 


Tl, = pl + 8, Pul + Si. (3) 

The equation of continuity (1) is the same for both motions. For a plane motion, (1) 

implies the existence of a stream function y such that 

V — cur! (ey) (4) 

where e is a unit vector perpendicular to the plane of motion. The assumptions on the 

region mentioned in the beginning ensure that y is single valued. The equation of motion 
(2) for I is 


div Si + pf = pv (5) 


In the rotating frame, however, we must add to f a field g which represents the centrif- 
ugal and Coriolis forces, so that we get 
div Sy 1 p(f -f- g£) = py (6) 
(It is understood that the extraneous forces f are to be rotated also). It follows from 
(3), (5), and (6) that motion II will be dynamically possible if and only if p,; can be 
chosen in such a way that 
pg = grad (pir — pr) (7) 
holds, i. e., if and only if g has a single valued potential. This is not the case for a general 
three-dimensional motion. But for a plane motion an easy analysis (see [1], p. 101) 


shows that 
g = grad (Qwy + 4w’r’), (8) 


where r is the distance from the axis of rotation. It follows from (8), (7) and (3) that 
the pressure p,; and the stress S,,; of motion II are given by 
Pu = Di + pl(2wy + 4w'r*), (9) 
Si = Sy oe p(2wy + Ler) I. (10) 


‘i : : ‘ . J . - : 
From (10), the forces exerted on the boundaries in motion II_can be obtained in the 


same way as was done by Taylor [1]. 
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BOOK REVIEWS 


(Continued from p. 282) 


The theory of functions of real variables. By Lawrence M. Graves. McGraw-Hill Book 
Co., Inc., New York, Toronto, London, 1956. xii + 375 pp. $7.50. 


Published ten years ago, the first edition of this book is now widely known and has been reviewed 
in detail: see, e.g., Math. Rev., 8 (1947), 319-20, or Bull. Amer. Math. Soc., 54 (1948), 487-9. Except 
for a few rewritten passages, the second edition differs from the first only in that two new chapters have 
been appended, and we shall confine comment to these. 

Chapter XIII, on sets and transfinite numbers, begins with an amplified account of sets and functions 
(= relations) from the non-formal ‘‘common-sense”’ point of view adopted by the author in Chapter I. 
In the same spirit there follows a discussion of order leading up to definitions of ordinal and cardinal 
number which are essentially von Neumann’s: an ordinal is a set which (i) is well-ordered by the class- 
membership relation and (ii) has each of its members as a subset; and a cardinal is an initial ordinal 
in a cardinal-similarity class of ordinals. Next, the maximality principle (Zorn’s Lemma) is stated and 
studied in a variety of equivalent forms. The chapter closes with an introduction to ordinal and cardinal 
algebra. Chapter XIV discusses, for metric spaces, basic facts concerning convergence, completeness, 
connectedness, separability, compactness, and category. Here the references for further study might 
well have included such recent items as J. L. Kelley’s General Topology, 1955, as well as the older works 
listed. Compared with the rest of the book both of the added chapters are notable for their profusion 


of well selected exercises. 
TRUMAN Bort 


Automata studies. Edited by C. E. Shannon and J. McCarthy. Princeton University 
Press, New Jersey, 1956. viii + 285 pp. $4.00. 


The recent development of high-speed digital computers and of electronic control equipment has 
revived interest in questions on the ultimate limits of a ‘‘mechanical brain” and on how much of animal 


or human mental activity can be described in physical terms. The present collection of studies represents 


some of the contributions of 
The articles on the properties and limitations of computers are based on the work of Turing, in which 
Tt ticl the pro] 1 limitations of computers are | 1 on the work of T . hicl 


mathematics to these perennially interesting fields of speculation. 


he showed that the action of most digital computers is equivalent to an appropriate sequence of actions 
5 1 I 1 
of a machine with simply describable properties, the limitations of which could be studied by known 


mathematical methods. The 
carry further the analysis of their behavior. 

The papers more specifically related to the analysis of mental action are concerned with the properties 
lilt up of ‘‘neurons” having a simple, all-or-none behavior analogous to living 


rticles in this collection presuppose a knowledge of Turing machines and 


of “automata,” machines | 
nerve cells, The first paper of the collection, by S. C. Keene, gives a readable introduction to the subject. 
Some of the later papers discuss limitations and extensions of the basic ideas. 

Of course an automaton, built up of infallible “neurons” could perhaps duplicate some of the 
human brain’s abilities in the development of formal logic. To duplicate actual learning and memory 
we have to introduce the possibility of failure of an element from time to time. The study of automata 
with unreliable elements is of considerable interest, for it may point the way to methods of improving 
is a whole, whether it be one of our all-too-fallible present day computers 


the reliability of the machine 
urticle on this subject, in the collection, is by von Neumann. One can hope 


or our own brains. The basi: 
that others will explore further the ideas he outlines here. 

The articles are organized into three sections: the first on Finite Automata, the next on Turing 
Machines and the last on Synthesis of Automata. The papers are recommended to anyone interested 
in this general field, though some of them make pretty hard reading. If the interesting Introduction, by 
John McCarthy and Claude Shannon, had been expanded by a dozen pages or so to provide a review of 
of previous work and a bibliography of important papers, this volume would have constituted a much- 
As it is, the collection is useful and thought-provoking. 


needed review of the field 
Puitip M. Morse 
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Elementary differential equations. By William Ted Martin and Eric Reissner. Addison- 
Wesley Publishing Co., Inc., Cambridge, Mass., 1956. xi + 260 pp. $5.50. 


This book is intended primarily as a text book for advanced undergraduate science and engineering 
students. The book contains six chapters relating to ordinary differential equations, one chapter on 
finite difference equations and one chapter on partial differential equations. Emphasis is on methods of 
solution for which the authors attempt to give some motivation although this is clearly rather difficult 
at times. There are numerous examples and problems including several examples to illustrate that the 
method being discussed does not always succeed. These latter examples show to the student better than 


words could describe that the subject goes much deeper. 
Gorpon F,. NEWELL 


Advances in mechanics: Volume 1V. Editors: H. L. Dryden, Th. von Karman. Managing 
iditor: G. Kuerti. Associated [ditors: F. H. Van den Dungen, L. Howarth, J. 
Pérés. Academic Press, Inc., New York, 1956. x + 413 pp. $10.00. 


volume is the first to appear since the death of Richard von Mises, the initiator of 
this series; it contains the following articles: The Turbulent Boundary Layer, by Francis H. Clauser 
51 pp Nonlinear Elasticity, by T. C. Doyle and J. L. Ericksen (63 pp.)—Physical and Statistical 
Aspects of Fatigue, by A. M. Freudenthal and E. J. Gumbel (42 pp.)—Three-Dimensional Boundary 
Layer Theory, by Franklin K. Moore (70 pp.)—Dislocation Theory of Plasticity of Methals, by G. 
Schoeck (51 pp.)—The Poinearé-Lighthill-Kuo Method, by H. S. Tsien (69 pp.)—On the Concept of 


Elastic Stability, by Hans Ziegler (53 pp.). 


The prest nt 


W. PRAGER 


Les Epreuves sur Iechantillon. By Maurice Dumas. Centre d’Etudes Mathematiques 
en vue des applications, of the Centre National de la Récherche Scientifique, Paris, 


1955. 170 pp. $2.90. 


This is a short monograph on sampling inspection intended for scientists and engineers who are 
supposed to have the mathematical background, but not the available time, to read up on the relevant 
theory. The presentation is thought-provoking, and perhaps a little complicated. Certainly it is far more 
sophisticated than that of the numerous American books and manuals on practical statistical quality 
control which have been appearing in profusion lately. Incidentally, the preface is by no less a celebrity 
than Professor Maurice Fréchet. 

The author arranges his book in this way: First there is a summary in a telegraphic style of the 
main formulas of the book and procedures for applying them. Then there is a careful chapter on the 
terminology, techniques and rationale of sampling inspection. This chapter includes a brief but clear 
non-mathematical discussion of the relations between statistics, practical probabilities, and mathe- 
matical theory of probability. After this comes a general description, with numerical examples, of the 
four categories of problems which the author intends to discuss: Estimation, acceptance inspection, 


comparison of lots with each other and with fixed standards, and statistical quality control in the 


Shewhart sense. 

In the next chapter (we have now arrived at Chapter III), detailed procedures for attacking these 
four types of problems are given for the case of sampling by attributes. In the last chapter, Chapter IV, 
methods for the case of sampling by variables are given. The idea of the arrangement is that a reader 
who wants only to get a general idea of what sampling inspection is all about, can do it by reading only 
the first part of the book, omitting the last two chapters. 

Mathematical notations and formulas are used liberally, but there are no proofs. For these, the 
reader is referred chiefly to two large works which stand behind this one: a six-hundred and thirty page 
work on industrial statistical methods written by the author in collaboration with P. Maheu, and a 
work called Technique modernes de contréle des fabrications, by J. Mothes. There is a paucity of ac- 
knowledgments to those whom the auther calls “the Anglo-Saxons’’, although much of the terminology 
is drawn directly from Fisher, Dodge and Romig, Neyman and Pearson, and Shewhart. (By the way, 
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bo 
tb 


Dr. Shewhart’s name is apparently never mentioned.) This is certainly excusable in view of the nature 
of the book, but Dumas indicates in his brief bibliography that he has another reason for not sensing 
a debt to the Anglo-Saxons too keenly. He says that anyhow, in a 1925 Mémorial de l’Artillerie Francaise 
entitled Sur une interprétation des conditions de recette, he anticipated the relevant portions of the Dodge- 
Romig and Neyman-Pearson theories. 

The author generally limits himself to strictly classical procedures and formulas. There are short 
sections on sequential analysis, and some of the tests of Barnard are given, but otherwise this book 
st twenty-five years ago. There is, for example, nothing about non-para- 


4 


could have been written at k 
metric analysis and nothing about sampling by variables when lot quality is measured by a percentag 


of non-conforming items. Perhaps the author should read more Anglo-Saxon articles on mathematical 


statistics. 
However there is an interesting emphasis on a posteriori probability, which this reviewer considers 
inspection from a pragmatic rather than from a philosophical point of 





to be of importance in samplir 
view. Almost all sampling inspection is performed on lots which themselves can be considered as prob- 
ability samples from a larger population, and one way to take this into account is to use Bayes’ Theroem. 
The essentially stochastic nat 
current Anglo-Saxon literature of the subject. 

In the case of sam] by variables, the author deals at some length with the technique of plotting 


observed frequency distributions on normal probability paper. He calls such a plot a “graphique de 


yf inspection lots does not seem to be sufficiently emphasized in the 





rea 


Henry’’. In this connection, readers of the monograph under review might do well to look at an article 
entitled ‘‘Use of Normal Probability Paper”, by Chernoff and Lieberman, in the Journal of the American 
Statistical Association, vol. 49, 1954, pp. 778-785, which discusses the question of optimal choices of 
1 
the ordinates for various problems. 
All in all, this book should be quite successful in meeting its aims, especially as regards French 
readers. In fact, some Anglo-Saxons who meet the specifications of the intended audience, and who can 


read the rather turgid prose which seems to be in favor nowadays in French scientific circles (many 


semi-colons and colons) might find the book useful. 


J. H. Curtiss 


Zur Ermittlung von Unterschallstrémungen mit der Transformationsmethode bet quad- 
ratischer Approximation der Adiabate. By Hans Schubert and Erich Schincke. 
Akademie-Verlag, Berlin, 1955. 32 pp. $.55. 

treatment of the problem of finding the compressible inviscid flow past 


} 


This monograph contains 
| object when the pressure-density relation is approximated by such a 


a symmetrical two-dimensiona 

function that the equation for the dimensionless stream function y becomes Ay + k*y = 0. Here k isa 
constant and the independent variables are distorted hodograph variables. It is shown that this p-v 
law agrees with the usual isentropic relation to quadratic terms in their series in powers of the velocity. 
Solutions are then expected as the sum of a singular term closely related to the singular part of the 
hodograph solution of the corresponding incompressible problem and a series of non-singular terms 
whose coefficients are determined by the boundary conditions in the physical plane and which are 
therefore to be extracted as the solution of a non-linear algebraic problem. The authors note that the 
last step is very laborious 

G. F. CARRIER 


Asymptotic expansions. By A. Erdélyi. Dover Publications, New York, 1956. vi + 108 
pp. $1.35. 
This is an unaltered republication of Technical Report 3, prepared under Contract Nonr-220(11) 
for the Office of Naval Research. 


An Essay on the foundations of geometry. B. A. W. Russell. Dover Publications, New York, 
1956. xvili + 201 pp. $3.25. 


Except for an added foreword by M. Kline (7 pp.), this is an unaltered edition of the 1897 edition. 
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Mechanik der Kontinua. By G. Hamel. Edited by I. Szabé. B. G. Teubner, Stuttgart, 
1956. 210 pp. $7.10. 

In spite of its sweeping title, the book is primarily concerned with mechanics of fluids. The first 
part (114 pp.) is devoted to perfect fluids and contains chapters on foundations, one-dimensional motions 
of air (sound, detonation), potential flows without free surfaces, vortex motion, steady two-dimensional 
flows of compressible fluids, and potential flows of compressible fluids. The second part (66 pp.) deals 
with viscous fluids and has chapters on the Navier-Stokes equations, slow flow, exact solutions, boundary 
layer theory and turbulence. The third part (19 pp.) has the heading “On more general deformable 
systems’; its only chapter has the title “Elastic vibrations” but is mostly devoted to the plane flow of 
a perfectly plastic solid. Though statements in the editor’s preface would seem to contradict this im- 
pression, the last part strongly suggests that the author’s death interrupted work on a more ambitious 


treatise, leaving us with two completed parts and miscellaneous notes. 


W. PRAGER 


Rayleigh’s principle and its applications to engineering. By G. Temple and W. G. Bickley. 
Dover Publications, Inc., New York, 1956. 152 pp. $1.50 (paper bound). 
This is an unabridged and unaltered republication of the book originally published by the Oxford 
University Press in 1913. 


Proceedings of the conference on high-speed aeronautics. Edited by Antonio Ferri, Nicholas 
J. Hoff, and Paul A. Libby. Polytechnic Institute of Brooklyn, Brooklyn, New York, 
1955. vi + 392 pp. $5.25. 
ollection of papers by distinguished aeronautical scientists given as part of the 

ions of the Polytechnic Institute of Brooklyn. Most contributions are on aerody- 

mportant structural fields are also covered. About one half of the papers are in the 


This volume i 


while the remainder describe new developments. 

imical side von Karman gives a very broad but connected criticism of High Speed 

test stage of development. Of particular interest are his presentation of linearized 
theory in a logical order and his lucid explanation of transonie and hypersonic flow phenomena. A 
second review paper by Dryden covers the complete development of the theory of transition; this 
of spots of turbulence in causing the breakdown of laminar flow. In the group of 
research papers Crocco gives a theory of shock boundary layer interaction applicable to the base pressure 
problem, Oswatitsch and Keune describe a new method of treating transonic flow part bodies and 
Ferri discusses conical fields with unusual mixed flow patterns in transverse planes. There are three papers 
on Hypersonic Flow; two, by Libby and Bloom and by Smelt respectively, describing experimental 
facilities and a third, by Kantrowitz, giving an all too brief survey of physical phenomena. The re- 
maining papers on fluid mechanics are by Busemann on the relation between noise and wave drag, by 
Jones on the efficient operation of high speed aircraft, by Schlicting on Cascade Theory, and by Michel 


] 


€ mphasize s the role 


and Sirieix on Transonic Lift. 
On the structural side Broglio gives a complete account of the balance method in thermal stress 
analysis while Hoff develops a long and clear discussion of creep phenomena. Two further papers, by 


Duberg and Horton respectively, deal with heating effects. 
M. Hour 


Numerical integration of differential equations. By Albert A. Bennett, William E. Milne, 
and Harry Bateman. Dover Publications, Inc., New York, 1956. 108 pp. $1.35 
(paper bound 
This is an unabridged and unaltered republication of a report of the Division of Physical Sciences 

of the National Research Council originally published in 1931. 
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Integral transforms in mathematical physics. By C. J. Tranter. John Wiley & Sons, Inc., 
New York, and Methuen & Co., Ltd., London, 1956. ix + 133 pp. $2.00. 
This is a new edition of a member of the Methuen Monograph Series, first published in 1951 and 
reviewed in Vol. X, No. 2 of this Quarterly. 
The present edition differs from the first in the addition of a chapter on dual integral equations, the 
correction of misprints and an enlargement of the bibliography. 
pter is largely of recent origin and contains many of the author’s 





The work given in the new cha 





own contributions. The equations discussed arise in the solution of boundary-value problems in which 
the condition on one boundary i “mixed” one. It is welcome that an account of transform methods 
for such problems is thus made more widely available. : en 

WaLTER F. FREIBERGER 


The theory of ordinary differential equations. By J. C. Burkill. Oliver and Boyd, Ltd., 
Edinburgh and London, and Interscience Publishers, Inc., New York, 1956. ix GE 101 


pp. $1.55. 


This little text-book is iseful addition to the I niversity Mathematical Texts series. The book 
outlines the theory of ordir differential equations in a concise but lucid manner, with the intention 
e basic results in the most useful rather than the most general form. 


of providing th 
The basic existence theory for first order equations is proved in Chapter I, and Chapter IT treats 


the general linear equation of order n. Chapter III deals with the oscillation theorems on the zeros of 


solutions of the homogeneous second order equation and briefly introduces eigenfunctions and ex- 
pansions. The remainder of the book considers ordinary differential equations, almost exclusively linear 





and of the second order, in complex variables. Chapters IV, V and VI are concerned with solutions in 


series, singularities of equations, and solutions expressed in integrals along a contour. The main results 
and formulae for Legendre and Bessel functions are deve lope din Ch ipters VII and VIII. The last chapter 
gives a brief introduction to asymptotic series and asymptotic solutions of differential equations. A 
ples and problems with answers are given throughout the text 


number of illustrative ex 


Although by no means intended to replace more comprehensive books on the title subject, the 
text should provide a good basis for a second course on differential equations. a 
R. T. SHrevp 


Random processes in automatic control. By J. Halcombe Laning, Jr., and Richard H. 
Battin. McGraw-Hill Book Co., Inc., New York, Toronto, London, 1956. ix + 434 


pp. $10.00 


That this book is an outgrowth of a course given to aeronautical engineers is rather surprising. The 
kground for study of automatic control of aircraft, of fire control, etc. and 


objective is to present a ba 
this direction. The main part of the test, however, describes the subject 


many of the examples lean ir 


trols in quite general terms. Approximately half of the book is devoted to 


of random processes and c 
probability theory and the description of random processes. Although this material is quite concen- 
trated and would require careful reading by someone 
statistics, the presentation is clear, concise and well organized. It is indeed remarkable how the authors 
manage to give such an extensive survey of probability theory in such a short space. The mathematical 
rigor does not approach that of advanced mathematical texts, but the exposition certainly could not 





lacking preliminary training in probability or 


be classed as a popularized version. 
The last half of the book is more applied. There is a chapter on analysis of effects of time-invariant 
cesses (filtering and prediction with noise) and another chapter on 


linear systems on stationary pro 
both of which are well developed subjects. Two 


optimum smoothing and prediction (Wiener theo: 
other chapters deal with generalizations to non-stationary systems and optimum operation with finite 
data. Much of the material in these latter two chapters is original, but the algebra is rather involved 
in places. The last half of the book does not, unfortunately, measure up to the standards of the first 
part either as to subject matter or presentation. As a whole, however, the book is well written and a very 


worthwhile contribution to the literature. . ats 
Gorpon F. NEWELL 
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Applied analysis. By Cornelius Lanczos, Prentice Hall, Inc., Englewood Cliffs, N. J., 
1956. xx + 539 pp. $9.00. 


Some ten years ago, satisfactory textbooks on numerical analysis in English could have been counted 
on the fingers of one hand. Since then, the spectacular development of high-speed computing equipment 
has revived the interest of mathematicians in this field. The reappearance of numerical analysis in the 
curricula of many colleges and universities was accompanied by the publication of several outstanding 
texts. These range in treatment from the strictly practical to the highly sophisticated, and in scope 
from collections of the authors’ favorite methods to comparative surveys of methods suitable for various 

pes of computing devices. In these terms of reference, the present book would have to be described 
as striking an excellent balance between the practical and the sophisticated points of view, but restrict- 
ing the discussion to a rather small number of methods. This statement, however, would fail to convey 
idea of the stimulating character of the book. As is indicated by their hyphenated names, many 


ati 
methods of numerical analysis have a long history of development by small improvements. If such a 
method were merely described and proved to furnish the desired result, the reader might be left wonder- 
ing how the complex procedure could ever have been invented. By concentrating on principal points 
and letting the reader fill in the details, the present author not only gains the space needed for a more 
adequate coverage of the development of the various methods, but actually makes the reader feel that 
he has participated in this development. 


Chapter I deals with algebraic equations, a ‘‘movable strip” version of Horner’s scheme and Bernoulli’s 
1ethod of moments being at the center of the discussion. Chapter II treats matrices and eigenvalue 


problems with weleome emphasis on geometrical aspects. Chapter III is devoted to the large-scale 
linear systems that result from the algebraisation of boundary value problems. In particular, the author’s 


spectroscopic method” of determining the real eigenvalues of large matrices and the corresponding 
method of solving large systems of linear equations are developed. Chapter IV is concerned with har- 
monic analysis, the emphasis being on the interpolation aspects of the Fourier series. Among the special 


tures of this chapter are discussions of the use of ‘‘e factors” for smoothing out the Gibbs oscillations, 


numerical aspects of the Fourier and Laplace transforms with applications to circuit design, and an 
improved method of searching for hidden periodicities. Under the title “Data Analysis’, Chapter V 
deals with interpolation and smoothing. Chapter VI treats quadrature methods, more than half of 
the chapter being devoted to Gaussian quadrature and several useful modifications of this. In particular, 


a quadrature formula is developed that employs only the values of the function and its derivatives up 
to a certain order at the two endpoints of the interval; the use of this formula in boundary value prob- 
lems is illustrated. The final chapter ‘‘Power Expansions” emphasizes the use of Chebyshev polynomials, 
iscussing the author’s ‘‘7-method” among other applications. 

W. PRAGER 


Momentum transfer in fluids. By William H. Corcoran, J. B. Opfell and B. H. Sage. 
Academic Press, New York, 1956. xi + 394 pp. $9.00. 


This is a book mainly written for the Chemical Engineer who does not have or need much previous 
vledge in the other branches of Fluid Mechanies. It consequently stresses the theoretical treatment 
the phenomena encountered in internal shear flows. In Chapter I the concept of momentum transfer 


is introduced and then applied to a number of cases involving laminar flow. It is concluded with dis- 
cussion of transition to turbulent flow and friction coefficients. Chapters II to IV deal with the properties, 
macroscopic characteristics and some analysis of turbulent flows, treating such topics as the similarity 

pothesis, the momentum transfer hypothesis, the vorticity transfer hypothesis and their applications 





to flows in channels. In Chapter V the basic equations of general Fluid Mechanics are developed. Chapter 
VI defines and discusses such terms as correlations, spectrum, etc. in the treatment of turbulence and 
the final chapter introduces the reader to the concept of the boundary layer and its various applications. 


While the main body of the book does not, make use of either vector or tensor analysis, an exhaustive 
and rigorous introduction to tensor analysis and the statistical theory of turbulence is added as an 


Appendix. It is the reviewers opinion that, while the book seems to be designed mainly for chemical 
ts chapters strike such a happy compromise between mathematical treatment and 


engineers, some ol 


physical discussion, as to warrant its use by other branches of engineering dealing with fluid mechanics. 
P. MAEDER 
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Elements of pure and applied mathematics. By Harry Lass. MeGraw-Hill Book Co., Inc., 
New York, Toronto, London, 1957. xi + 491 pp. $7.50. 


This is another addition to the growing collection of text books that give a very brief introduction 
to numerous topics. The chapter headings are: 1. linear equations, determinants, and matrices, 2. vector 
analysis, 3. tensor analysis, 4. complex-variable theory, 5. differential equations, 6. orthogonal poly- 
nomials, Fourier series and Fourier integrals, 7. the Stieltjies integral, Laplace transform, and calculus 
of variations, 8. group theory and algebraic equations, 9. probability theory and statistics, 10. real- 
variable theory. Each of these is about 40 pages except for the last chapter that is about 100 pages. 

As compared h most other books intended as text books for an advanced undergraduate or 
graduate level, this book covers more subjects in more condensed form and in a more abstract way 
although not in an extreme way. The book is partly “pure” mostly by virtue of the choice of topics; 





the presentation is not consistently more rigorous than other books. One annoying feature is the frequent 
ne chapters of mathematical concepts discussed in the tenth chapter. Aside 





use throughout the first 
from this, the book is fairly well organized and well written. It is perhaps just the type of book that 


some instructors would like to use, this being largely a matter of personal preference. 


G. F. NEWELL 


The the f hydrodynamic stability. By C. C. Lin. Cambridge University Press, New 


if rk, 1955. XI T 155 pp. $4.25 


This monograph presents a unified survey of the theory of stability with respect to infinitesimal 








distur I S of laminar fiows OI a homogen US viscous fluid. This is a subject to which the author has 
made many fundamental contributions and to which he has been a continual source of stimulation 
for furt by th | students and colleagues. In addition to the presentation of the general 
theor reful « { 1 of the mathematical questions raised by this theory there ar separate 
chapters on ( tte, I e Poiseuille, and Boundary Layer flows, and on the use of stability theory to 
problems in Astro- and Geophysics. Thx rious contributions of the many workers in the field are specifi- 
call l t I it ex! tive i ' , complete through 1955. 

\ : n perfor tion throughout the book of what 
parts of nsidered as known and what parts are controversial. This is particularly 
useful in vie i ther ysical and delicate mathematical questions that have been raised over the 
various u 1 mathemati t have been used in the development 
of the I should serve as I] ‘ce bor nd stimulus to anyone wish- 
ing to e! ictive nt esti gx field. The € re man sections which end with a ren suci is 
“‘adequat hods till lacking” or “‘a complete theory.....is as yet far from being completely 
develo 


PETER CHIARULLI 


General elat? wLy and ¢ »smology. By (. Cl. MeVittie. Joh Wiley « Son : Ine.. New 


York, 1956. x + 198 pp. $9.00. 


This is the fourth volume of the International Astrophy sics Serie 8, which aims to provi le authori- 
tative volumes, suitable for both specialists and students, dealing with the main branches of astrophysics 
and radio stronomy. On le afing through the book, two things catch the eye. First, a welcome rarity 
in technical literature, th ithor explains at some length his philosophical views about science S« ondly, 


he presents details which one is glad to have readily available. Thus, we find Dingle’s formulae for the 


Z 


Einstein tensor calculated for any orthogonal line element, and tables of observed results for the bending 
of a light ray grazing the sun and for red-shifts, solar and cosmological. The introductory paragraphs 
to the various chapters are excellent, and it was only when I started to study the mathematics of the 
book that I found myself out of tune with the author. 

Perhaps this is due to two related views I hold about the general theory of relativity: (1) it is a 
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thing of beauty, the essential simplicity of which must at all costs be made to shine through those 
formal manipulations which cannot be avoided in applying it, and (2) we must develop and use to the 


full our geometrical intuitions about space-time in the same way as mathematical physicists, through 


the ages, developed and used their kinematical intuitions about space and time in the Newtonian scheme. 
[ can best describe the book in a negative sense by saying that those views are not the author’s. He 
has no use for geometrical intuition in space-time, but, in supporting his own scepticism, he should 
not mislead the reader by telling him (p. 31) that in a Riemannian 2-space it is impossible to write 
down the curvature tensor, the Ricci tensor and the curvature invariant in terms of the Gaussian curva- 


ture; the very simple formulae which do this will be found on pp. 89, 96 of Synge and Schild, Tensor 
Calculus 


The following list of chapter headings will show the scope of the book: I. Introduction. II. The 
tensor calculus and Riemannian geometry. III. Newtonian mechanics and special relativity. IV. The 
principles of general relativity. V. The Schwarzschild space-time. VI. Approximations to Einstein’s 
equation and Newtonian gas-dynamics. VII. Special cases in Newtonian gas-dynamics. VIII. Theory 
of uniform model universes. IX. Model universes and the system of galaxies. There are references and 
notes at the end, and a good index. 

This is not, in my opinion, a book for a beginner in relativity. But for one who understands the theory 
already in its general outlines, and seeks additional special information, particularly in relation to 
cosmological problems, the last two chapters will provide much food for thought. 

Here are some minor criticisms. The derivation of the equations of a null geodesic (p. 22) seems 
insound, because it supposes falsely that all curves adjacent to a null geodesic are null to the first order. 
The view (p. 35) that the cosmological constant is essentially a constant of integration stems from the 
vuthor’s curious way of looking at the matter. In eq. (3.204) read mp for p, where m is molecular weight. 
[It is a pity (p. 51) to ery down Minkowski’s great achievement by referring to him as “‘one of the mathe- 
maticians who first employed” the idea of flat space-time; or is there another claimant to priority? 


J. L. SYNGE 


Introduction to statistical analysis. By Wilfrid J. Dixon and Frank J. Massey, Jr. MeGraw- 
Hill Book Co., Ine., New York, London, Toronto, 1957. xiii + 488 pp. $6.00. 


As the title might suggest, this is an undergraduate text for the non-specialist. The authors (who 
rk in biostatistics) have been careful to select varied material from industry, educational measure- 
nts, and commercial fields, as well as from their specialty; and every problem and exercise is in 








nguage comprehensible to the average undergraduate. So many topics are touched upon, and often 
dismissed with so few sentences, that the “introduction” is sometimes little more than a distant glimpse. 
There are few demands upon mathematical proficiency. In other respects the section headings pretty 
well cover the range of current statistical interest, save for the notable omission of the specialized topics 
relating to business forecasting and index numbers. Almost a quarter of the book consists of well-chosen 
tables. This second edition inclu des a brief terminal chapter on a posteriori probability. A useful feature 
of the form of presentation is the separate listing of ‘‘Discussion Questions,” “Class Exercises,” and 

Problems.” A chapter on “Microstatics” deals with the theory of small samples. Many exercises 
carry on in succession with data found by the individual student in sampling either by use of random- 
number tables or from populations of tags, disks, or beads. No new principles are developed in this text. 


ALBERT A, BENNETT 
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THIRD U. S. NATIONAL CONGRESS OF APPLIED MECHANICS 
BROWN UNIVERSITY, PROVIDENCE, RHODE ISLAND 
June 11-14, 1958 


All research workers in the field are cordially invited to submit papers constituting 
original experimental or theoretical contributions to Applied Mechanics, including 
mechanics of rigid bodies and deformable solids, mechanics of fluids and gases, thermo- 
dynamics and heat transfer. Instructions to authors of papers may be obtained from the 
Organizing Committee. It is expected that papers accepted by the Editorial Committee 
with the advice of recognized authorities and presented at the Congress, will be published 
in full in the Proceedings of the Congress. 

To be considered for presentation at the Congress, complete papers and manuscripts 
must be submitted to the Chairman of the Editorial Committee before January 1, 1958; 
to be scheduled for presentation the final manuscript of a paper must have been accepted 
before May 1, 1958. To avoid delays caused by overburdening reviewers and editorial 
staff, authors are urged to submit manuscripts well ahead of the deadline of January 
1, 1958. 


The papers will be grouped by subject and 30 minutes will be allotted for presentation 
and discussion of each paper. Arrangements will also be made for general lectures by 
outstanding authorities on subjects of general interest to members of the Congress. 
Facilities will be provided for informal discussions and social contact. 


THE ORGANIZING COMMITTEE 

W. Prager 
Secretary Daniel C. Drucker 
Treasurer... A eo ee ee oe sy sis E. H. Lee 


Chairman, Editorial Committee ....... . .R.M. Haythornthwaite 
Committee on Arrangements. ...... . .E. T. Onat and R. T. Shield 


Chairman 


INQUIRIES REGARDING THE CONGRESS SHOULD BE ADDRESSED TO THE SEC- 
RETARY OF THE CONGRESS AT BROWN UNIVERSITY, BOX D, PROVIDENCE 12, 
RHODE ISLAND. 





























